Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical continuation”

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 199 records · Page 11

Reliability of CHAMP Anomaly Continuations

CHAMP is recording state-of-the-art magnetic and gravity field observations at altitudes ranging over roughly 300 - 550 km. However, anomaly continuation is severely limited by the non-uniqueness of the process and satellite anomaly errors. Indeed, our numerical anomaly simulations from satellite to airborne altitudes show that effective downward continuations of the CHAMP data are restricted to within approximately 50 km of the observation altitudes while upward continuations can be effective over a somewhat larger altitude range. The great unreliability of downward continuation requires that the satellite geopotential observations must be analyzed at satellite altitudes if the anomaly details are to be exploited most fully. Given current anomaly error levels, joint inversion of satellite and near- surface anomalies is the best approach for implementing satellite geopotential observations for subsurface studies. We demonstrate the power of this approach using a crustal model constrained by joint inversions of near-surface and satellite magnetic and gravity observations for Maude Rise, Antarctica, in the southwestern Indian Ocean. Our modeling suggests that the dominant satellite altitude magnetic anomalies are produced by crustal thickness variations and remanent magnetization of the normal polarity Cretaceous Quiet Zone.

vonFrese, Ralph R. B.↗

A singularity free approach to post glacial rebound calculations

Calculating the post glacial response of a viscoelastic Earth model using the exponential decay normal mode technique leads to intrinsic singularities if viscosity varies continuously as a function of radius. We develop a numerical implementation of the Complex Real Fourier transform (CRFT) method as an accurate and stable procedure to avoid these singularities. Using CRFT, we investigate the response of a set of Maxwell Earth models to surface loading. We find that the effect of expanding a layered viscosity structure into a continuously varying structure is to destroy the modes associated with the boundary between layers. Horizontal motion is more sensitive than vertical motion to the viscosity structure just below the lithosphere. Horizontal motion is less sensitive to the viscosity of the lower mantle than the vertical motion is. When the viscosity increases at 670 km depth by a factor of about 60, the response of the lower mantle is close to its elastic limit. Any further increase of the viscosity contrast at 670 km depth or further increase of viscosity as a continuous function of depth starting from 670 km depth is unlikely to be resolved.

Fang, Ming↗

Continuation Methods for Qualitative Analysis of Aircraft Dynamics

A class of numerical methods for constructing bifurcation curves for systems of coupled, non-linear ordinary differential equations is presented. Foundations are discussed, and several variations are outlined along with their respective capabilities. Appropriate background material from dynamical systems theory is presented.

Cummings, Peter A.↗

Analysis of two-dimensional incompressible flows by a subsurface panel method

A new approach to panel methods is explored for two-dimensional steady incompressible flows. The method uses linear distributions of sources and vortices on straight-line panels, but satisfies boundary conditions on the actual body surface, at nodes that are also end points of the panels. The result is continuity in body-surface velocity distribution, without recourse to numerical quadrature for the velocity influence coefficients. The method is unusually sensitive to the distribution of the nodes. For example, it almost always fails to give acceptable results when the nodes are distributed randomly. However, the continuity of the velocity distribution makes possible a unique node redistribution scheme, which may be iterated to give accurate results reliably.

Moran, J.↗

The Prediction of Ducted Fan Engine Noise Via a Boundary Integral Equation Method

A computationally efficient Boundary Integral Equation Method (BIEM) for the prediction of ducted fan engine noise is presented. The key features of the BIEM are its versatility and the ability to compute rapidly any portion of the sound field without the need to compute the entire field. Governing equations for the BIEM are based on the assumptions that all acoustic processes are linear, generate spinning modes, and occur in a uniform flow field. An exterior boundary value problem (BVP) is defined that describes the scattering of incident sound by an engine duct with arbitrary profile. Boundary conditions on the duct walls are derived that allow for passive noise control treatment. The BVP is recast as a system of hypersingular boundary integral equations for the unknown duct surface quantities. BIEM solution methodology is demonstrated for the scattering of incident sound by a thin cylindrical duct with hard walls. Numerical studies are conducted for various engine parameters and continuous portions of the total pressure field are computed. Radiation and duct propagation results obtained are in agreement with the classical results of spinning mode theory for infinite ducts.

Boundary integral equation↗

Numerical simulation of the reactive flow in advanced (HSR) combustors using KIVA-2

Recent work has been done with the goal of establishing ultralow emission aircraft gas turbine combustors. A significant portion of the effort is the development of three dimensional computational combustor models. The KIVA-II computer code which is based on the Implicit Continuous Eulerian Difference mesh Arbitrary Lagrangian Eulerian (ICED-ALE) numerical scheme is one of the codes selected by NASA to achieve these goals. This report involves a simulation of jet injection through slanted slots within the Rich burn/Quick quench/Lean burn (RQL) baseline experimental rig. The RQL combustor distinguishes three regions of combustion. This work specifically focuses on modeling the quick quench mixer region in which secondary injection air is introduced radially through 12 equally spaced slots around the mixer circumference. Steady state solutions are achieved with modifications to the KIVA-II program. Work currently underway will evaluate thermal mixing as a function of injection air velocity and angle of inclination of the slots.

