Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linearized distribution”

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

Nonperturbative methods in HZE ion transport

A nonperturbative analytic solution of the high charge and energy (HZE) Green's function is used to implement a computer code for laboratory ion beam transport. The code is established to operate on the Langley Research Center nuclear fragmentation model used in engineering applications. Computational procedures are established to generate linear energy transfer (LET) distributions for a specified ion beam and target for comparison with experimental measurements. The code is highly efficient and compares well with the perturbation approximations.

Wilson, John W.↗

Approximate Green's function methods for HZE transport in multilayered materials

A nonperturbative analytic solution of the high charge and energy (HZE) Green's function is used to implement a computer code for laboratory ion beam transport in multilayered materials. The code is established to operate on the Langley nuclear fragmentation model used in engineering applications. Computational procedures are established to generate linear energy transfer (LET) distributions for a specified ion beam and target for comparison with experimental measurements. The code was found to be highly efficient and compared well with the perturbation approximation.

Wilson, John W.↗

LET spectra of fragmented 600-MeV/A iron beam on aluminum and polyethylene: A preliminary analysis

An iron beam experiment recently conducted at the Lawrence Berkeley Laboratory Bevalac by Benton et al. provides an opportunity for verifying the new Green's function computer code (GRNTRN) and assessing the related nuclear database. The iron beam with 600 MeV/A at extraction traversed a series of beam transport elements, lead foil, and several triggering devices before impacting the target. Of these, only the 2.24 g/sq cm lead foil and target are considered in the transport analysis with an assumed (inferred from calibration) beam energy of 557 MeV/A. A thin layer of CR-39 plastic nuclear track detector (PNTD) was placed in front of the target to monitor the incident flux and a stack of four PNTD's placed behind to measure the linear energy transfer (LET) distribution of the transported beam. Test data are analyzed for three separate targets: 2 g/sq cm aluminum, 5-cm polyethylene, and 8-cm polyethylene. The two-layer GRNTRN results were mapped into the detector response function for comparison with the measured LET spectra. Reasonable agreement is obtained. Future research and analysis can be improved by using a more accurate isotope set or including other important media which significantly alter the beam. The assumption that the fragmentation cross sections are too small for aluminum is inconclusive because the production of fragments by 2 g/sq cm aluminum can be more significantly affected by the fragment contribution originating from the media in front of the target than that currently considered in the analysis.

Shinn, Judy L.↗

HZE beam transport in multilayered materials

A nonperturbative analytic solution of the high charge and energy (HZE) Green's function is used to implement a computer code for laboratory ion beam transport in multiple-layered materials. The code is established to operate on the Langley nuclear fragmentation model used in space engineering applications. Computational procedures are established to generate linear energy transfer (LET) distributions for a specified ion beam and target for comparison with experimental measurement. Comparison with Fe-56 ion with Pb-Al and Pb-(CH2)(x) targets shows reasonable agreement.

Shinn, J. L.↗

Active Sonic Boom Control

A theory and simulation code are developed to study nonsteady sources as means to control sonic booms of supersonic aircraft. A key result is that the source of sonic boom pressure is not confined to the length of the aircraft but occupies an extensive segment of flight path. An aircraft in nonsteady flight functions as a synthetic aperture antenna, generating complex acoustic waves with no simple relation to instantaneous volume or lift distributions. The theory applies linear acoustics to slender nonsteady sources but requires no far field approximation. The solution for pressure contains a term not seen in Whitham's theory for sonic booms of distant supersonic aircraft. The term describes a pressure field that decays algebraically behind the Mach cone and, in the case of steady flight, integrates to a ground load equal to the weight of the aircraft. The algebraic term is separate from those that describe the sonic boom. Two nonsteady source phenomena are evaluated: periodic velocity changes (surge), and periodic longitudinal lift redistribution (slosh). Surge can attenuate a sonic boom and convert it into prolonged weak reverberation, but accelerations needed to produce the phenomenon seem too large for practical use. Slosh may be practical and can alter sonic booms but does not, on the average, result in boom attenuation. The conclusion is that active sonic boom abatement is possible in theory but maybe not practical

Crow, Steven C.↗

