Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “characteristic equations”

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 487 records · Page 27

Magnetohydrodynamic MACH Code Used to Simulate Magnetoplasmadynamic Thrusters

The On-Board Propulsion program at the NASA Glenn Research Center is utilizing a state of-the-art numerical simulation to model the performance of high-power electromagnetic plasma thrusters. Such thrusters are envisioned for use in lunar and Mars cargo transport, piloted interplanetary expeditions, and deep-space robotic exploration of the solar system. The experimental portion of this program is described in reference 1. This article describes the numerical modeling program used to guide the experimental research. The synergistic use of numerical simulations and experimental research has spurred the rapid advancement of high-power thruster technologies for a variety of bold new NASA missions. From its inception as a U.S. Department of Defense code in the mid-1980's, the Multiblock Arbitrary Coordinate Hydromagnetic (MACH) simulation tool has been used by the plasma physics community to model a diverse range of plasma problems--including plasma opening switches, inertial confinement fusion concepts, compact toroid formation and acceleration, z-pinch implosion physics, laser-target interactions, and a variety of plasma thrusters. The MACH2 code used at Glenn is a time-dependent, two-dimensional, axisymmetric, multimaterial code with a multiblock structure. MACH3, a more recent three-dimensional version of the code, is currently undergoing beta tests. The MACH computational mesh moves in an arbitrary Lagrangian-Eulerian (ALE) fashion that allows the simulation of diffusive-dominated and dispersive-dominated problems, and the mesh can be refined via a variety of adaptive schemes to capture regions of varying characteristic scale. The mass continuity and momentum equations model a compressible viscous fluid, and three energy equations are used to simulate nonthermal equilibrium between electrons, ions, and the radiation field. Magnetic fields are modeled by an induction equation that includes resistive diffusion, the Hall effect, and a thermal source for magnetic fields. Various models of plasma resistivity are included, along with ablation models and multiport circuit solvers. The set of equations is closed using either an ideal gas or real equation of state.

Mikellides, Pavlos G.↗

Effect of Inductance and Requirements for Surge Current Testing of Tantalum Capacitors

Surge current testing is considered one of the most important techniques to evaluate reliability and/or screen out potentially defective tantalum capacitors for low-impedance applications. Analysis of this test, as it is described in the MIL-PRF-55365 document, shows that it does not address several issues that are important to assure adequate and reproducible testing. This work investigates the effect of inductance of the test circuit on voltage and current transients and analyzes requirements for the elements of the circuit, in particular, resistance of the circuit, inductance of wires and resistors, type of switching devices, and characteristics of energy storage bank capacitors. Simple equations to estimate maximum inductance of the circuit to prevent voltage overshooting and minimum duration of charging/discharging cycles to avoid decreasing of the effective voltage and overheating of the parts during surge current testing are suggested.

Teverosky, Alexander↗

A Piloted Simulator Evaluation of Transport Aircraft Rudder Pedal Force/Feel Characteristics

A piloted simulation study has been conducted in a fixed-base research simulator to assess the directional handling qualities for various rudder pedal feel characteristics for commercial transport airplanes. That is, the effects of static pedal force at maximum pedal travel, breakout force, and maximum pedal travel on handling qualities were studied. An artificial maneuver with a severe lateral wind shear and requiring runway tracking at an altitude of 50 feet in a crosswind was used to fully exercise the rudder pedals. Twelve active airline pilots voluntarily participated in the study and flew approximately 500 maneuvers. The pilots rated the maneuver performance with various rudder pedal feel characteristics using the Cooper- Harper rating scale. The test matrix had 15 unique combinations of the 3 static pedal feel characteristics. A 10-term, second-order equation for the Cooper-Harper pilot rating as a function of the 3 independent pedal feel parameters was fit to the data. The test matrix utilized a Central Composite Design that is very efficient for fitting an equation of this form. The equation was used to produce contour plots of constant pilot ratings as a function of two of the parameters with the third parameter held constant. These contour plots showed regions of good handling qualities as well as regions of degraded handling qualities. In addition, a numerical equation solver was used to predict the optimum parameter values (those with the lowest pilot rating). Quantitative pilot performance data were also analyzed. This analysis found that the peak values of the cross power spectra of the pedal force and heading angle could be used to quantify the tendency toward directional pilot induced oscillations (PIO). Larger peak values of the cross power spectra were correlated with larger (degraded) Cooper-Harper pilot ratings. Thus, the subjective data (Cooper-Harper pilot ratings) were consistent with the objective data (peak values of the cross power spectra).

Stewart, Eric C.↗

