Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Approximation methods”

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 1,369 records · Page 76

Thin shell formation in radiative shocks. 1: Supernova remnants in low-density media

This paper explores the onset of thin-shell formation in interstellar shocks associated with supernova explosions. We outline a simple but useful scheme that indicates the time at which thin shell formation begins for supernova remnants (SNRs) evolving in a range of interstellar environments, extending the previous analytical models to arbitrary power-law density media. The result depends on the gas cooling properties and the shock velocity and radius. This is then applied to the specific case of SNRs in low-density media. The procedure for defining the time for the onset of shell formation, t(sub sf), equates the value of the adiabat, kappa = p/rho(exp gamma), to zero using the known time dependence of the shock radius and velocity. For the case of a power-law density ambient medium of the form rho(r) = Br(exp -omega), it is found that shell formation can be prevented when the ambient density drops faster than a critical rate. For a cooling function of the form Lambda = Lambda(sub 0) tau(exp beta), with beta = -0.5 (appropriate for line cooling), shell formation never occurs for omega greater than or equal to 9/5. The shell formation time is then computed for spherical shocks in a power-law density medium. For omega = 0, the onset of shell formation is found to be at t(sub sf) approx. equal to 2.87 x 10(exp 4) E51(exp 3/14) n(sub 0 exp -4.7) yr, which agrees well with previous estimates derived by other means. We compare the analytical shell formation time with the results of detailed numerical models for omega = 0 and three different ambient densities and find good agreement. The extension of the criterion for the onset of thin shell formation using the ratio of cooling to swept-up column density is also described. This method provides a useful approximation for cases when the exact solution is not known.

Franco, Jose↗

The Orion OB1 association. 1: Stellar content

Walfraven photometry of established and probable members of the Orion OB1 association is presented. Effective temperature, surface gravity, luminosity and mass are derived for all stars, using atmosphere model by Kurucz (1979). Absolute magnitudes are calculated using the Straizys and Kuriliene (1981) tables. Distance moduli and visual extinctions are determined. A comparison of the visual extinctions to IRAS 100 micrometers data shows that the near edge of the Orion A and B clouds lies at a distance of approximately 320 pc, while the far edge is at approximately 500 pc. A method for deriving the ages of the subgroups by comparing theoretical isochrones to the observations in the log g, log T(sub eff) plane is presented. The derived ages suggest, contrary to earlier studies, that subgroup 1b is younger than 1c, which can possibly be explained by past geometries of the system of stars and gas. The initial mass function for Orion OB1 is derived with the aid of the Kolmogorov-Smirnoff test. Through extensive simulations, we show that it is very difficult to derive accurately the Initial Mass Function (IMF) from the available data. To within somewhat weak limits the IMF is found to be of the form xi(log M) = AM(exp -1.7 +/- 0.2) for all subgroups. The energy output of the subgroups in the form of stellar winds and supernovae is calculated and compared to the observed size and expansion velocity of the Orion-Eridanus bubble. It is shown that the energy output of the association can account for the morphology and kinematics of the interstellar medium (ISM).

Brown, A. G. A.↗

Spectral Density of Laser Beam Scintillation in Wind Turbulence: Theory - Part 1

The temporal spectral density of the log-amplitude scintillation of a laser beam wave due to a spatially dependent vector-valued crosswind (deterministic as well as random) is evaluated. The path weighting functions for normalized spectral moments are derived, and offer a potential new technique for estimating the wind velocity profile. The Tatarskii-Klyatskin stochastic propagation equation for the Markov turbulence model is used with the solution approximated by the Rytov method. The Taylor 'frozen-in' hypothesis is assumed for the dependence of the refractive index on the wind velocity, and the Kolmogorov spectral density is used for the refractive index field.

Balakrishnan, A. V.↗

Orbital Debris Impact Damage to Reusable Launch Vehicles

In an effort by the National Aeronautics and Space Administration (NASA), hypervelocity impact tests were performed on thermal protection systems (TPS) applied on the external surfaces of reusable launch vehicles (RLV) to determine the potential damage from orbital debris impacts. Three TPS types were tested, bonded to composite structures representing RLV fuel tank walls. The three heat shield materials tested were Alumina-Enhanced Thermal Barrier-12 (AETB-12), Flexible Reusable Surface Insulation (FRSI), and Advanced Flexible Reusable Surface Insulation (AFRSI). Using this test data, predictor equations were developed for the entry hole diameters in the three TPS materials, with correlation coefficients ranging from 0.69 to 0.86. Possible methods are proposed for approximating damage occurring at expected orbital impact velocities higher than tested, with references to other published work.