Land Surface Temperature Measurements from EOS MODIS Data

We made modifications to the linear kernel bidirectional reflectance distribution function (BRDF) models from Roujean et al. and Wanner et al. that extend the spectral range into the thermal infrared (TIR). With these TIR BRDF models and the IGBP land-cover product, we developed a classification-based emissivity database for the EOS/MODIS land-surface temperature (LST) algorithm and used it in version V2.0 of the MODIS LST code. Two V2.0 LST codes have been delivered to the MODIS SDST, one for the daily L2 and L3 LST products, and another for the 8-day 1km L3 LST product. New TIR thermometers (broadband radiometer with a filter in the 10-13 micron window) and an IR camera have been purchased in order to reduce the uncertainty in LST field measurements due to the temporal and spatial variations in LST. New improvements have been made to the existing TIR spectrometer in order to increase its accuracy to 0.2 C that will be required in the vicarious calibration of the MODIS TIR bands.

Wan, Zhengming↗

Numerical Study of the Generation of Linear Energy Transfer Spectra for Space Radiation Applications

In analyzing charged particle spectra in space due to galactic cosmic rays (GCR) and solar particle events (SPE), the conversion of particle energy spectra into linear energy transfer (LET) distributions is a convenient guide in assessing biologically significant components of these spectra. The mapping of LET to energy is triple valued and can be defined only on open energy subintervals where the derivative of LET with respect to energy is not zero. Presented here is a well-defined numerical procedure which allows for the generation of LET spectra on the open energy subintervals that are integrable in spite of their singular nature. The efficiency and accuracy of the numerical procedures is demonstrated by providing examples of computed differential and integral LET spectra and their equilibrium components for historically large SPEs and 1977 solar minimum GCR environments. Due to the biological significance of tissue, all simulations are done with tissue as the target material.

Badavi, Francis F.↗

Guidance and Control System for a Satellite Constellation

A distributed guidance and control algorithm was developed for a constellation of satellites. The system repositions satellites as required, regulates satellites to desired orbits, and prevents collisions. 1. Optimal methods are used to compute nominal transfers from orbit to orbit. 2. Satellites are regulated to maintain the desired orbits once the transfers are complete. 3. A simulator is used to predict potential collisions or near-misses. 4. Each satellite computes perturbations to its controls so as to increase any unacceptable distances of nearest approach to other objects. a. The avoidance problem is recast in a distributed and locally-linear form to arrive at a tractable solution. b. Plant matrix values are approximated via simulation at each time step. c. The Linear Quadratic Gaussian (LQG) method is used to compute perturbations to the controls that will result in increased miss distances. 5. Once all danger is passed, the satellites return to their original orbits, all the while avoiding each other as above. 6. The delta-Vs are reasonable. The controller begins maneuvers as soon as practical to minimize delta-V. 7. Despite the inclusion of trajectory simulations within the control loop, the algorithm is sufficiently fast for available satellite computer hardware. 8. The required measurement accuracies are within the capabilities of modern inertial measurement devices and modern positioning devices.

Bryson, Jonathan Lamar↗

Multiple Velocity Profile Measurements in Hypersonic Flows Using Sequentially-Imaged Fluorescence Tagging

Nitric-oxide planar laser-induced fluorescence (NO PLIF) was used to perform velocity measurements in hypersonic flows by generating multiple tagged lines which fluoresce as they convect downstream. For each laser pulse, a single interline, progressive scan intensified CCD (charge-coupled device) camera was used to obtain two sequential images of the NO molecules that had been tagged by the laser. The CCD configuration allowed for sub-microsecond acquisition of both images, resulting in sub-microsecond temporal resolution as well as sub-mm spatial resolution (0.5-mm horizontal, 0.7-mm vertical). Determination of axial velocity was made by application of a cross-correlation analysis of the horizontal shift of individual tagged lines. A numerical study of measured velocity error due to a uniform and linearly-varying collisional rate distribution was performed. Quantification of systematic errors, the contribution of gating/exposure duration errors, and the influence of collision rate on temporal uncertainty were made. Quantification of the spatial uncertainty depended upon the signal-to-noise ratio of the acquired profiles. This velocity measurement technique has been demonstrated for two hypersonic flow experiments: (1) a reaction control system (RCS) jet on an Orion Crew Exploration Vehicle (CEV) wind tunnel model and (2) a 10-degree half-angle wedge containing a 2-mm tall, 4-mm wide cylindrical boundary layer trip. The experiments were performed at the NASA Langley Research Center's 31-Inch Mach 10 Air Tunnel.