Assessing the Age of an Asteroid's Surface with Data from the International Rosetta Mission

Rosetta is an international mission led by the European Space Agency (ESA) with key support and instrumentation from the National Aeronautics and Space Administration (NASA). Rosetta is currently on a ten-year mission to catch comet 67P/Churyumov-Gerasimenko (C-G); throughout its voyage, the spacecraft has performed flybys of two main belt asteroids (MBA): Steins and Lutetia. Data on the physical, chemical, and geological properties of these asteroids are currently being processed and analyzed. Accurate interpretation of such data is fundamental in the success of Rosetta's mission and overall objectives. Post-flyby data analyses strive to correlate the size, shape, volume, and rotational rate of Lutetia, in addition to interpreting its multi-color imagining, albedo, and spectral mapping. Although advancements in science have contributed to the examination of celestial bodies, methods to analyze asteroids remain largely empirical, not semi-empirical, nor ab initio. This study aims to interpret and document the scientific methods currently utilized in the characterization of asteroid (21) Lutetia in order to render these processes and methods accessible to the public. Examples include a standardized technique for assessing the age of an asteroid surface, complete with clickable reference maps, methodology of grouping surface characteristics together, and a standardized power law equation for the age. Other examples include determining the density of an object. Context for what both density and age mean is a bi-product of this study. Results of the study will aid in the development of pedagogical material on asteroids for public use, and in creation of an academic database for selected targets that might be used as a reference.

asteriods↗

A Transient Hydrodynamic Model of Screen Channel Liquid Acquisition Devices for In-Space Cryogenic Propellant Management

Screen channel liquid acquisition devices (LADs) will play a crucial role in future deep space travel. It is essential that vapor-free delivery of propellants during tank-to-tank transfer is ensured to maximize yield from storage tanks and prevent potential combustion instabilities. The screen channel LAD utilizes a fine screen wire mesh that can separate phases in a low Bond number (i.e. microgravity) environment using surface tension forces. This study presents the development and verification of a new model for transient screen compliance, one of the influential factors for screen channel LAD design. Screen compliance is crucial during LAD channel outflow transients because the slight deflection of the screen can provide needed mass to satisfy rapid outflow demands and reduce the pressure difference across the screen. The model is successfully verified against CFD simulations. In addition, the characteristic speed for the governing screen compliance equations is derived which allows for numerical stability criteria to be established. As shown in this study, the transient maximum pressure difference across the screen can greatly exceed the steady state maximum pressure difference across the screen in many cases.

Hydrodynamics Simulations↗

Design of a 3rd Order Lumped-Element Filtering Attenuator

This paper presents a novel filtering attenuator (filtenuator) concept, which allows for adjustable passband attenuation and filtering characteristics in a single design. The mathematical equations relating the desired passband loss to required resistance values are provided and validated with ideal circuit simulations. A third-order Chebyshev static filtenuator is designed at 1 GHz with multiple values of attenuation using lumped-element components as a proof of concept. Simulated and measured results for a 3 dB filtenuator show good agreement and several attenuation levels are shown to validate the concept. The proposed static filtenuator has advantages of ease of modification, compact size, and low cost, which are necessities in meeting the needs of reducing size, weight, power, and cost (SWaP-C) in future radio frequency transceivers.

Filter, Attenuator, Band-pass Filter (BPF), Lumped↗

Analysis of some fully-discrete algorithms for the one-dimensional heat equation

The present investigation is concerned with a fully discrete accuracy and stability analysis of the one-dimensional heat equation, taking into account the evaluation of two-pass explicit schemes which simultaneously employ lumped and coupled capacity matrices. Schemes of the considered characteristics are not amenable to uncoupled semidiscrete and ordinary differential equation analyses. The obtained results illustrate that superior behavior may be achieved by schemes of the employed type when compared with the performance of the standard one-pass explicit schemes. The key idea in the considered approach is related to the utilization of a reduced-quadrature capacity matrix in the evaluation of the right-hand-side residual.

Hughes, T. J. R.↗

Finite-element solution for a combined radiative-conductive analysis with mixed diffuse-specular surface characteristics

The finite-element method has been applied to solve a combined radiative-conductive heat transfer problem for a large space telescope similar to those used in orbiting satellites. The derivation of the underlying matrices and associated solution algorithm for a 2-D triangular element is presented in detail. The resulting expressions for this triangular element typify such an analysis, which yields constitutive matrices when the heat equation is cast in the matrix form. The relevant matrices include those pertaining to thermal conductance, internal heat generation, radiative exchanges, and all possible external thermal loadings. Emphasis is placed on the treatment of non-linear radiative interchange between surfaces in an enclosure having mixed diffuse-specular surface characteristics. Essential differences in governing equations describing these distinctive surface characteristics are identified. Concluding remarks are drawn from an example simulating a Cassegrainian space telescope.