Winowich, Nicholas S.↗

Perturbation theory for charged-particle transport in one dimension

Perturbation theory, when applied to charged-particle transport, generates a series solution that requires a double quadrature per term. The continuity of higher-order terms leads to numerical evaluation of the series. The high rate of convergence of the series makes the method a practical tool for charged-particle transport problems. The coupling of the neutron component in the case of proton transport in tissue does not greatly alter the rate of convergence. The method holds promise for a practical high-energy proton transport theory.

Wilson, J. W.↗

Composition of the nighttime ionospheric F 1 region near the magnetic equator

The effects of vertical E x B transport on NO(+), O2(+) and O(+) densities in the nighttime equatorial ionospheric F 1 region are investigated. Ion densities are calculated as functions of altitude, latitude and local time by the numerical solution of coupled, time-dependent ion continuity equations, taking into account production, loss by charge exchange and dissociative recombination and transport by diffusion and E x B drift. The results of the calculations are compared with measurements of NO(+), O2(+) and O(+) ion densities obtained at low altitudes by a mass spectrometer on board the Atmospheric Explorer C satellite, and are found to be consistent with the observations, suggesting that in the equatorial region, vertical transport by E x B drift is primarily responsible for producing the observed NO(+), O2(+) and O(+) density profiles. In addition, the reaction of O2(+) with N(4S) is found to be an important sink for O2(+) and a source of NO(+) ions. Implications of the observed and calculated near constancy of electron and ion densities with altitude when NO(+) is the dominant ion on the growth of large-scale irregularities are also considered

Anderson, D. N.↗

A refined shear deformation theory for the analysis of laminated plates

A refined, third-order plate theory that accounts for the transverse shear strains is presented, the Navier solutions are derived for certain simply supported cross-ply and antisymmetric angle-ply laminates, and finite-element models are developed for general laminates. The new theory does not require the shear correction factors of the first-order theory (i.e., the Reissner-Mindlin plate theory) because the transverse shear stresses are represented parabolically in the present theory. A mixed finite-element model that uses independent approximations of the generalized displacements and generalized moments, and a displacement model that uses only the generalized displacements as degrees of freedom are developed. The displacement model requires C sup 1-continuity of the transverse deflection across the inter-element boundaries, whereas the mixed model requires a C sup 0-element. Also, the mixed model does not require continuous approximations (between elements) of the bending moments. Numerical results are presented to show the accuracy of the present theory in predicting the transverse stresses. Numerical results are also presented for the nonlinear bending of plates, and the results compare well with the experimental results available in the literature.

Reddy, J. N.↗

Solar transition region and coronal response to heating rate perturbations

Observations of Doppler shifts in UV emission lines formed in the solar transition region show continual plasma downflows and impulsive plasma upflows. Using numerical simulations, the authors examine the conjecture that areas of downflowing plasma are the base regions of coronal loops in which the heating is gradually decreasing and that areas of upflowing plasma are the base regions of coronal loops in which the heating rate is gradually increasing. Beginning with a coronal loop in equilibrium, the heating rate is reduced on time scales of 100, 1000, and 2000 s to 10 percent and 1 percent of the initial value, and the loop is allowed to evolve to a new equilibrium. The heating rate for the cooled models is then increased back to the initial value on the same time scales. While significant mass motions do develop in the simulations, both the emission measure and the velocity at 100,000 K do not show the characteristics present in UV observations.

Mariska, John T.↗

The application of TSIM software to act design and analysis on flexible aircraft

The TSIM software is described. This is a package which uses an interactive FORTRAN-like simulation language for the simulation on nonlinear dynamic systems and offers facilities which include: mixed continuous and discrete time systems, time response calculations, numerical optimization, automatic trimming of nonlinear aircraft systems, and linearization of nonlinear equations for eigenvalues, frequency responses and power spectral response evaluation. Details are given of the application of TSIM to the analysis of aeroelastic systems under the RAE Farborough extension FLEX-SIM. The aerodynamic and structural data for the equations of motion of a flexible aircraft are prepared by a preprocessor program for incorporation in TSIM simulations. Within the simulation, the flexible aircraft model may then be selected interactively for different flight conditions and modal reduction techniques applied. The use of FLEX-SIM is demonstrated by an example of the flutter prediction for a simple aeroelastic model. By utilizing the numerical optimization facility of TSIM, it is possible to undertake identification of required parameters in the TSIM model within the simulation.

Kaynes, Ian W.↗

Time-dependent shock acceleration of energetic electrons including synchrotron losses

The present investigation of the time-dependent particle acceleration problem in strong shocks, including synchrotron radiation losses, solves the transport equation analytically by means of Laplace transforms. The particle distribution thus obtained is then transformed numerically into real space for the cases of continuous and impulsive injections of particles at the shock. While in the continuous case the steady-state spectrum undergoes evolution, impulsive injection is noted to yield such unpredicted features as a pile-up of high-energy particles or a steep power-law with time-dependent spectral index. The time-dependent calculations reveal varying spectral shapes and more complex features for the higher energies which may be useful in the interpretation of outburst spectra.

Fritz, Klaus-Dieter↗