Bathel, Brett F.↗

Factorization and the synthesis of optimal feedback gains for distributed parameter systems

An approach based on Volterra factorization leads to a new methodology for the analysis and synthesis of the optimal feedback gain in the finite-time linear quadratic control problem for distributed parameter systems. The approach circumvents the need for solving and analyzing Riccati equations and provides a more transparent connection between the system dynamics and the optimal gain. The general results are further extended and specialized for the case where the underlying state is characterized by autonomous differential-delay dynamics. Numerical examples are given to illustrate the second-order convergence rate that is derived for an approximation scheme for the optimal feedback gain in the differential-delay problem.

Milman, Mark H.↗

Linear and nonlinear properties of the ULF waves driven by ring-beam distribution functions

The problem of the exitation of obliquely propagating magnetosonic waves which can steepen up (also known as shocklets) is considered. Shocklets have been observed upstream of the Earth's bow shock and at comets Giacobini-Zinner and Grigg-Skjellerup. Linear theory as well as two-dimensional (2-D) hybrid (fluid electrons, particle ions) simulations are used to determine the properties of waves generated by ring-beam velocity distributions in great detail. The effects of both proton and oxygen ring-beams are considered. The study of instabilities excited by a proton ring-beam is relevant to the region upstream of the Earth's bow shock, whereas the oxygen ring-beam corresponds to cometary ions picked up by the solar wind. Linear theory has shown that for a ring-beam, four instabilities are found, one on the nonresonant mode, one on the Alfven mode, and two along the magnetosonic/whistler branch. The relative growth rate of these instabilities is a sensitive function of parameters. Although one of the magnetosonic instabilities has maximum growth along the magnetic field, the other has maximum growth in oblique directions. We have studied the competition of these instabilities in the nonlinear regime using 2-D simulations. As in the linear limit, the nonlinear results are a function of beam density and distribution function. By performing the simulations as both initial value and driven systems, we have found that the outcome of the simulations can vary, suggesting that the latter type simulations is needed to address the observations. A general conclusion of the simulation results is that field-aligned beams do not result in the formation of shocklets, whereas ring-beam distributions can.

Killen, K.↗

A PV view of the zonal mean distribution of temperature and wind in the extratropical troposphere

The dependence of the temperature and wind distribution of the zonal mean flow in the extratropical troposphere on the gradient of pontential vorticity along isentropes is examined. The extratropics here refer to the region outside the Hadley circulation. Of particular interest is whether the distribution of temperature and wind corresponding to a constant potential vorticity (PV) along isentropes resembles the observed, and the implications of PV homogenization along isentropes for the role of the tropics. With the assumption that PV is homogenized along isentropes, it is found that the temperature distribution in the extratropical troposphere may be determined by a linear, first-order partial differential equation. When the observed surface temperature distribution and tropical lapse rate are used as the boundary conditions, the solution of the equation is close to the observed temperature distribution except in the upper troposphere adjacent to the Hadley circulation, where the troposphere with no PV gradient is considerably colder. Consequently, the jet is also stronger. It is also found that the meridional distribution of the balanced zonal wind is very sensitive to the meridional distribution of the tropopause temperature. The result may suggest that the requirement of the global momentum balance has no practical role in determining the extratropical temperature distribution. The authors further investigated the sensitivity of the extratropical troposphere with constant PV along isentropes to changes in conditions at the tropical boundary (the edge of the Hadley circulation). It is found that the temperature and wind distributions in the extratropical troposphere are sensitive to the vertical distribution of PV at the tropical boundary. With a surface distribution of temperature that decreases linearly with latitude, the jet maximum occurs at the tropical boundary and moves with it. The overall pattern of wind distribution is not sensitive to the change of the position of the tropical boundary. Finally, the temperature and wind distributions of an extratropical troposphere with a finite PV gradient are calculated. It is found that the larger the isentropic PV gradient, the warmer the troposphere and the weaker the jet.