Lee, H.-P.↗

The characteristic Rossby frequency

The characteristic Rossby frequency is defined for a fixed zonal wavenumber perturbation as the variational integral of the Rayleigh-Ritz method. It is a measure of the time scale of the disturbance. For a disturbance which locally has the shape of an eigenfunction but is not global in extent, the characteristic Rossby frequency is very close to the true eigenvalue, and additionally remains unchanged under linear inviscid dynamics. Results are presented for the shallow water equations, both with and without a mean zonal wind. The characteristic Rossby frequency of a wavenumber 1 perturbation having the shape of the second symmetric Rossby mode but confined to the Northern Hemisphere is close to the corresponding Rossby frequency. This finding is helpful in understanding the behavior of the observed wavenumber 1 pattern of January 1979, which propagated westward with nearly the pure Rossby frequency but was discernible only in the Northern Hemisphere (as discussed by Daley and Williamson).

Straus, David M.↗

Parallel-in-Time Solution of Hyperbolic PDE Systems via Characteristic-Variable Block Preconditioning

We consider the parallel-in-time solution of both linear and nonlinear hyperbolic partial differential equation (PDE) systems in one spatial dimension. In the nonlinear setting, the discretized equations are solved with a preconditioned residual iteration based on a global linearization. The linear(ized) equation systems are approximately solved parallel-in-time using a block preconditioner applied in the characteristic variables of the underlying linear(ized) hyperbolic PDE. This change of variables is motivated by the observation that intervariable coupling between characteristic variables is weak, at least locally where spatio-temporal variations in the eigenvectors of the associated flux Jacobian are sufficiently small, while that between the original variables is not. For an ℓ-dimensional system of PDEs, applying the preconditioner consists of solving a sequence of ℓ scalar linear(ized)-advection-like problems, each associated with a different characteristic wave-speed in the underlying linear(ized) PDE. Furthermore, we approximately solve these linear advection problems using multigrid reduction-in-time (MGRIT); however, any other suitable parallel-in-time method could be used. Numerical examples are shown for the (linear) acoustics equations in heterogeneous media and for the (nonlinear) shallow water equations and Euler equations of gas dynamics with shocks and rarefactions. For many test problems, the solver converges in just a handful of iterations and with mesh-independent convergence rates.

97 MATHEMATICS AND COMPUTING↗

Detection of Moving Targets Using Soliton Resonance Effect

The objective of this research was to develop a fundamentally new method for detecting hidden moving targets within noisy and cluttered data-streams using a novel "soliton resonance" effect in nonlinear dynamical systems. The technique uses an inhomogeneous Korteweg de Vries (KdV) equation containing moving-target information. Solution of the KdV equation will describe a soliton propagating with the same kinematic characteristics as the target. The approach uses the time-dependent data stream obtained with a sensor in form of the "forcing function," which is incorporated in an inhomogeneous KdV equation. When a hidden moving target (which in many ways resembles a soliton) encounters the natural "probe" soliton solution of the KdV equation, a strong resonance phenomenon results that makes the location and motion of the target apparent. Soliton resonance method will amplify the moving target signal, suppressing the noise. The method will be a very effective tool for locating and identifying diverse, highly dynamic targets with ill-defined characteristics in a noisy environment. The soliton resonance method for the detection of moving targets was developed in one and two dimensions. Computer simulations proved that the method could be used for detection of singe point-like targets moving with constant velocities and accelerations in 1D and along straight lines or curved trajectories in 2D. The method also allows estimation of the kinematic characteristics of moving targets, and reconstruction of target trajectories in 2D. The method could be very effective for target detection in the presence of clutter and for the case of target obscurations.

Kulikov, Igor K.↗

Compressive behavior of unidirectional fibrous composites

The longitudinal compressive behavior of unidirectional fiber composites was investigated by using the Illinois Institute of Technology Research Institute (IITRI) test method with thick and thin test specimens. The test data obtained are interpreted by means of stress/strain curves from back-to-back strain gages, examination of fracture surfaces by scanning electron microscope, and predictive equations for distinct failure modes including fiber compression failure. Euler buckling, delamination, and flexure. The results show that longitudinal compressive fracture is induced by a combination of delamination, flexure, and fiber tier breaks. No distinct fracture surface characteristics can be associated with unique failure modes. An equation is described that can be used to extract the longitudinal compressive strength from the longitudinal tensile and flexural strengths of the same composite system.

