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 379 records · Page 21

Numerical Characterization of Astronaut CaOx Renal Stone Incidence Rates to Quantify In-flight and Post-flight Relative Risk

Changes in urine chemistry potentially alter the risk of renal stone formation in astronauts. Quantifying spaceflight renal stone incidence risk compared to pre-flight levels remains a significant challenge for assessing the appropriate vehicle, mission, and countermeasure design. A computational biochemistry model representing CaOx crystal precipitation, growth, and agglomeration is combined with a probabilistic analysis to predict the in- and post-flight CaOx renal stone incidence risk ratio (IRR) relative to pre-flight values using 1517 astronaut 24-hour urine chemistries. Our simulations predict that in-flight fluid intake alone would need to increase from current prescriptions of 2.0 - 2.5 L/day to approximately 3.2 L/day to approach CaOx IRR of the pre-flight population. Bone protective interventions would reduce CaOx risk to pre-flight levels if Ca excretion alone is reduced to < 150 mg/day or if current levels are diminished to 190mg/day in combination with increasing fluid intake to 2.5 - 2.7 L/day. This analysis provides a quantitative risk assessment that can influence the critical balance between engineering and astronaut health requirements.

renal stones↗

The Surface Chemical Composition of Lunar Samples and Its Significance for Optical Properties

The surface iron, titanium, calcium, and silicon concentration in numerous lunar soil and rock samples was determined by Auger electron spectroscopy. All soil samples show a large increase in the iron to oxygen ratio compared with samples of pulverized rock or with results of the bulk chemical analysis. A solar wind simulation experiment using 2 keV energy alpha -particles showed that an ion dose corresponding to approximately 30,000 years of solar wind increased the iron concentration on the surface of the pulverized Apollo 14 rock sample 14310 to the concentration measured in the Apollo 14 soil sample 14163, and the albedo of the pulverized rock decreased from 0.36 to 0.07. The low albedo of the lunar soil is related to the iron + titanium concentration on its surface. A solar wind sputter reduction mechanism is discussed as a possible cause for both the surface chemical and optical properties of the soil.

Gold, T.↗

The evolution and discharge of electric fields within a thunderstorm

An analysis of the present three-dimensional thunderstorm electrical model and its finite-difference approximations indicates unconditional stability for the discretization that results from the approximation of the spatial derivatives by a box-schemelike method and of the temporal derivative by either a backward-difference or Crank-Nicholson scheme. Lightning propagation is treated through numerical techniques based on the inverse-matrix modification formula and Cholesky updates. The model is applied to a storm observed at the Kennedy Space Center in 1978, and numerical comparisons are conducted between the model and the theoretical results obtained by Wilson (1920) and Holzer and Saxon (1952).

Hager, William W.↗

An analysis of carrier phase jitter in an MPSK receiver utilizing map estimation

The use of 8 and 16 PSK TCM to support satellite communications in an effort to achieve more bandwidth efficiency in a power-limited channel has been proposed. This project addresses the problem of carrier phase jitter in an M-PSK receiver utilizing the high SNR approximation to the maximum aposteriori estimation of carrier phase. In particular, numerical solutions to the 8 and 16 PSK self-noise and phase detector gain in the carrier tracking loop are presented. The effect of changing SNR on the loop noise bandwidth is also discussed. These data are then used to compute variance of phase error as a function of SNR. Simulation and hardware data are used to verify these calculations. The results show that there is a threshold in the variance of phase error versus SNR curves that is a strong function of SNR and a weak function of loop bandwidth. The M-PSK variance thresholds occur at SNR's in the range of practical interest for the use of 8 and 16-PSK TCM. This suggests that phase error variance is an important consideration in the design of these systems.

Osborne, William P.↗

Modeling and numerical simulation of microwave pulse propagation in air breakdown environment

