SEARCH · Engineering Papers

Results for “numerical approximation & analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 361 records · Page 20

All sources NASA NTRS DOE OSTI

NASA NTRS1983

Voyager 2 observations of energetic particle variations in the Ganymede wake region - A possible acceleration mechanism

Voyager 2's passage through the downstream corotation region of Ganymede found disturbances in the field and particle environment. Large fluctuations in the intensities and energy spectra of ions in the 0.1-to 4.0-MeV interval were also observed with the low energy charged particle experiment. All ion species were evidently affected, up through medium (C, N, O) nuclei. Electrons intensities, while disturbed, did not behave similarly to the ions. These effects are probably associated with the wake produced by the absorption of Jovian magnetospheric corotating plasma by Ganymede. A two-dimensional theoretical model for perturbed field in the wake region is proposed and it is shown by direct numerical simulation of exact particle trajectories that such perturbations are capable of modulating the charged particle energies to about the observed amount for the observed magnetic and plasma parameters. Analysis shows that the energization process is dependent on the pitch angles of the energetic ions and the process itself is nonadiabatic (the Larmor radius is approximately the size of the disturbance region).

Tariq, G.F.↗

NASA NTRS1988

Interaction and merging of vortex filaments

The asymptotic solutions of Navier-Stokes equations for vortex filaments of finite strength with small effective vortical cores are summarized with special emphasis placed on the physical meaning and the practical limit to the applicability of the asymptotic solution. Finite-difference solutions of Navier-Stokes equations for the marging of the filament(s) are described with a focus on the development of the approximate boundary conditions for the computational domain. An efficiency study employing a model problem is used to assess the advantages of the present approximate boundary condition method over previously used techniques. Applications of the present method are presented for the motion and decay of a 3:1 elliptic vortex ring, and for the merging process of a pair of coaxial vortex rings. A numerical procedure for the problem of local merging of vortex filaments, which requires the asymptotic analysis as well as the numerical Navier-Stokes solver, is also presented.

Liu, C. H.↗

NASA NTRS1944

The flow of a compressible fluid past a circular arc profile

The Ackeret iteration process is utilized to obtain higher approximations than that of Prandtl and Glauert for the flow of a compressible fluid past a circular arc profile. The procedure is to expand the velocity potential in a power series of the camber coefficient. The first two terms of the development correspond to the Prandtl-Glauert approximation and yield the well-known correction to the circulation about the profile. The second approximation, involving the square of the camber coefficient, improves the velocity and pressure fields but yields no new results with regard to the circulation, since the circulation about the profile is an odd function of the camber coefficient. The third approximation, involving the cube of the camber coefficient, permits the use of higher values of the camber coefficient and furthermore yields an improvement to the Prandtl-Glauert rule with regard to the effect of compressibility on the circulation of the circular arc profile. Numerical examples with tables and graphs illustrate the results of the analysis.

Kaplan, Carl↗

NASA NTRS1994

Gradient drift eigenmodes in the equatorial electrojet

The problem of kilometer-scale irregularities in the daytime equatorial electrojet is revisited by means of an eigenmode analysis of the gradient drift instability. Realistic physical parameters are used, including the modeled altitude variations of ion and electron collision frequencies and mobilities. The full fourth-order system of two coupled differential equations (each of second order) for the denisty and electrostatic potential perturbations is solved numerically by a relaxation technique. Under some approximations, the fourth-order system can be shown to reduce to a second-order differential equation for the perturbed potential or density. The latter is solved using a shooting technique and provides initial guesses for numerical solutions to the full problem. It is shown that the linear growth rate peaks for kilometer-scale waves, contrary to the findings of recent initial-value studies. This occurs because the equilibrium velocity shear is much more effective as a damping mechanism for short-wavelength modes than it is for the longer, kilometer-scale modes. These results provide a natural qualitative explanation for the observed dominance of kilometer-scale structures in the daytime electrojet spectrum.

Wang, X.-H.↗

NASA NTRS2024