Sun, De-Zheng↗

Thin-plate spline quadrature of geodetic integrals

Thin-plate spline functions (known for their flexibility and fidelity in representing experimental data) are especially well-suited for the numerical integration of geodetic integrals in the area where the integration is most sensitive to the data, i.e., in the immediate vicinity of the evaluation point. Spline quadrature rules are derived for the contribution of a circular innermost zone to Stoke's formula, to the formulae of Vening Meinesz, and to the recursively evaluated operator L(n) in the analytical continuation solution of Molodensky's problem. These rules are exact for interpolating thin-plate splines. In cases where the integration data are distributed irregularly, a system of linear equations needs to be solved for the quadrature coefficients. Formulae are given for the terms appearing in these equations. In case the data are regularly distributed, the coefficients may be determined once-and-for-all. Examples are given of some fixed-point rules. With such rules successive evaluation, within a circular disk, of the terms in Molodensky's series becomes relatively easy. The spline quadrature technique presented complements other techniques such as ring integration for intermediate integration zones.

Vangysen, Herman↗

Boundary-Layer Receptivity and Integrated Transition Prediction

The adjoint parabold stability equations (PSE) formulation is used to calculate the boundary layer receptivity to localized surface roughness and suction for compressible boundary layers. Receptivity efficiency functions predicted by the adjoint PSE approach agree well with results based on other nonparallel methods including linearized Navier-Stokes equations for both Tollmien-Schlichting waves and crossflow instability in swept wing boundary layers. The receptivity efficiency function can be regarded as the Green's function to the disturbance amplitude evolution in a nonparallel (growing) boundary layer. Given the Fourier transformed geometry factor distribution along the chordwise direction, the linear disturbance amplitude evolution for a finite size, distributed nonuniformity can be computed by evaluating the integral effects of both disturbance generation and linear amplification. The synergistic approach via the linear adjoint PSE for receptivity and nonlinear PSE for disturbance evolution downstream of the leading edge forms the basis for an integrated transition prediction tool. Eventually, such physics-based, high fidelity prediction methods could simulate the transition process from the disturbance generation through the nonlinear breakdown in a holistic manner.

Chang, Chau-Lyan↗

Distribution function of continuously created newborn and pickup ions in outer cometary exospheres

The time evolution of the distribution function of newborn ions in the solar wind is investigated using a quasi-linear-type diffusion equation. The initial distribution is taken to be a ring beam, which is approximated by delta function in pitch angle and velocity, and it is assumed that the ions are created at a constant rate with a similar distribution. A long-time asymptotic form of ion distribution is obtained, which is a mixture of newborn ions and ions generated throughout the entire process. It is shown that the time asymptotic distribution function exists even in the presence of a continuous ionization process. The stability of the long-time asymptotic distribution was examined for the case of parallel propagation, and the results show that the distribution function can be unstable to low-frequency hydromagnetic waves. The results of the analysis were found to agree with recent satellite observations.

Gaffey, J. D., Jr.↗

Proton dosimeter design for distributed body organs

A simple dosimeter design has been developed by NASA to monitor the space proton dose (in rads) to a distributed body organ as a linear combination of ion chambers with varying wall thickness. Estimated doses are given for ion chambers of thickness 2, 3, 4, and 5 g/sq cm. The analytical equation used to calculate the dose distribution factor is also given.

Wilson, J. W.↗

Stability of electron distributions within the earth's bow shock

The present examination of the linear Vlasov stability of a class of electron velocity distributions modeling those observed within the earth's bow shock is restricted to electrostatic waves propagating parallel to the ambient magnetic field B. Two instabilities are identified as driven by free energy in the direction parallel to B: an ion acoustic wave with real frequency below the ion plasma frequency, and an electron acoustic wave whose real frequency is several times the ion plasma frequency. Unstable wave characteristics are in accord with the trend of the observed electrostatic waves toward polarization parallel to the magnetic field. The instabilities identified may contribute to electron dissipation in collisionless shocks.

Thomsen, M. F.↗