Plate and butt-weld stresses beyond elastic limit, material and structural modeling

Ultimate safety factors of high performance structures depend on stress behavior beyond the elastic limit, a region not too well understood. An analytical modeling approach was developed to gain fundamental insights into inelastic responses of simple structural elements. Nonlinear material properties were expressed in engineering stresses and strains variables and combined with strength of material stress and strain equations similar to numerical piece-wise linear method. Integrations are continuous which allows for more detailed solutions. Included with interesting results are the classical combined axial tension and bending load model and the strain gauge conversion to stress beyond the elastic limit. Material discontinuity stress factors in butt-welds were derived. This is a working-type document with analytical methods and results applicable to all industries of high reliability structures.

Verderaime, V.↗

Physics of heavy ions (1989-1990)

The results from studies on polar wind ion heating due to kinetic ion beam instabilities and the effects of such ion heating on the outflow of O(+) in the polar wind are presented and discussed. First, the linear instabilities associated with an O(+) and H(+) polar wind plasma in the presence of O(+) and H(+) beams for a range of O(+)/H(+) beam densities, T(sub e)/T(sub i), and ion beam speeds were examined. Then, nonlinear heating of the polar wind ions was studied, using numerical simulations. The O(+) and H(+) polar wind ions were modeled by isotropic Maxwellian distributions, and the electrons, O(+) beams, and H(+) beams were modeled by drifting Maxwellian distributions. The effects of the kinetic ion heating on the outflow of the polar wind ions were examined from the ionosphere, using a time-dependent hydrodynamic model. A numerical code to solve the O(+) and H(+) continuity and momentum equations in a flux tube from ionospheric to magnetospheric altitudes were developed. The effects of ion heating were included by allowing for the altitudinal variation of the ion temperatures in the momentum equation. The ion temperature profiles were specified based on the ion heating characteristics found from previous kinetic simulations. It was assumed that heating occurred above 1500 km and increased to a saturated value of temperature that was obtained directly from the kinetic simulation study. The characteristics of the dynamical polar wind without ion heating were studied, and a flux tube on closed field lines that suddenly became open at t = 0 was simulated. Then, the effects of ion heating were included. To gain some physical insight, two limiting cases were considered: preferential H(+) heating and preferential O(+) heating. How O(+) heating can lead to enhanced polar wind O(+) fluxes in the polar magnetosphere is shown.

Ashour-Abdalla, M.↗

Improved finite-element methods for rotorcraft structures

An overview of the research directed at improving finite-element methods for rotorcraft airframes is presented. The development of a modification to the finite element method which eliminates interelement discontinuities is covered. The following subject areas are discussed: geometric entities, interelement continuity, dependent rotational degrees of freedom, and adaptive numerical integration. This new methodology is being implemented as an anisotropic, curvilinear, p-version, beam, shell, and brick finite element program.

Hinnant, Howard E.↗

A Computational and Experimental Investigation of Shear Coaxial Jet Atomization

The instability and subsequent atomization of a viscous liquid jet emanated into a high-pressure gaseous surrounding is studied both computationally and experimentally. Liquid water issued into nitrogen gas at elevated pressures is used to simulate the flow conditions in a coaxial shear injector element relevant to liquid propellant rocket engines. The theoretical analysis is based on a simplified mathematical formulation of the continuity and momentum equations in their conservative form. Numerical solutions of the governing equations subject to appropriate initial and boundary conditions are obtained via a robust finite difference scheme. The computations yield real-time evolution and subsequent breakup characteristics of the liquid jet. The experimental investigation utilizes a digital imaging technique to measure resultant drop sizes. Data were collected for liquid Reynolds number between 2,500 and 25,000, aerodynamic Weber number range of 50-500 and ambient gas pressures from 150 to 1200 psia. Comparison of the model predictions and experimental data for drop sizes at gas pressures of 150 and 300 psia reveal satisfactory agreement particularly for lower values of investigated Weber number. The present model is intended as a component of a practical tool to facilitate design and optimization of coaxial shear atomizers.

Ibrahim, Essam A.↗

"Fly-by-Wireless": A Revolution in Aerospace Vehicle Architecture for Instrumentation and Control

Aerospace vehicle programs have always counted on the cables and connectors to provide power, grounding, data and time synchronization throughout a vehicle's life-cycle. Even with numerous improvements, wiring and connector problems and sensors continue to be key failure points, causing many hours of troubleshooting and replacement. Costly flight delays have been precipitated by the need to troubleshoot cables/connections, and/or repair a sensor. Wiring continues to be too expensive to remove once it is installed, even with the weight penalties. Miles of test instrumentation and low flight sensor wires still plague the aerospace industry. New technology options for data connectivity, processing and micro/nano manufacturing are making it possible to retrofit existing vehicles, like the Space Shuttle. New vehicles can now develop architectures that provide for and take advantage of alternatives to wired connectivity. This project motivates the aerospace industry and technology providers to establish: (1) A new emphasis for system engineering approaches to reduce cables and connectors. (2) Provisions for modularity and accessibility in the vehicle architecture. (3) A set of technologies that support alternatives to wired connectivity.

Studor, George↗