SEARCH · Engineering Papers
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.
A quantitative model is presented, which describes the ocean surface as an ensemble of flat facets with a normal distribution of slopes. The variance of the slope distribution is linearly related to frequency up to 35 GHz and constant at higher frequencies. These facets are partially covered with an absorbing nonpolarized foam layer. Experimental evidence is presented for this model.
Finite element analysis is applied to study the large deflection and stress of a thin-walled pressurized torus loaded by normal contact with a plane. Experimentally determined stresses are compared to those resulting from finite element analysis at various meridional and circumferential coordinates of the torus, including the footprint area. Circumferential strains compare favorably while meridional strains are higher in the finite element analysis, probably due to slippage of the boundary at the rim. The STAGS finite element computer program was used. The purpose of the study was to evaluate the various program options for structural loading, for material modeling, for stress analysis, and for grid refinement. The experimental model, a thin-walled rubber tube mounted on a steel cylindrical rim, provided measured results to compare with various analytical trade-offs. It was found that there was almost no difference in predicted deflections or stress distributions between linear and nonlinear material description. However, the difference between linear analysis and that of geometric nonlinearities utilizing incremental loading was marked.
Numerical solutions of the instability wave pattern behind a harmonic point source in a Blasius boundary layer are used to form line sources by superposition. For infinite-length spanwise line sources of constant amplitude and phase, the result is just the two-dimensional normal mode of the same frequency; for a sinusoidal amplitude or linear phase distribution, the result is an oblique normal mode of the same spanwise wavenumber. A finite-length spanwise source simulates a vibrating ribbon. In a study of the influence of the source length on the downstream amplitude, the effective field of view is found to have a half-angle of about 16 deg. If the source tips are within this field, the amplitude may be either greater than or less than the comparable normal-mode amplitude, depending on the distance from the source and the spanwise location. For an oblique line source, the downstream wave development at each spanwise location is found to be close to, but not identical with, that of an oblique normal mode which originates at the source with the initial wave angle of the source and satisfies the irrotationality condition on the wavenumber vector.
The fully differential cross section for photon-electron pair production is integrated numerically over phase space. Results are obtained for the astrophysically interesting case in which the interaction between an ultrarelativistic electron and a soft photon results in electron-positron pair production. The positron spectrum is a function of the energies of both the photon and the electron, as well as the angle of interaction. It is found that the energy at which the positron distribution peaks is inversely proportional to the photon energy and independent of the electron energy. The positron spectrum is integrated once more over initial electron energies for a power-law energy distribution of primary electrons. The same procedure is repeated for the recoil particle; it is shown that the peak of the recoil energy distribution depends linearly on the energy of the primary electron. Finally, semianalytical expressions are obtained for the energy losses of the primary electrons.
The thermocapillary motion of spherical bubbles present in an unbounded liquid with a linear temperature distribution, when the Reynolds number and the Marangoni number are large is analyzed. Previous calculations of the terminal velocity performed for this parametric range did not take into complete consideration the thermal boundary layer present near the surface of the bubble. A scaling analysis is presented for this problem. The thermal boundary layer is analyzed by an integral method. The resulting terminal velocity is lower than the one previously calculated, though it is of the same order of magnitude.
For the recapture debugging design introduced by Nayak (1988) the problem of estimating the hitting rates of the faults remaining in the system is considered. In the context of a conditional likelihood, moment estimators are derived and are shown to be asymptotically normal and fully efficient. Fixed sample properties of the moment estimators are compared, through simulation, with those of the conditional maximum likelihood estimators. Properties of the conditional model are investigated such as the asymptotic distribution of linear functions of the fault hitting frequencies and a representation of the full data vector in terms of a sequence of independent random vectors. It is assumed that the residual hitting rates follow a log linear rate model and that the testing process is truncated when the gaps between the detection of new errors exceed a fixed amount of time.
The classical linear theory of electromagnetic wave growth in a warm plasma is considered for waves propagating parallel to a uniform ambient magnetic field. Wave-growth rates are calculated for ion-driven right-hand mode waves for Kappa and Maxwellian particle distribution functions and for various values of the spectral index, the temperature anisotropy, and the ratio of plasma pressure to magnetic pressure appropriate to the solar wind. When the anisotropy is low the wave growth is limited to frequencies below the proton gyrofrequency and the growth rate increases dramatically as the spectral index is reduced. The growth rate for any Kappa distribution greatly exceeds that for a Maxwellian with the same bulk properties. For large thermal anisotropy the growth rate from either distribution is greatly enhanced. The growth rates from a Kappa distribution are generally larger than for a Maxwellian distribution, and significant wave growth occurs over a broader range of frequencies.
The flyby missions of Voyagers 1 and 2 at Jupiter, Saturn and Uranus revealed intense waves above the electron gyrofrequency. Observation of waves at the upper hybrid frequency is often accompanied by power at adjacent electron Bernstein harmonics, and the relative power in these modes depends both on the density and temperature ratios of the cold background electron population and the hot magnetospheric electrons which drive the instability. A model of electron distributions which is consistent with observations is used for analysis of the excited waves, their dependence upon plasma parameters, and the time scales of the saturation processes. It is shown that in the presence of two-temperature electron distributions the linear excitation is due to a fluidlike coupling of two eigenmodes for perpendicular propagation and to kinetic destalization of oblique modes. The dependence of linear growth rates on propagation angle is presented, along with results from particle simulations. A quasi-linear diffusion time for relaxing the hot electron loss cone is calculated and compared with simulation results. This time scale is faster than for local saturation by heating the cold population, and also the convective amplification time scale, suggesting that the waves saturate at quasi-linear levels, while being convectively localized to the equatorial regions of the outer planetary magnetospheres.
The linear size distribution of 218 galaxies in the IRAS infrared bright galaxy sample is studied by examining new measurements of the 1.49 Hz radio continuum emitting regions. We find that the radio surface brightness varies over a wide range at a given luminosity, but its median scales approximately as luminosity to the power of 6/5. Consequently, the median effective radio size of galaxies in a given luminosity range actually decreases with increasing luminosity, especially among galaxies with LFIR greater than 10 exp 11 solar luminosities. The close correlation between FIR and radio fluxes of these star-forming galaxies then leads to a picture contrary to earlier theoretical predictions that the sizes of FIR emitting regions should increase with luminosity. We propose that such a steep increase of surface brightness in luminous galaxies is a direct result of a high interstellar medium density and an enhanced radiation density of the heating photons driven by more active star formation.
We conducted a historical cohort study of mortality among 195 astronauts who were exposed to space and medical sources of radiation between 1959 and 1991. Cumulative occupational and medical radiation exposures were obtained from the astronaut radiation exposure history data base. Causes of death were obtained from obligatory death certificates and autopsy reports that were on file in the medical records. A total of 18 deaths occurred during the 32-year follow-up period for which the all-cause standardized mortality ratio (SMR) was 142 (95 percent confidence interval 84 225). There was one cancer death in the buccal cavity and pharynegeal ICD-9 rubric whose occurrence was significantly beyond expectation. Mortality for coronary disease was 59 percent lower than expected (2 deaths; SMR = 41; 95 percent confidence limit 5 147). The crude death rate for 10 occupationally related accidents was 400 deaths per 100,000 person-years, which is an order of magnitude greater than accidental death rates in mining industries. The SMR of 1027 for fatal accidents was significantly beyond expectation (14 deaths; 95 percent confidence limit 561 1723) and was similar to SMRs for accidents among aerial pesticide applications. The 10-year cumulative risk of occupational fatalities based on the exponential, Weibull, Gompertz, and linear-exponential distributions was 10 percent. Mortality from motor vehicle accidents was slightly higher than expected but was not significant (1 death; SMR = 145; 95 percent confidence limit 2 808). Radiation exposures from medical procedures accounted for a majority of cumulative dose when compared with space radiation exposures. The results of the study do not confirm the impression that astronauts are at increased risk of cancer, but this does not obviate the need for further study. Overall, it was found that astronauts are at a health disadvantage as a result of catastrophic accidents.
The specific attraction and, in large part, the significance of solar magnetograms lie in the fact that they give the most important data on the electric currents and the nonpotentiality of active regions. Using the vector magnetograms from the Marshall Space Flight Center (MSFC), we employ a unique technique in the area of data analysis for resolving the 180 deg ambiguity in order to calculate the spatial structure of the vertical electric current density. The 180 deg ambiguity is resolved by applying concepts from the nonlinear multivariable optimization theory. The technique is shown to be of particular importance in very nonpotential active regions. The characterization of the vertical electric current density for a set of vector magnetograms using this method then gives the spatial scale, locations, and magnitude of these current systems. The method, which employs an intermediate parametric function which covers the magnetogram and which defines the local `preferred' direction, minimizes a specific functional of the observed transverse magnetic field. The specific functional that is successful is the integral of the square of the vertical current density. We find that the vertical electric current densities have common characteristics for the extended bipolar (beta) (gamma) (delta)-regions studied. The largest current systems have j(sub z)'s which maximizes around 30 mA/sq m and have a linear decreasing distribution to a diameter of 30 Mn.
This monograph provides an extensive list of formulas for airfoil polynomials. These polynomials provide convenient expansion functions for the description of the downwash and pressure distributions of linear theory for airfoils in both steady and unsteady subsonic flow.
The solar-particle event (SPE) of September 29, 1989, produced an iron-rich spectrum with energies approaching 1 A GeV with an approximate spectral slope parameter of 2.5. These high charge and energy (HZE) ions challenge conventional methods of shield design and assessment of astronaut risks. In the past, shield design and risk assessment have relied on proton shielding codes and biological response models derived from X-ray and neutron exposure data. Because the HZE spectra decline rapidly with energy and HZE attenuation in materials is limited by their penetration power, details of the mass distributions about the sensitive tissues (shielding materials and the astronaut's body) are important determining factors of the exposure levels and distributions of linear energy transfer. Local tissue environments during the SPE of September 29, 1989, with its f= components are examined to analyze the importance of these ions to human SPE exposure. Typical space suit and lightly shielded structures leave significant contributions from HZE components to certain critical body tissues and have important implications on the models for risk assessment. A heavily shielded equipment room of a space vehicle or habitat requires knowledge of the breakup of these ions into lighter components, including neutrons, for shield design specifications.
Magnetic suspensions (MS) satisfy the long life and low loss conditions demanded by satellite and ISS based flywheels used for Energy Storage and Attitude Control (ACESE) service. This paper summarizes the development of a novel MS that improves reliability via fault tolerant operation. Specifically, flux coupling between poles of a homopolar magnetic bearing is shown to deliver desired forces even after termination of coil currents to a subset of failed poles . Linear, coordinate decoupled force-voltage relations are also maintained before and after failure by bias linearization. Current distribution matrices (CDM) which adjust the currents and fluxes following a pole set failure are determined for many faulted pole combinations. The CDM s and the system responses are obtained utilizing 1D magnetic circuit models with fringe and leakage factors derived from detailed, 3D, finite element field models. Reliability results are presented vs. detection/correction delay time and individual power amplifier reliability for 4, 6, and 7 pole configurations. Reliability is shown for two success criteria, i.e. (a) no catcher bearing contact following pole failures and (b) re-levitation off of the catcher bearings following pole failures. An advantage of the method presented over other redundant operation approaches is a significantly reduced requirement for backup hardware such as additional actuators or power amplifiers.
The specific attraction and, in large part, the significance of solar vector magnetograms lie in the fact that they give the most important data on the electric currents and the nonpotentiality of active regions. Using the vector magnetograms from the Marshall Space Flight Center (MSFC), we employ a unique technique in the area of data analysis for resolving the 180 degree ambiguity in order to calculate the spatial structure of the vertical electric current density. The 180 degree ambiguity is resolved by applying concepts from the nonlinear multivariable optimization theory. The technique is shown to be of particular importance in very nonpotential active regions. The characterization of the vertical electric current density for a set of vector magnetograms using this method then gives the spatial scale, locations, and magnitude of these current systems. The method, which employs an intermediate parametric function which covers the magnetogram and which defines the local "preferred" direction, minimizes a specific functional of the observed transverse magnetic field. The specific functional that is successful is the integral of the square of the vertical current density. We find that the vertical electric current densities have common characteristics for the extended bipolar beta gamma delta-regions studied. The largest current systems have j(sub z)'s which maximizes around 30 mA per square meter and have a linear decreasing distribution to a diameter of 30 Mm.
The Refined Zigzag Theory (RZT) for homogeneous, laminated composite, and sandwich plates is presented from a multi-scale formalism starting with the inplane displacement field expressed as a superposition of coarse and fine contributions. The coarse kinematic field is that of first-order shear-deformation theory, whereas the fine kinematic field has a piecewise-linear zigzag distribution through the thickness. The condition of limiting homogeneity of transverse-shear properties is proposed and yields four distinct sets of zigzag functions. By examining elastostatic solutions for highly heterogeneous sandwich plates, the best-performing zigzag functions are identified. The RZT predictive capabilities to model homogeneous and highly heterogeneous sandwich plates are critically assessed, demonstrating its superior efficiency, accuracy ; and a wide range of applicability. The present theory, which is derived from the virtual work principle, is well-suited for developing computationally efficient CO-continuous finite elements, and is thus appropriate for the analysis and design of high-performance load-bearing aerospace structures.
Wind spacecraft measurements are analyzed to obtain a current sheet (CS) normal width dcs-distribution of 3374 confirmed magnetic reconnection exhausts in the ecliptic plane of the solar wind at 1 AU. The dcs-distribution displays a nearly exponential decay from a peak at dcs=25 di to a median at dcs=85 di and a 95th percentile at dcs=905 di with a maximum exhaust width at dcs=8077 di. A magnetic field ϴ-rotation angle distribution increases linearly from a relatively few high-shear events toward a broad peak at 35o<ϴ<65o. The azimuthal φ-angles of the CS normal directions of 430 thick dcs≥500 di exhausts are consistent with a dominant Parker-spiral magnetic field and a CS normal along the ortho-Parker direction. The CS normal orientations of 370 kinetic-scale dcs<25 di exhausts are isotropic in contrast, and likely associated with Alfvénic solar wind turbulence. We propose that the alignment of exhaust normal directions from narrow dcs~15-25 di widths to well beyond dcs~500 di with an ortho-Parker azimuthal direction of a large-scale HCS is a consequence of CS bifurcation and turbulence within the HCS exhaust that may trigger reconnection of the adjacent pair of bifurcated CSs. The proposed HCS-31 avalanche scenario suggests that the underlying large-scale parent HCS closer to the Sun evolves with heliocentric distance to fracture into many, more or less aligned, secondary current sheets due to reconnection. A few wide exhaust-associated HCS-like CSs could represent a population of HCSs that failed to reconnect as frequently between the Sun and 1 AU as other HCSs.
This study presents a method for derivation of effective properties at the interface and in-depth of porous materials. The method defines a Representative Elementary Volume (REV) and applies filtering techniques to computer effective properties such as porosity and flow quantities, such as velocity and pressure. The script, developed to process the data was tested on the VTK type files that contain the mesh information and the flow solution. The method allows to choose between two types of filters, such as cellular and top-hat and define the size of the REV and number of samples along the domain. Extraction of the REV from the domain is performed to exact boundaries requested for the user. This is done using a triangulation technique and cutting through the cells to comply to the requested boundaries of the volume. The method can be applied to both structured and unstructured meshes. Filtering the material porosity and flow quantities involves integration of the numerical data. The algorithm provides three integration methods, such as Riemann sum, Monte Carlo and Quadrature rule to perform the integration. The Monte-Carlo technique permits the use of either uniform or linearly spaced distribution of points. The Quadrature rule is currently applicable to tetrahedral element types. The Monte Carlo and Quadrature rule methods require interpolation of the flow quantities at the sample points. For interpolation, two methods were tested and are readily available, Gaussian interpolation and re-sampling. It has been shown that re-sampling method has better consistency and acceptable accuracy in interpolation of the data. The algorithm was written in Python language and uses a number of modules. The major module besides numpy is PyVista. It is used to process the computational domain, clip the REV and interpolate the data. Quadrature rule integration was performed using a quadpy module. ParaView software was used externally to convert the flow solution to the VTK (or more specifically VTU) format. Integration of ParaView in the same environment with PyVista encountered problems and could not be implemented in this work. The developed algorithm is expected to be applicable to unstructured meshes and more complex porous structures as soon as the data can be passed in VTK type format. With the report is provided Python script for filtering the solution and a Matlab script for simple generation and processing of 2-D and 3-D porous channel geometries. The two scripts don't communicate.