Robinson, Jennifer H.↗

Method Improvements in Thermal Analysis of Mach 10 Leading Edges

Several improvements have recently been made in the thermal analysis methods for leading edges of a hypersonic vehicle. The leading edges of this vehicle undergo exceptionally high heat loads that incorporate extreme spatial gradients as well as severe transients. Due to the varying flight conditions, complex geometry, and need for thermal loads at many points along the trajectory, full computational fluid dynamics (CFD) analysis of the aeroheating; loads is not feasible. Thus, engineering methods must be used to determine the aeroheating on the vehicle surfaces, and that must be utilized in the thermal analysis. Over the last year, the thermal analysis of a hypersonic vehicle has been enhanced in several ways. Two different engineering codes are used to predict aeroheating loads: one over the curve near the stagnation point, and the other on flat surfaces downstream of the leading edge. These two are matched together at the intersection point using a method that allows closer approximation of CFD results. User-developed FORTRAN, which is part of the thermal solver PATRAN Thermal, is used to accomplish this. The customizable FORTRAN code also allows use of many different time- and space-dependent factors, interpolation of the heat load in time and space, and inclusion of both highly swept and unswept grid structures. This FORTRAN is available to other PATRAN users who may want to accomplish a similar objective in analysis. Flux, rather than convective coefficient, is used to define heat loads, which allows more accurate analysis as well as better application of margins. Improvements have also been made in more efficient utilization of imported CAD geometry, by creating faces on solids to facilitate load application.

Amundsen, Ruth M.↗

Calculation of Scattering Amplitude Without Partial Analysis: Inclusion of Exchange - II

There was a method for calculating the whole scattering amplitude, f(Omega(sub k)), directly. The idea was to calculate the complete wave function Psi numerically, and use it in an integral expression for f, which can be reduced to a 2 dimensional quadrature. The original application was for e-H scattering without exchange. There the Schrodinger reduces a 2-d partial differential equation (pde), which was solved using the finite element method (FEM). Here we extend the method to the exchange approximation. The S.E. can be reduced to a pair of coupled pde's, which are again solved by the FEM. The formal expression for f(Omega(sub k)) consists two integrals, f+/- = f(sub d) +/- f(sub e); f(sub d) is formally the same integral as the no-exchange f. We have also succeeded in reducing f(sub e) to a 2-d integral. Results will be presented at the meeting.

Temkin, Aaron↗

Rooting the archaebacterial tree: the pivotal role of Thermococcus celer in archaebacterial evolution

The sequence of the 16S ribosomal RNA gene from the archaebacterium Thermococcus celer shows the organism to be related to the methanogenic archaebacteria rather than to its phenotypic counterparts, the extremely thermophilic archaebacteria. This conclusion turns on the position of the root of the archaebacterial phylogenetic tree, however. The problems encountered in rooting this tree are analyzed in detail. Under conditions that suppress evolutionary noise both the parsimony and evolutionary distance methods yield a root location (using a number of eubacterial or eukaryotic outgroup sequences) that is consistent with that determined by an "internal rooting" method, based upon an (approximate) determination of relative evolutionary rates.

Non-NASA Center↗

An Experimental Investigation of Sting-Support Effects on Drag and a Comparison with Jet Effects at Transonic Speeds

This paper presents the results of an investigation of sting-support interference on afterbody drag at transonic speeds. Stings with varying diameter, cone angle, and cylindrical length were tested at the rea r of a model with various afterbody shapes. The data were obtained at an angle of attack of O deg. and at Mach numbers from 0.80 to 1.10. It was found that, in general, the addition of a sting caused a drag reduction. A method is presented whereby approximate sting-interferen ce corrections may be made to models with afterbodies and sting suppo rts of similar size and scale to those of this paper provided the bou ndary layer is turbulent at the model base and the Reynolds numbers a re of the same order of magnitude. Reynolds number of the tests prese nted varied from 15.0 x 10 (exp 6) to 17.4 x 10 (exp 6) based on body length. Sting effects from this investigation are compared with data of jet effects on the same afterbodies. The results of this comparis on indicate that for the more gradually contoured afterbodies, a stin g shape can be found which will duplicate the jet effects, but that f or blunt afterbodies no solid sting shape will duplicate the jet effe cts.

Cahn, Maurice S.↗