Adiabatic Gravitational Waveform Model for Compact Objects Undergoing Quasi-Circular Inspirals Into Rotating Massive Black Holes

We present bhpwave: a new python-based, open-source tool for generating the gravitational waveforms of stellar-mass compact objects undergoing quasicircular inspirals into rotating massive black holes. These binaries, known as extreme-mass-ratio inspirals (EMRIs), are exciting mHz gravitational wave sources for future space-based detectors such as the Laser Interferometer Space Antenna (LISA). Relativistic models of EMRI gravitational wave signals are necessary to unlock the full scientific potential of mHz detectors, yet few open-source EMRI waveform models exist. Thus we built bhpwave, which uses the adiabatic approximation from black hole perturbation theory to rapidly construct gravitational waveforms based on the leading-order inspiral dynamics of the binary. In this work, we present the theoretical and numerical foundations underpinning bhpwave. We also demonstrate how bhpwave can be used to assess the impact of EMRI modeling errors on LISA gravitational wave data analysis. In particular, we find that for retrograde orbits and slowly spinning black holes we can mismodel the gravitational wave phasing by as much as ∼10 radians without significantly biasing EMRI parameter estimation.

Zachary Nasipak↗

NASA NTRS2008

Evolution of the Mobile Information SysTem (MIST)

The Mobile Information SysTem (MIST) had its origins in the need to determine whether commercial off the shelf (COTS) technologies could improve intervehicular activities (IVA) on International Space Station (ISS) crew maintenance productivity. It began with an exploration of head mounted displays (HMDs), but quickly evolved to include voice recognition, mobile personal computing, and data collection. The unique characteristic of the MIST lies within its mobility, in which a vest is worn that contains a mini-computer and supporting equipment, and a headband with attachments for a HMD, lipstick camera, and microphone. Data is then captured directly by the computer running Morae(TM) or similar software for analysis. To date, the MIST system has been tested in numerous environments such as two parabolic flights on NASA's C-9 microgravity aircraft and several mockup facilities ranging from ISS to the Altair Lunar Sortie Lander. Functional capabilities have included its lightweight and compact design, commonality across systems and environments, and usefulness in remote collaboration. Human Factors evaluations of the system have proven the MIST's ability to be worn for long durations of time (approximately four continuous hours) with no adverse physical deficits, moderate operator compensation, and low workload being reported as measured by Corlett Bishop Discomfort Scale, Cooper-Harper Ratings, and the NASA Total Workload Index (TLX), respectively. Additionally, through development of the system, it has spawned several new applications useful in research. For example, by only employing the lipstick camera, microphone, and a compact digital video recorder (DVR), we created a portable, lightweight data collection device. Video is recorded from the participants point of view (POV) through the use of the camera mounted on the side of the head. Both the video and audio is recorded directly into the DVR located on a belt around the waist. This data is then transferred to another computer for video editing and analysis. Another application has been discovered using simulated flight, in which, a kneeboard is replaced with mini-computer and the HMD to project flight paths and glide slopes for lunar ascent. As technologies evolve, so will the system and its application for research and space system operations.

Litaker, Harry L., Jr.↗

NASA NTRS1984

MMIC devices for active phased array antennas

The study of printed circuit discontinuities is necessary in order to design, for example, transitions between rectangular waveguides and printed circuits. New developments with respect to the analytical approaches to this problem are discussed. A summary of the progress in the experimental approach is presented. The accurate solution for the modes in various millimeter-wave waveguides is essential in the analysis of many integrated circuit components, such as filters and impedance transformers. Problems associated with the numerical computation of these modes in two frequently used waveguide forms, namely, the finline and microstrip, are presented. The spectral domain method of formulation, with a moment method solution, is considered. This approach can be readily extended to analyze an arbitrary configuration of dielectric and metallized regions in a shielded enclosure. Galerkin's method is used, where the testing and basic functions are the same. It is shown that the mode functions, or eigenfunctions, are more sensitive to errors than the phase constants, or eigenvalues. The approximate mode functions do not satisfy the orthogonality relationship well, resulting in difficulties when these modal solutions are used to form an approximate Green's function or are used in a mode matching analysis.

