Engineering Papers⌕ Search

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.

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

At least 451 records · Page 25
All sourcesNASA NTRSDOE OSTI
NASA NTRS1976

Estimates of nitric oxide production for lifting spacecraft reentry

An approximate analysis to estimate the quantity of nitric oxide that could be formed in the wake of a reentering lifting spacecraft is reported. Three different approaches are undertaken: two simplified analytical models utilizing the sudden-freezing concept named the 'trailing-edge-freezing' and 'wake-freezing' approximations, and a computer calculation involving numerical integration of chemical rate equations. The three methods predict a maximum nitric oxide production equivalent to 6, 9.5 and 8% of the mass of the spacecraft, respectively. Thus, the amount of nitric oxide expected to be produced by future space activities is negligibly small compared with that produced by the natural processes.

Park, C.↗
NASA NTRS1983

Helicopter rotor loads using discretized matched asymptotic expansions

The numerical practicality of a matched asymptotic expansion approach for the computation of unsteady three dimensional airloads on a helicopter rotor was improved. This effort utilizes a discretized repesentation of the doublet strength distribution and helical streamlines to decrease the computational requirements of the original analysis. The continuous variation of the doublet strength was approximated by piecewise constant or piecewise quadratic distributions, and the helical trajectory of a fluid particle was approximated by connected straight line segments. As a direct result of these simplified representations the computational time required for the execution of a typical flight condition was reduced by an order of magnitude with respect to the requirements of the original analysis. Airloads which were computed using the discretized method for a two bladed model rotor and a full scale four bladed rotor are in close agreement with measured results and airloads from the original asymptotic analysis. For conditions characterized by significant rotor/wake interaction the piecewise constant representation requires a reduced azimuth spacing to maintain acceptable accuracy.

Pierce, G. A.↗
NASA NTRS1990

Effect of interleaves on fracture of laminated composites. I - Analysis. II - Solution and results

The effect of interleaves placed between the fiber-reinforced plies of multilayered composite laminates is investigated analytically, with the composite approximated by linearly elastic isotropic media consisting of half-planes separated from a finite-width layer by thin interleaves of low extensional and shear moduli. Integral equations are derived which describe the effect of interleaves on damage growth. Numerical solutions are then developed which are based on Gaussian integration techniques and Hadamard's concept of differentiation of Cauchy integrals. Critical stresses and displacements are calculated to investigate the influence of relative elastic properties and geometry of the interleafed composite on interface debonding, extent and suppression and delamination, and stress relieving in the undamaged plies.

Kaw, A. K.↗
NASA NTRS1991

An analysis of the Hubble Space Telescope fine guidance sensor fine lock mode

There are two guiding modes of the Hubble Space Telescope (HST) used for the acquisition of astronomical data by one of its six scientific instruments. The more precise one is called Fine Lock. Command and control problems in the onboard electronics has limited Fine Lock to brighter stars, V less than 13.0 mag, instead of fulfilling its goal of V = 14.5 mag. Consequently, the less precise guiding mode of Coarse Track (approximately 40 milli-arc seconds) has to be used fairly frequently. Indeed, almost half of the scientific observations to have been made with the HST will be compromised. The only realistic or extensive simulations of the Fine Lock guidance mode are reported. The theoretical analysis underlying the Monte Carlo experiments and the numerical computations clearly show both that the control electronics are severely under-engineered and how to adjust the various control parameters to successfully extend Fine Lock guiding performance back to V = 14.0 mag and sometimes beyond.

Taff, L. G.↗
NASA NTRS1993

DSMC and continuum analyses of low-density nozzle flow

Two different approaches, the direct-simulation Monte Carlo (DSMC) method based on molecular gas dynamics and a finite-volume approximation of the Navier-Stokes equations, which are based on continuum gas dynamics, are employed in the analysis of a low-density gas flow in a small converging-diverging nozzle. The fluid experiences various kinds of flow regimes including continuum, slip, transition, and free-molecular. Results from the two numerical methods are compared with Rothe's experimental dam, in which density and rotational temperature variations along the centerline and at various locations inside a low density nozzle were measured by the electron-beam fluorescence technique. The continuum approach showed good agreement with the experimental data as far as density is concerned. The results from the DSMC method showed good agreement with the experimental data both in the density and the rotational temperature. It is also shown that the simulation parameters, such as the gas/surface interaction model, the energy exchange model between rotational and translational modes, and the viscosity temperature exponent, have substantial effects on the results of the DSMC method.