Multiparallel Three-Dimensional Optical Microscopy

Multiparallel three-dimensional optical microscopy is a method of forming an approximate three-dimensional image of a microscope sample as a collection of images from different depths through the sample. The imaging apparatus includes a single microscope plus an assembly of beam splitters and mirrors that divide the output of the microscope into multiple channels. An imaging array of photodetectors in each channel is located at a different distance along the optical path from the microscope, corresponding to a focal plane at a different depth within the sample. The optical path leading to each photodetector array also includes lenses to compensate for the variation of magnification with distance so that the images ultimately formed on all the photodetector arrays are of the same magnification. The use of optical components common to multiple channels in a simple geometry makes it possible to obtain high light-transmission efficiency with an optically and mechanically simple assembly. In addition, because images can be read out simultaneously from all the photodetector arrays, the apparatus can support three-dimensional imaging at a high scanning rate.

Nguyen, Lam K.↗

Blocking Losses With a Photon Counter

It was not known how to assess accurately losses in a communications link due to photodetector blocking, a phenomenon wherein a detector is rendered inactive for a short time after the detection of a photon. When used to detect a communications signal, blocking leads to losses relative to an ideal detector, which may be measured as a reduction in the communications rate for a given received signal power, or an increase in the signal power required to support the same communications rate. This work involved characterizing blocking losses for single detectors and arrays of detectors. Blocking may be mitigated by spreading the signal intensity over an array of detectors, reducing the count rate on any one detector. A simple approximation was made to the blocking loss as a function of the probability that a detector is unblocked at a given time, essentially treating the blocking probability as a scaling of the detection efficiency. An exact statistical characterization was derived for a single detector, and an approximation for multiple detectors. This allowed derivation of several accurate approximations to the loss. Methods were also derived to account for a rise time in recovery, and non-uniform illumination due to diffraction and atmospheric distortion of the phase front. It was assumed that the communications signal is intensity modulated and received by an array of photon-counting photodetectors. For the purpose of this analysis, it was assumed that the detectors are ideal, in that they produce a signal that allows one to reproduce the arrival times of electrons, produced either as photoelectrons or from dark noise, exactly. For single detectors, the performance of the maximum-likelihood (ML) receiver in blocking is illustrated, as well as a maximum-count (MC) receiver, that, when receiving a pulse-position-modulated (PPM) signal, selects the symbol corresponding to the slot with the largest electron count. Whereas the MC receiver saturates at high count rates, the ML receiver may not. The loss in capacity, symbol-error-rate (SER), and count-rate were numerically computed. It was shown that the capacity and symbol-error-rate losses track, whereas the count-rate loss does not generally reflect the SER or capacity loss, as the slot-statistics at the detector output are no longer Poisson. It is also shown that the MC receiver loss may be accurately predicted for dead times on the order of a slot.

Moision, Burce E.↗

Line Interference Effects Using a Refined Robert-Bonamy Formalism: the Test Case of the Isotropic Raman Spectra of Autoperturbed N2

A symmetrized version of the recently developed refined Robert-Bonamy formalism [Q. Ma, C. Boulet, and R. H. Tipping, J. Chem. Phys. 139, 034305 (2013)] is proposed. This model takes into account line coupling effects and hence allows the calculation of the off-diagonal elements of the relaxation matrix, without neglecting the rotational structure of the perturbing molecule. The formalism is applied to the isotropic Raman spectra of autoperturbed N2 for which a benchmark quantum relaxation matrix has recently been proposed. The consequences of the classical path approximation are carefully analyzed. Methods correcting for effects of inelasticity are considered. While in the right direction, these corrections appear to be too crude to provide off diagonal elements which would yield, via the sum rule, diagonal elements in good agreement with the quantum results. In order to overcome this difficulty, a re-normalization procedure is applied, which ensures that the off-diagonal elements do lead to the exact quantum diagonal elements. The agreement between the (re-normalized) semi-classical and quantum relaxation matrices is excellent, at least for the Raman spectra of N2, opening the way to the analysis of more complex molecular systems.

nitrogen↗

Rapid State Space Modeling Tool for Rectangular Wing Aeroservoelastic Studies