Sinclair, J. H.↗

Characteristics-based methods applied to infinite Prandtl number thermal convection in the hard turbulent regime

Characteristics-based methods for the advection-diffusion equation are presented and directly applied to study thermal convection with extremely large Rayleigh number (Ra). It is shown that the operator-splitting method for advection-diffusion problems is very accurate for determining the advected field at extremely high Peclet number (Pe). The technique presented is considered to have great potential for solving advection-dominated problems, while the Langrangian method is more accurate for lower Pe. It is noted that the accuracy of these characteristics-based methods strongly depends on the quality of interpolation. The computational time for the operator-splitting method grows with the number of time steps employed. The Langrangian method was used for simulations of convection at very high Ra, up to 3 x 10 to the 9th, and time-dependent, thermal convection solutions were obtained for infinite Prandtl number.

Malevsky, A. V.↗

Compression behavior of unidirectional fibrous composite

The longitudinal compression behavior of unidirectional fiber composites is investigated using a modified Celanese test method with thick and thin test specimens. The test data obtained are interpreted using the stress/strain curves from back-to-back strain gages, examination of fracture surfaces by scanning electron microscope, and predictive equations for distinct failure modes including fiber compression failure, Euler buckling, delamination, and flexure. The results show that the longitudinal compression fracture is induced by a combination of delamination, flexure, and fiber tier breaks. No distinct fracture surface characteristics can be associated with unique failure modes. An equation is described which can be used to extract the longitudinal compression strength knowing the longitudinal tensile and flexural strengths of the same composite system.

Sinclair, J. H.↗

Flapping response characteristics of hingeless rotor blades by a gereralized harmonic balance method

Linearized equations of motion for the flapping response of flexible rotor blades in forward flight are derived in terms of generalized coordinates. The equations are solved using a matrix form of the method of linear harmonic balance, yielding response derivatives for each harmonic of the blade deformations and of the hub forces and moments. Numerical results and approximate closed-form expressions for rotor derivatives are used to illustrate the relationships between rotor parameters, modeling assumptions, and rotor response characteristics. Finally, basic hingeless rotor response derivatives are presented in tabular and graphical form for a wide range of configuration parameters and operating conditions.

Peters, D. A.↗

Computer program for investigating effects of nonlinear suspension-system elastic properties on parachute inflation loads and motions

A computer program is presented by which the effects of nonlinear suspension-system elastic characteristics on parachute inflation loads and motions can be investigated. A mathematical elastic model of suspension-system geometry is coupled to the planar equations of motion of a general vehicle and canopy. Canopy geometry and aerodynamic drag characteristics and suspension-system elastic properties are tabular inputs. The equations of motion are numerically integrated by use of an equivalent fifth-order Runge-Kutta technique.

Poole, L. R.↗

Deconvolution Methods and Systems for the Mapping of Acoustic Sources from Phased Microphone Arrays

Mapping coherent/incoherent acoustic sources as determined from a phased microphone array. A linear configuration of equations and unknowns are formed by accounting for a reciprocal influence of one or more cross-beamforming characteristics thereof at varying grid locations among the plurality of grid locations. An equation derived from the linear configuration of equations and unknowns can then be iteratively determined. The equation can be attained by the solution requirement of a constraint equivalent to the physical assumption that the coherent sources have only in phase coherence. The size of the problem may then be reduced using zoning methods. An optimized noise source distribution is then generated over an identified aeroacoustic source region associated with a phased microphone array (microphones arranged in an optimized grid pattern including a plurality of grid locations) in order to compile an output presentation thereof, thereby removing beamforming characteristics from the resulting output presentation.

Brooks, Thomas F.↗

Numerical solution of the two-dimensional Euler equations by second-order upwind difference schemes

Two time-level, five-point explicit and implicit upwind difference schemes based on the characteristic flux difference splitting concept have been developed for the two-dimensional Euler equations. The method is conservative, second-order accurate in time and space, and general coordinate systems are used to treat complex geometries. Nonlinear flux limiters are employed to yield oscillation free sharp shock profiles. Upstream interpolation is used to yield a class of higher-order upwind schemes which closely mimic the locally one-dimensional method of characteristics (with fixed time intervals) through operator splitting. Numerical results have been obtained for a plane shock reflection and for flow over a circular arc in a channel. Characteristics of upwind TVD schemes, as applied to two-dimensional flows with embedded shocks are discussed.

Yang, J. Y.↗