Chung, Chan-Hong↗
NASA NTRS2004

MIT Haystack Observatory Analysis Center Report

The data from twenty years of the NCEP numerical weather model have been used to calculate the IMF hydrostatic mapping function for several sites distributed in latitude from -66 degrees to +78 degrees. Comparison of heights estimated with the NMF hydrostatic mapping function demonstrates that using NMFh results in height errors at annual and semi-annual periods with amplitudes as large as approximately 8 mm and 4 mm, respectively, when data down to 5 degrees are included. The errors are smallest at the equator and increase towards the poles.

Niell, Arthur↗
DOE OSTI2024

Variational Quantum Algorithms for Semidefinite Programming

A semidefinite program (SDP) is a particular kind of convex optimization problem with applications in operations research, combinatorial optimization, quantum information science, and beyond. In this work, we propose variational quantum algorithms for approximately solving SDPs. For one class of SDPs, we provide a rigorous analysis of their convergence to approximate locally optimal solutions, under the assumption that they are weakly constrained (i.e., N$\gg$M, where N is the dimension of the input matrices and M is the number of constraints). We also provide algorithms for a more general class of SDPs that requires fewer assumptions. Finally, we numerically simulate our quantum algorithms for applications such as MaxCut, and the results of these simulations provide evidence that convergence still occurs in noisy settings.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2015

Capillary Channel Flow (CCF) EU2-02 on the International Space Station (ISS): An Experimental Investigation of Passive Bubble Separations in an Open Capillary Channel

It would be signicantly easier to design fluid systems for spacecraft if the fluid phases behaved similarly to those on earth. In this research an open 15:8 degree wedge-sectioned channel is employed to separate bubbles from a two-phase flow in a microgravity environment. The bubbles appear to rise in the channel and coalesce with the free surface in much the same way as would bubbles in a terrestrial environment, only the combined effects of surface tension, wetting, and conduit geometry replace the role of buoyancy. The host liquid is drawn along the channel by a pump and noncondensible gas bubbles are injected into it near the channel vertex at the channel inlet. Control parameters include bubble volume, bubble frequency, liquid volumetric flow rate, and channel length. The asymmetrically confined bubbles are driven in the cross-flow direction by capillary forces until they at least become inscribed within the section or until they come in contact with the free surface, whereupon they usually coalesce and leave the flow. The merging of bubbles enhances, but does not guarantee, the latter. The experiments are performed aboard the International Space Station as a subset of the Capillary Channel Flow experiments. The flight hardware is commanded remotely and continuously from ground stations during the tests and an extensive array of experiments is conducted identifying numerous bubble flow regimes and regime transitions depending on the ratio and magnitude of the gas and liquid volumetric flow rates. The breadth of the publicly available experiments is conveyed herein primarily by narrative and by regime maps, where transitions are approximated by simple expressions immediately useful for the purposes of design and deeper analysis.

fluid dynamics↗
NASA NTRS1975

Numerical nonlinear inelastic analysis of stiffened shells of revolution. Volume 3: Engineer's program manual for STARS-2P digital computer program

Engineering programming information is presented for the STARS-2P (shell theory automated for rotational structures-2P (plasticity)) digital computer program, and FORTRAN 4 was used in writing the various subroutines. The execution of this program requires the use of thirteen temporary storage units. The program was initially written and debugged on the IBM 370-165 computer and converted to the UNIVAC 1108 computer, where it utilizes approximately 60,000 words of core. Only basic FORTRAN library routines are required by the program: sine, cosine, absolute value, and square root.

Svalbonas, V.↗
NASA NTRS1984

Numerical studies of unsteady transonic flow over an oscillating airfoil

A finite-difference solution to the Navier-Stokes equations combined with a time-varying grid-generation technique was used to compute unsteady transonic flow over an oscillating airfoil. These computations were compared with experimental data (obtained at Ames Research Center) which form part of the AGARD standard configuration for aeroelastic analysis. A variety of approximations to the full Navier-Stokes equations was used to determine the effect of frequency, shock-wave motion, flow separation, and airfoil geometry on unsteady pressures and overall air loads. Good agreement is shown between experiment and theory with the limiting factor being the lack of a reliable turbulence model for high-Reynolds-number, unsteady transonic flows.

Chyu, W. J.↗
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.↗