This paper introduces a modeling and simulation tool for aeroservoelastic analysis of rectangular wings with trailing-edge control surfaces. The inputs to the code are planform design parameters such as wing span, aspect ratio, and number of control surfaces. Using this information, the generalized forces are computed using the doublet-lattice method. Using Roger's approximation, a rational function approximation is computed. The output, computed in a few seconds, is a state space aeroservoelastic model which can be used for analysis and control design. The tool is fully parameterized with default information so there is little required interaction with the model developer. All parameters can be easily modified if desired. The focus of this paper is on tool presentation, verification, and validation. These processes are carried out in stages throughout the paper. The rational function approximation is verified against computed generalized forces for a plate model. A model composed of finite element plates is compared to a modal analysis from commercial software and an independently conducted experimental ground vibration test analysis. Aeroservoelastic analysis is the ultimate goal of this tool, therefore, the flutter speed and frequency for a clamped plate are computed using damping-versus-velocity and frequency-versus-velocity analysis. The computational results are compared to a previously published computational analysis and wind-tunnel results for the same structure. A case study of a generic wing model with a single control surface is presented. Verification of the state space model is presented in comparison to damping-versus-velocity and frequency-versus-velocity analysis, including the analysis of the model in response to a 1-cos gust.

Modeling and Simulation↗

Rapid State Space Modeling Tool for Rectangular Wing Aeroservoelastic Studies

This report introduces a modeling and simulation tool for aeroservoelastic analysis of rectangular wings with trailing-edge control surfaces. The inputs to the code are planform design parameters such as wing span, aspect ratio, and number of control surfaces. Using this information, the generalized forces are computed using the doublet-lattice method. Using Roger's approximation, a rational function approximation is computed. The output, computed in a few seconds, is a state space aeroservoelastic model which can be used for analysis and control design. The tool is fully parameterized with default information so there is little required interaction with the model developer. All parameters can be easily modified if desired. The focus of this report is on tool presentation, verification, and validation. These processes are carried out in stages throughout the report. The rational function approximation is verified against computed generalized forces for a plate model. A model composed of finite element plates is compared to a modal analysis from commercial software and an independently conducted experimental ground vibration test analysis. Aeroservoelastic analysis is the ultimate goal of this tool, therefore, the flutter speed and frequency for a clamped plate are computed using damping-versus-velocity and frequency-versus-velocity analysis. The computational results are compared to a previously published computational analysis and wind-tunnel results for the same structure. A case study of a generic wing model with a single control surface is presented. Verification of the state space model is presented in comparison to damping-versus-velocity and frequency-versus-velocity analysis, including the analysis of the model in response to a 1-cos gust.

Verification and Validation↗

Development of an Integrated Nonlinear Aeroservoelastic Flight Dynamic Model of the NASA Generic Transport Model

This paper describes a recent development of an integrated fully coupled aeroservoelastic flight dynamic model of the NASA Generic Transport Model (GTM). The integrated model couples nonlinear flight dynamics to a nonlinear aeroelastic model of the GTM. The nonlinearity includes the coupling of the rigid-body aircraft states in the partial derivatives of the aeroelastic angle of attack. Aeroservoelastic modeling of the control surfaces which are modeled by the Variable Camber Continuous Trailing Edge Flap is also conducted. The R.T. Jones' method is implemented to approximate unsteady aerodynamics. Simulations of the GTM are conducted with simulated continuous and discrete gust loads..

Aeroservoelasticit↗

Change of Inertia Tensor Due to a Severed Radial Boom for Spinning Spacecraft

Many spinning spacecraft have long, flexible, radial booms to carry science instrumentation. These radial booms often have low mass but contribute significantly to the spacecraft moment of inertia due to their length. There are historical cases where radial booms have been severed or have failed to deploy. This paper presents models for the center of mass (CM) and inertia tensor that account for variable boom geometry and investigates how the CM and inertia tensor change when a radial boom is severed.The CM and inertia tensor models presented here will be included in the Attitude Ground System (AGS) for the Magnetospheric Multiscale (MMS) mission. This work prepares the AGS to provide uninterrupted support in the event of a radial boom anomaly. These models will improve the AGS computations for spin-axis precession prediction, Kalman filter propagation for the definitive attitude, and mass property generation needed for the onboard control system. As an additional application, a method is developed for approximating the location on the boom where the break occurred based on the new models and readily observable attitude parameters.

Inertia↗

Change of Inertia Tensor Due to a Severed Radial Boom for Spinning Spacecraft