Mittra, R.↗

NASA NTRS1981

Numerical analysis of spectral properties of coupled oscillator Schroedinger operators. I - Single and double well anharmonic oscillators

Several methods for computing many eigenvalues and eigenfunctions of a single anharmonic oscillator Schroedinger operator whose potential may have one or two minima are described. One of the methods requires the solution of an ill-conditioned generalized eigenvalue problem. This method has the virtue of using a bounded amount of work to achieve a given accuracy in both the single and double well regions. Rigorous bounds are given, and it is proved that the approximations converge faster than any inverse power of the size of the matrices needed to compute them. The results of computations for the g:phi(4):1 theory are presented. These results indicate that the methods actually converge exponentially fast.

Isaacson, D.↗

NASA NTRS2022

Nonlinear Evolution of Instabilities in a Laminar Separation Bubble at a Hypersonic Mach Number

The development of both convective stationary perturbation as well as global instabilities in the vicinity of a laminar separation bubble above an axisymmetric compression corner in a hypersonic flow is investigated using numerical simulations. The flow configuration of primary interest corresponds to the cone-cylinder-flare model used in experimental measurements in the Boeing/AFOSR Mach-6 Quiet Tunnel at Purdue University. For a flare angle of 10 degrees and a unit Reynolds number of 11.5 x 10 6 m -1 , their surface flow visualizations identified the presence of streamwise elongated thermal streaks near the reattachment position and the dominant azimuthal spacing between the streaks was determined to be approximately 10 degrees (i.e., an azimuthal mode number of 36). Previous linear stability analyses predicted that the amplification characteristics of small amplitude, unsteady, convective instabilities within this flow were consistent with the surface pressure fluctuations measured in the experiment. However, their accompanying analysis of global instabilities had found the separation bubble to be weakly unstable at the 10 degrees flare angle, with the most unstable global mode corresponding to a stationary disturbance with an azimuthal wavenumber of m ≈ 5, which was well below the measured wavenumber of m = 36. Besides confirming the presence of the global instability at these flow conditions, the present numerical simulations quantify the details of the stationary equilibrium state associated with the supercritical bifurcation resulting from the nonlinear saturation of the unstable global mode. Although velocity perturbations associated with the saturated global mode are dominated by the fundamental spanwise wavelength associated with the linear global instability, the surface heat flux downstream of reattachment is dominated by m = 36 in agreement with experimental measurements. However, the shorter wavelength peak ( m = 36) in the azimuthal spectrum is not manifested as streaks with an azimuthal spacing of 10 degrees. An additional simulation of the combined evolution of both unsteady convective and global instabilities indicates that the nonlinear interactions between a broadband spectrum of oblique Mack’s first mode instabilities generate stationary vortex modes with a higher dominant wavenumber ( m = 60) as that observed in the experiments. Additional analyses are performed to predict the global mode behavior at other flare angles and how it differs from the convective centrifugal instability due to concave streamline curvature over the compression corner.

Hypersonic flow↗

NASA NTRS2022

Nonlinear Evolution of Instability in a Laminar Separation Bubble at Hypersonic Mach Number

High speed flow↗

NASA NTRS1986

Stability of semidiscrete approximations for hyperbolic initial-boundary-value problems: An eigenvalue analysis

A hyperbolic initial-boundary-value problem can be approximated by a system of ordinary differential equations (ODEs) by replacing the spatial derivatives by finite-difference approximations. The resulting system of ODEs is called a semidiscrete approximation. A complication is the fact that more boundary conditions are required for the spatially discrete approximation than are specified for the partial differential equation. Consequently, additional numerical boundary conditions are required and improper treatment of these additional conditions can lead to instability. For a linear initial-boundary-value problem (IBVP) with homogeneous analytical boundary conditions, the semidiscrete approximation results in a system of ODEs of the form du/dt = Au whose solution can be written as u(t) = exp(At)u(O). Lax-Richtmyer stability requires that the matrix norm of exp(At) be uniformly bounded for O less than or = t less than or = T independent of the spatial mesh size. Although the classical Lax-Richtmyer stability definition involves a conventional vector norm, there is no known algebraic test for the uniform boundedness of the matrix norm of exp(At) for hyperbolic IBVPs. An alternative but more complicated stability definition is used in the theory developed by Gustafsson, Kreiss, and Sundstrom (GKS). The two methods are compared.