It is shown that ionization occurs wherever the field intensity of the pulse exceeds the local breakdown threshold field of the background air. The produced plasma then attenuates the pulse and gives rise to a tail erosion phenomenon that plays the primary role in limiting the energy transfer of the pulse from source to destination. A theoretical model describing the propagation of an intense microwave pulse in an air breakdown environment is developed that includes the possible focusing effect introduced by either using phase array antennas or the other arrangements. The self-consistent description of the propagation process is provided by a set of two modal equations. These include a continuity equation (Poynting's equation) for the energy density of the pulse and a rate equation of the electron density. A forward wave approximation is used to simplify Poynting's equation, and a semiempirical formula is used for the ionization frequency, nu(sub i). This frequency provides the coupling between the two modal equations, and is used to express the electron rate equation explicitly. In terms of the relevant parameters of the atmosphere, these two equations are normalized for numerical analysis of pulse propagation in the atmosphere. The dependencies of the propagation characteristics of the pulse on intensity, frequency, width, and shape of the pulse are determined. The numerical simulations lead to a useful empirical relation p(exp 3)w = alpha = constant, where p and w are the incident power and width of the pulse and alpha depends on the percentage of pulse energy transferred from the source point to a destined position. The density distribution of the pulse's self-generated plasma is also evaluated. The results also show that for ionization caused by a single unfocused microwave pulse transmitted upwards from the ground, the maximum electron density produced at, for example, 50 km altitude is limited by the tail erosion effect to below 10(exp 6) cm(exp -3). Repetitive pulse and focused beam approaches are also examined. Both approaches can increase the maximum electron density by no more than an order of magnitude. A scheme using two obliquely propagating pulses intersecting at the destined height, e.g. 50 km, is considered. It is shown that the electron density generated at the lowest intersecting position can easily reach a value of 6.6 x 10(exp 8) cm(exp -3), which is considered to be high enough for artificial ionospheric mirror (AIM) application.

Kim, J.↗

Iterative Discrete Ordinates Solution of the Equation for the Surface-Reflected Radiance

This paper presents a new method of numerical solution of the integral equation for the radiance reflected from an anisotropic surface. The equation relates the radiance at the surface level with BRDF and solutions of the standard radiative transfer problems for a slab with no reflection on its surfaces. It is also shown that the kernel of the equation satisfies the condition of the existence of a unique solution and the convergence of the successive approximations to that solution. The developed method features two basic steps: discretization on a 2D quadrature, and solving the resulting system of algebraic equations with successive over-relaxation method based on the Gauss-Seidel iterative process. Presented numerical examples show good coincidence between the surface-reflected radiance obtained with DISORT and the proposed method. Analysis of contributions of the direct and diffuse (but not yet reflected) parts of the downward radiance to the total solution is performed. Together, they represent a very good initial guess for the iterative process. This fact ensures fast convergence. The numerical evidence is given that the fastest convergence occurs with the relaxation parameter of 1 (no relaxation). An integral equation for BRDF is derived as inversion of the original equation. The potential of this new equation for BRDF retrievals is analyzed. The approach is found not viable as the BRDF equation appears to be an ill-posed problem, and it requires knowledge the surface-reflected radiance on the entire domain of both Sun and viewing zenith angles.

Alexander Radkevich↗

An analysis of auroral E region neutral winds based on incoherent scatter radar observations at Chatanika

The physical origin of auroral E-region neutral winds during geomagnetic disturbances is investigated by comparing incoherent-backscatter radar observations made at Chatanika, Alaska, during such disturbances with detailed theoretical calculations of neutral velocities for these conditions. In the analysis, an ion-drag collisional-coupling model is adopted as the primary driving force for the neutral winds. Electric fields and electron number densities observed by the radar facility are used in the calculations, and the coupled ion and neutral momentum equations are solved numerically as functions of altitude and time in the ion-drag approximation. The altitude structure of the velocity field is examined, and the detailed response of the neutral gas to the ion-drag mechanism is delineated. Weighted altitude averages of the calculated velocities are obtained, the accuracy of the procedure for deriving neutral winds from incoherent-scatter radar observations is evaluated, and the comparison is carried out. It is concluded that both ion drag and pressure gradients appear to be required to explain the neutral winds observed in the auroral E-region during geomagnetic disturbances.

Comfort, R. H.↗

Three-dimensional turbulent boundary layer development and separation in V/STOL engine inlets at incidence with small-cross flow and curvature influences

The study concerns the influence of the three-dimensional cross flows on the compressible turbulent boundary layer development and flow separation prediction around V/STOL engine inlets at high incidence. The governing equations for the three-dimensional boundary layer flow with small-cross approximation are solved numerically on an intrinsic streamline coordinate system. Results are presented to illustrate the effects of small cross-flow, compressibility and streamline curvatures on the flow. Comparisons of the results with the wind tunnel data for scaled model and with data obtained from another existing compressible axisymmetric turbulent boundary layer scheme are included in the analysis.

Chou, D. C.↗

Design and calibration of a rocket-borne electron spectrometer for investigation of particle ionization in the nighttime midlatitude E region

An explanation was developed for the formation, near midnight at midlatitudes, of a broad electron density layer extending approximately from 120 to 180 km and usually referred to as the intermediate E layer. The responsible mechanism is believed to be the converging vertical ion drifts resulting from winds of the solar semidiurnal tide. Numerical solutions of the continuity equation appropriate to the intermediate layer is described for particular models of ion drift, diffusion coefficents, and ionization production. Analysis of rocket observations of the layer show that the ionization rate is highly correlated with the planetary geomagnetic index, K sub p. Particle flux measurements support the idea that energetic electrons are the principal source of this ionization. A semiconductor spectrometer experiment for investigation of the particle flux, spectrum, and angular properties was designed and successfully flown on a Nike Apache rocket. A detailed description of the theory, design, and calibration of the experiment and some preliminary results presented.

Voss, H. D.↗

Finite element representations for thin shell instability analysis

The development of finite element calculational procedures for thin-shell instability analysis has involved the definition of appropriate element representations (i.e. the geometric form of the element and the approximation of displacement and/or stress), constitutive expressions, and computational algorithms. Among these, the status of thin-shell finite element representations remains unsettled. This paper, therefore, emphasizes recent developments in the basic aspects of thin-shell finite element analysis and discusses one simplified approach in more detail. Formulative and computational procedures for elastic instability analysis are then summarized, and numerical results are shown for the simplified shell element formulation.

Gallagher, R. H.↗

On damping mechanisms in beams

A partial differential equation model of a cantilevered beam with a tip mass at its free end is used to study damping in a composite. Four separate damping mechanisms consisting of air damping, strain rate damping, spatial hysteresis and time hysteresis are considered experimentally. Dynamic tests were performed to produce time histories. The time history data is then used along with an approximate model to form a sequence of least squares problems. The solution of the least squares problem yields the estimated damping coefficients. The resulting experimentally determined analytical model is compared with the time histories via numerical simulation of the dynamic response. The procedure suggested here is compared with a standard modal damping ratio model commonly used in experimental modal analysis.

Banks, H. T.↗

Numerical simulation of 3-D Benard convection with gravitational modulation

In this numerical study, randomly and sinusoidally modulated gravitational fields imposed on three-dimensional Rayleigh-Benard convection are investigated in an effort to understand the effects of vibration (G-Jitter) on fluid systems. The time-dependent, Navier-Stokes equations and the energy equation with Boussinesq approximations are solved by a semi-implicit, pseudospectral procedure. An analysis of energy balances indicates that with increasing modulation amplitude, transition from synchronous to relaxation oscillation goes through the subharmonic response. Random modulations are found to be less stabilizing than sinusoidal and are shown to impose three-dimensionality on the flow for some parameter ranges both at terrestrial and zero base gravity conditions.

Biringen, S.↗

On damping mechanisms in beams

A partial differential equation model of a cantilevered beam with a tip mass at its free end is used to study damping in a composite. Four separate damping mechanisms consisting of air damping, strain rate damping, spatial hysteresis and time hysteresis are considered experimentally. Dynamic tests were performed to produce time histories. The time history data is then used along with an approximate model to form a sequence of least squares problems. The solution of the least squares problem yields the estimated damping coefficients. The resulting experimentally determined analytical model is compared with the time histories via numerical simulation of the dynamic response. The procedure suggested here is compared with a standard modal damping ratio model commonly used in experimental modal analysis.

Banks, H. T.↗

Small-Body Extensions for the Satellite Orbit Analysis Program (SOAP)

An extension to the SOAP software allows users to work with tri-axial ellipsoid-based representations of planetary bodies, primarily for working with small, natural satellites, asteroids, and comets. SOAP is a widely used tool for the visualization and analysis of space missions. The small body extension provides the same visualization and analysis constructs for use with small bodies. These constructs allow the user to characterize satellite path and instrument cover information for small bodies in both 3D display and numerical output formats. Tri-axial ellipsoids are geometric shapes the diameters of which are different in each of three principal x, y, and z dimensions. This construct provides a better approximation than using spheres or oblate spheroids (ellipsoids comprising two common equatorial diameters as a distinct polar diameter). However, the tri-axial ellipsoid is considerably more difficult to work with from a modeling perspective. In addition, the SOAP small-body extensions allow the user to actually employ a plate model for highly irregular surfaces. Both tri-axial ellipsoids and plate models can be assigned to coordinate frames, thus allowing for the modeling of arbitrary changes to body orientation. A variety of features have been extended to support tri-axial ellipsoids, including the computation and display of the spacecraft sub-orbital point, ground trace, instrument footprints, and swathes. Displays of 3D instrument volumes can be shown interacting with the ellipsoids. Longitude/latitude grids, contour plots, and texture maps can be displayed on the ellipsoids using a variety of projections. The distance along an arbitrary line of sight can be computed between the spacecraft and the ellipsoid, and the coordinates of that intersection can be plotted as a function of time. The small-body extension supports the same visual and analytical constructs that are supported for spheres and oblate spheroids in SOAP making the implementation of the more complex algorithms largely transparent to the user.

Carnright, Robert↗

Theoretical analysis of rough-surface shadowing from point-source radiation.

The shadowing function, defined as the probability that a point on a rough surface will be illuminated by radiation from a source at fixed height above the surface, is considered. An exact expression for this probability is derived and evaluated approximately when the surface elevation is a normal random process. The effect of surface correlation is considered and shown to be appreciable. Several numerical computations of the shadowing function are presented.

Hardin, J. C.↗

Incremental analysis of large elastic deformation of a rotating cylinder

The effect of finite deformation upon a rotating, orthotropic cylinder was investigated using a general incremental theory. The incremental equations of motion are developed using the variational principle. The governing equations are derived using the principle of virtual work for a body with initial stress. The governing equations are reduced to those for the title problem and a numerical solution is obtained using finite difference approximations. Since the problem is defined in terms of one independent space coordinate, the finite difference grid can be modified as the incremental deformation occurs without serious numerical difficulties. The nonlinear problem is solved incrementally by totaling a series of linear solutions.

Buchanan, G. R.↗

Approximation methods for the solution of inverse problems in lake and sea sediment core analysis

A theoretical model employing one-dimensional (depth) transport equations to describe vertical redistribution of ocean-floor and lake-floor sediment (particulates, volcanic ash, microtektites, or radioactive tracers) by episodic and nonepisodic events including bioturbation is developed analytically and demonstrated. The principles underlying the model are explained; the model equations are derived; the inverse problem of identifying the depth-dependent bioturbation coefficient is addressed; two approximation theorems are presented; and numerical results for two sample problems are presented graphically. It is suggested that compatification, porosity effects, and depth-dependent sedimentation be taken into account when formulating future models.

Banks, H. T.↗

An investigation of the transient thermal analysis of spur gears

A finite element computer program is developed for evaluating the transient behavior of surface temperature in high performance spur gears. The time dimension is implemented using two and three point finite difference schemes. The different schemes are provided for the purpose of numerical stability and convergence studies. A detailed explanation of the gear cooling process leading to the establishment of a modified Blok model is also included. Other conventional models for approximating the heat transfer coefficients are available for comparison. Preliminary results are given showing snap shots of gear temperature contours at the initial stages of tooth engagement.

El-Bayoumy, L. E.↗