Many spinning spacecraft have long, flexible, radial booms to carry science instrumentation. These radial booms often have low mass but contribute significantly to the spacecraft moment of inertia due to their length. There are historical cases where radial booms have been severed or have failed to deploy. This paper presents models for the center of mass (CM) and inertia tensor that account for variable boom geometry and investigates how the CM and inertia tensor change when a radial boom is severed.The CM and inertia tensor models presented here will be included in the Attitude Ground System (AGS) for the Magnetospheric Multiscale (MMS) mission. This work prepares the AGS to provide uninterrupted support in the event of a radial boom anomaly. These models will improve the AGS computations for spin-axis precession prediction, Kalman filter propagation for the definitive attitude, and mass property generation needed for the onboard control system. As an additional application, a method is developed for approximating the location on the boom where the break occurred based on the new models and readily observable attitude parameters.

Sedlak, Joseph E.↗

Oro-Nasal Mask Versus Two-Way Non-Rebreathing Valves for Maximal Aerobic Capacity Testing in Astronauts

Astronauts complete maximal aerobic capacity (VO2pk) testing as part of their annual fitness assessment (AFA) as well as several times once assigned to an International Space Station mission. Historically, the 2-Way T-Shape Non-Rebreathing valve with a mouthpiece and nose clip (“Mouthpiece”) has been used in these tests. The testing procedure was updated to use the oro-nasal mask (“Mask”) for the AFA starting in June 2017. Astronauts who used the mask during their AFA requested it be certified to be used for all mission associated tests. Considering the criticality of the data and the schedule constraints of astronauts, it is imperative that the requested hardware change provide data with equivalent reliability and repeatability as provided by the mouthpiece. PURPOSE: To assess the reliability and validity of mask vs. mouthpiece by comparing submaximal and VO2pkdata within subjects (approximately 1 year apart). METHODS: Each of 17 active astronauts completed a VO (sub 2pk) test with the mouthpiece (first) and the mask (second) for their AFA. The VO (sub 2pk) test was conducted on a cycle ergometer with a metabolic cart. The nominal protocol started with a 3-minute warm-up at 50 Watts (W) and increased 25W every minute until volitional fatigue (Light: 45W start; 15W increase). The VO (sub 2pk)s were compared between tests and the expected day-to-day variation (plus or minus 5 percent) was used as the threshold for determining agreement between tests. Submaximal values were plotted and evaluated visually for deviations between mask and mouthpiece. RESULTS: VO (sub 2pk) values were more than 5 percent different, despite similar test times, between mouthpiece and mask in 6 of 17 comparisons, 3 of which were higher with the mask (9.0 plus or minus 5.9 percent) while 3 were lower (minus10.8 plus or minus 2.0 percent) with the mask. The submaximal data did not indicate a leak in either apparatus during these tests. An Astronaut Strength & Conditioning Rehabilitation specialist confirmed that the measured differences in VO (sub 2pk) of these 6 astronauts was consistent with observed changes in exercise habits during the year that separated the two tests. CONCLUSION: After being presented with the results of this data mining effort the mask was accepted for use in all tests, excepting that, if a leak is detected without resolve, the test will be repeated (if schedule allows) and remaining tests will be completed with the mouthpiece.

Buxton, R.↗

Analysis of Arc Welding Process in Space

This work is motivated by NASA plans to conduct welding experiments on the ISS. It is expected that deployment of welding and additive manufacturing technologies in the space environment has the potential to revolutionise how orbiting platforms are designed, manufactured, and assembled. However, the structure, composition and quality of a weld is extremely dependent on the environment and can be difficult to control in space. Shielding gases would also be tough to manage as gases behave differently in zero gravity and airless environments. Additional points of concern are related to the spatter and sparks dynamics in space. Therefore, there is a need for a more basic understanding of welding processes by computational modelling. To provide such an insight we developed state of art models of the ARC torch, droplet detachment and transfer, and the melt pool build up using magnetohydrodynamics approximation and level set method for two-phase liquid metal/-gas flow modelling. Two finite element models were built using 2D axisymmetric geometry in COMSOL Multiphysics®: (i) stationary model of the ARC torch and (ii) dynamical model of droplet detachment and transfer. Both models demonstrate reasonable agreement with earlier experimental observations and high sensitivity to the temperature dependence of the thermophysical parameters on the system materials. The models were used to provide physical insight into ARC welding in various environments and geometries. It was shown, in particular, that there is a significant probability of gas bubbles trapping in the meltpool and a possibility of unbounded wondering of sparks in the welding chamber in zero gravity. In addition, we estimated metal evaporation rate that may be hazardous in the confined environment of the ISS. These issues raise concern of quality of the weld and safety for ISS applications.

ARC welding↗