Warming, Robert F.↗

NASA NTRS1976

Heat transfer evaluation in a plasma core reactor

Numerical evaluations of heat transfer in a fissioning uranium plasma core reactor cavity, operating with seeded hydrogen propellant, was performed. A two-dimensional analysis is based on an assumed flow pattern and cavity wall heat exchange rate. Various iterative schemes were required by the nature of the radiative field and by the solid seed vaporization. Approximate formulations of the radiative heat flux are generally used, due to the complexity of the solution of a rigorously formulated problem. The present work analyzes the sensitivity of the results with respect to approximations of the radiative field, geometry, seed vaporization coefficients and flow pattern. The results present temperature, heat flux, density and optical depth distributions in the reactor cavity, acceptable simplifying assumptions, and iterative schemes. The present calculations, performed in cartesian and spherical coordinates, are applicable to any most general heat transfer problem.

Smith, D. E.↗

NASA NTRS1979

Stability of sequences generated by nonlinear differential systems

A local stability analysis is presented for both the analytic and numerical solutions of the initial value problem for a system of ordinary differential equations. It is shown that, using a proper choice of Liapunov function, a connected region of stable initial values of both the analytic solution and the one-leg k-step numerical solution can be approximated. Attention is given to the example of the two-dimensional problem involving the stability of the longitudinal equations of motion of a gliding jet aircraft.

Brown, R. L.↗

NASA NTRS1987

Factorization and the synthesis of optimal feedback kernels for differential-delay systems

A combination of ideas from the theories of operator Riccati equations and Volterra factorizations leads to the derivation of a novel, relatively simple set of hyperbolic equations which characterize the optimal feedback kernel for the finite-time regulator problem for autonomous differential-delay systems. Analysis of these equations elucidates the underlying structure of the feedback kernel and leads to the development of fast and accurate numerical methods for its computation. Unlike traditional formulations based on the operator Riccati equation, the gain is characterized by means of classical solutions of the derived set of equations. This leads to the development of approximation schemes which are analogous to what has been accomplished for systems of ordinary differential equations with given initial conditions.

Milman, Mark M.↗

NASA NTRS1994

Anisotropic Alfven-ballooning modes in Earth's magnetosphere

We have carried out a theoretical analysis of the stability and parallel structure of coupled shear Alfven and slow magnetosonic waves in Earth's inner magnetopause (i.e., at equatorial distances between about five and ten Earth radii) including effects of finite anisotropic Grad-Shafranov equation yields an approximate self-consistent magnetohydrodynamic (MHD) equilibrium. This MHD equilibrium is used in the numerical solution of a set of eigenmode equations which describe the field line eigenfrequency, linear stability, and parallel eigenmode structure. We call these modes anisotropic Alfven-ballooning modes. The main results are: (1) The field line eigenfrequency can be significantly lowered by finite pressure effects. (2) The parallel mode structure of the transverse wave components is fairly insensitive to changes in the plasma pressure, but the compressional magnetic component can become highly peaked near the magnetic equator as a result of increased pressure, especially when P(sub perpendicular to) is greater than P(sub parallel) (here P(sub perpendicular to) and P(sub parallel) are the perpendicular and parallel plasma pressure). (3) For the isotropic (P(sub parallel) = P(sub perpendicular to) = P) case ballooning instability can occur when the ratio of the plasma presure to the magnetic pressure, beta = P/(B squared/8 pi), exceeds a critical value beta(sup B)(sub 0) is approximately equal to 3.5 at the equator. (4) Compared to the isotropic case the critical beta value is lowered by anisotropy, either due to decreased field line bending stabilization when P(sub parallel) is greater than P(sub perpendicular to) or due to increased ballooning-mirror destabilization when P(sub perpendicular to) is greater than P(sub parallel). (5) We use a beta-delta stability diagram to display the regions of instability with respect to the equatorial values of the parameters bar beta and delta, where bar beta = (1/3)(beta(sub parallel) + 2 beta(sub perpendicular to)) is an average beta value and delta = 1 - P(sub parallel)/P(sub perpendicular to) is a measure of the plasma anisotropy. The diagram is divided into regions corresponding to the firehose, mirror and ballooning instabilities. It appears that observed values of the plasma pressure are below the critical value for the isotropic ballooning instability but it may be possible to approach a ballooning-mirror instability when P(sub perpendicular to)/P(sub parallel) is greater than or approximately 2.

Chan, Anthony A.↗

NASA NTRS1995

An overview of engineering concepts and current design algorithms for probabilistic structural analysis

The article begins by examining the fundamentals of traditional deterministic design philosophy. The initial section outlines the concepts of failure criteria and limit state functions two traditional notions that are embedded in deterministic design philosophy. This is followed by a discussion regarding safety factors (a possible limit state function) and the common utilization of statistical concepts in deterministic engineering design approaches. Next the fundamental aspects of a probabilistic failure analysis are explored and it is shown that deterministic design concepts mentioned in the initial portion of the article are embedded in probabilistic design methods. For components fabricated from ceramic materials (and other similarly brittle materials) the probabilistic design approach yields the widely used Weibull analysis after suitable assumptions are incorporated. The authors point out that Weibull analysis provides the rare instance where closed form solutions are available for a probabilistic failure analysis. Since numerical methods are usually required to evaluate component reliabilities, a section on Monte Carlo methods is included to introduce the concept. The article concludes with a presentation of the technical aspects that support the numerical method known as fast probability integration (FPI). This includes a discussion of the Hasofer-Lind and Rackwitz-Fiessler approximations.

Duffy, S. F.↗

NASA NTRS1988

Transonic aeroelastic analysis of launch vehicle configurations

A numerical study of the aeroelastic stability of typical launch vehicle configurations in transonic flight is performed. Recent computational fluid dynamics techniques are used to simulate the transonic aerodynamic flow fields, as opposed to relying on experimental data for the unsteady aerodynamic pressures. The flow solver is coupled to an appropriate structural representation of the vehicle. The aerodynamic formulation is based on the thin layer approximation to the Reynolds-Averaged Navier-Stokes equations, where the account for turbulent mixing is done by the two-layer Baldwin and Lomax algebraic eddy viscosity model. The structural-dynamic equations are developed considering free-free flexural vibration of an elongated beam with variable properties and are cast in modal form. Aeroelastic analyses are performed by integrating simultaneously in the two sets of equations. By tracing the growth or decay of a perturbed oscillation, the aeroelastic stability of a given constant configuration can be ascertained. The method is described in detail, and results that indicate its application are presented. Applications include some validation cases for the algorithm developed, as well as the study of configurations known to have presented flutter programs in the past.

Filgueirasdeazevedo, Joao Luiz↗

NASA NTRS2002

Transport of Space Environment Electrons: A Simplified Rapid-Analysis Computational Procedure

A computational procedure for describing transport of electrons in condensed media has been formulated for application to effects and exposures from spectral distributions typical of electrons trapped in planetary magnetic fields. The procedure is based on earlier parameterizations established from numerous electron beam experiments. New parameterizations have been derived that logically extend the domain of application to low molecular weight (high hydrogen content) materials and higher energies (approximately 50 MeV). The production and transport of high energy photons (bremsstrahlung) generated in the electron transport processes have also been modeled using tabulated values of photon production cross sections. A primary purpose for developing the procedure has been to provide a means for rapidly performing numerous repetitive calculations essential for electron radiation exposure assessments for complex space structures. Several favorable comparisons have been made with previous calculations for typical space environment spectra, which have indicated that accuracy has not been substantially compromised at the expense of computational speed.

Nealy, John E.↗