Search NASA⌕ Search

SEARCH · Search NASA

Results for “calculation method”

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,387 records · Page 77

Heating of the solar corona by the resonant absorption of Alfven waves

An improved method for calculating the resonance absorption heating rate is discussed and the results are compared with observations in the solar corona. To accomplish this, the wave equation for a dissipative, compressible plasma is derived from the linearized magnetohydrodynamic equations for a plasma with transverse Alfven speed gradients. For parameters representative of the solar corona, it is found that a two-scale description of the wave motion is appropriate. The large-scale motion, which can be approximated as nearly ideal, has a scale which is on the order of the width of the loop. The small-scale wave, however, has a transverse scale much smaller than the width of the loop, with a width of about 0.3-250 km, and is highly dissipative. These two wave motions are coupled in a narrow resonance region in the loop where the global wave frequency equals the local Alfven wave frequency. Formally, this coupling comes about from using the method of matched asymptotic expansions to match the inner and outer (small and large scale) solutions. The resultant heating rate can be calculated from either of these solutions. A formula derived using the outer (ideal) solution is presented, and shown to be consistent with observations of heating and line broadening in the solar corona.

Davila, Joseph M.↗

Strength of adhesive-bonded hybrid structures

Structures prepared from materials with different thermal and mechanical properties by means of fiber-strengthened binders can fail in a number of ways. The present lecture is focused on failures through debonding at the metal or at the fiber-reinforced plastic. A method for calculating the stress distribution in adhesive layers as a function of the load is outlined, and its usefulness in providing insight into the behavior of bonds in hybrid structures is noted. Means of eliminating the unfavorable effects of temperature, humidity, creep and relaxation on the bonds in the manufacture of hybrid structures are examined, along with test methods developed for such structures.

Kirschke, L.↗

Development of a phonon-based sampling method for thermal neutron scattering data

Simulations of reactor systems require access to accurate nuclear data. For many systems, thermal neutron scattering data can have large effects on the eigenvalue and neutron flux distributions. Inelastic thermal neutron scattering can excite or de-excite vibrational, rotational, and translational modes in a material, so thermal scattering evaluations are often obtained by summing over the number of phonons created/destroyed by a scattering event. In recent years, the thermal scattering cross sections and angular distributions have greatly improved in accuracy, but the format in which this data is delivered to simulation codes has remained virtually unchanged. Thermal scattering data is typically either compiled into large tables and sorted by incoming neutron energy, outgoing neutron energy, scattering angle, and material temperature, or represented as cumulative distribution functions of momentum exchange or energy exchange. Either method can be quite memory intensive when fine bins are used. In an effort to decrease the amount of space that processed thermal scattering data requires, an alternate format is proposed. The phonon-based sampling method introduced here can sample the number of phonons excited for each collision, the change in neutron energy, and the scattering angle while avoiding pre-computed angular bins and limiting the amount of data that is dependent on incoming energy. Through this method, the generation and storage of large interpolation tables is avoided, which could have benefits in both memory storage and accuracy. While the initial implementation of this method is slower than current alternatives, it is significantly more resistant to grid coarseness errors and has good potential for improvement. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Lunar physical librations and laser ranging

The analysis of lunar laser ranging data requires very accurate calculations of the lunar physical librations. Libration terms are given which arise from the additive and planetary terms in the lunar theory. The large size of the recently discovered terms due to third degree gravitational harmonics will allow some of these harmonics to be measured by laser ranging to the moon. Numerical integration promises to be an effective method of calculating librations. Comparison of numerical integrations with analytic series indicates that the calculation of the series due to third and fourth degree harmonics is not yet as accurate as the more extensively developed second degree terms.-

Williams, J. G.↗

An Improved Magnetic Field Model for Jupiter's Inner Magnetosphere Using a Microsignature of Amalthea

Observation of a particle absorption microsignature of the Jovian satellite Amalthea during the Pioneer II close flyby of Jupiter on December 3, 1974, has been described by McKibben et al. The microsignature was also observed by the University of Iowa/Pioneer II instrument in the distribution of protons in the kinetic energy range 0.61 less than T(sub p) less than 3.41 MeV but has not been previously reported. The finer time resolution and superior data quality of the latter observations provide a fresh basis for assessing the accuracy with which various published magnetic field models of Jupiter's magnetic field models describe Jupiter's magnetic field for radial distances of the order of or less than several planetary radii. The expected time of occurrence of the minimum of Amalthea's microsignature was calculated for each model and compared to the time of the observed minimum. The discrepancy between these two times was significant for each of the published models, but it could be reduced to zero by interpolating between two of the closely related models. The spherical harmonic coefficients of this combination model are tabulated. Also included is a novel method for calculating the intensity-time profile of the microsignature. In a later paper, the improved model is utilized to analyze the energetic particle measurements in Jupiter's innermost magnetosphere by the Galileo entry probe in December 1995.

Randall, Bruce A.↗

Advanced Modeling Techniques to Study Anthropogenic Influences on Atmospheric Chemical Budgets

This research work is a collaborative effort between research groups at MCNC and the University of North Carolina at Chapel Hill. The overall objective of this research is to improve the level of understanding of the processes that determine the budgets of chemically and radiatively active compounds in the atmosphere through development and application of advanced methods for calculating the chemical change in atmospheric models. The research performed during the second year of this project focused on four major aspects: (1) The continued development and refinement of multiscale modeling techniques to address the issue of the disparate scales of the physico-chemical processes that govern the fate of atmospheric pollutants; (2) Development and application of analysis methods utilizing process and mass balance techniques to increase the interpretive powers of atmospheric models and to aid in complementary analysis of model predictions and observations; (3) Development of meteorological and emission inputs for initial application of the chemistry/transport model over the north Atlantic region; and, (4) The continued development and implementation of a totally new adaptive chemistry representation that changes the details of what is represented as the underlying conditions change.

Mathur, Rohit↗

Calculation of Flight Deck Interval Management Assigned Spacing Goals Subject to Multiple Scheduling Constraints

The Federal Aviation Administration's Next Generation Air Transportation System will combine advanced air traffic management technologies, performance-based procedures, and state-of-the-art avionics to maintain efficient operations throughout the entire arrival phase of flight. Flight deck Interval Management (FIM) operations are expected to use sophisticated airborne spacing capabilities to meet precise in-trail spacing from top-of-descent to touchdown. Recent human-in-the-loop simulations by the National Aeronautics and Space Administration have found that selection of the assigned spacing goal using the runway schedule can lead to premature interruptions of the FIM operation during periods of high traffic demand. This study compares three methods for calculating the assigned spacing goal for a FIM operation that is also subject to time-based metering constraints. The particular paradigms investigated include: one based upon the desired runway spacing interval, one based upon the desired meter fix spacing interval, and a composite method that combines both intervals. These three paradigms are evaluated for the primary arrival procedures to Phoenix Sky Harbor International Airport using the entire set of Rapid Update Cycle wind forecasts from 2011. For typical meter fix and runway spacing intervals, the runway- and meter fix-based paradigms exhibit moderate FIM interruption rates due to their inability to consider multiple metering constraints. The addition of larger separation buffers decreases the FIM interruption rate but also significantly reduces the achievable runway throughput. The composite paradigm causes no FIM interruptions, and maintains higher runway throughput more often than the other paradigms. A key implication of the results with respect to time-based metering is that FIM operations using a single assigned spacing goal will not allow reduction of the arrival schedule's excess spacing buffer. Alternative solutions for conducting the FIM operation in a manner more compatible with the arrival schedule are discussed in detail.

ATD-1↗

Transport of energetic electrons in a fully ionized hydrogen plasma

A Monte Carlo method for calculating energetic electron transport in a plasma is presented. The energy loss and angular deflection due to Coulomb collisions as a function of travel distance and energy are derived for energetic electrons in a fully ionized plasma with a uniform magnetic field. Formulas which include the effect of nonuniform B fields on the angular deflection are derived. The Monte Carlo method is applied to the thick-target beam model in which the energetic electrons are injected vertically downward.

Bai, T.↗

Interstellar polycyclic aromatic hydrocarbons - The infrared emission bands, the excitation/emission mechanism, and the astrophysical implications

A comprehensive study of the PAH hypothesis is presented, including the interstellar, IR spectral features which have been attributed to emission from highly vibrationally excited PAHs. Spectroscopic and IR emission features are discussed in detail. A method for calculating the IR fluorescence spectrum from a vibrationally excited molecule is described. Analysis of interstellar spectrum suggests that the PAHs which dominate the IR spectra contain between 20 and 40 C atoms. The results are compared with results from a thermal approximation. It is found that, for high levels of vibrational excitation and emission from low-frequency modes, the two methods produce similar results. Also, consideration is given to the relationship between PAH molecules and amorphous C particles, the most likely interstellar PAH molecular structures, the spectroscopic structure produced by PAHs and PAH-related materials in the UV portion of the interstellar extinction curve, and the influence of PAH charge on the UV, visible, and IR regions.

Allamandola, L. J.↗

Method of Obtaining High Resolution Intrinsic Wire Boom Damping Parameters for Multi-Body Dynamics Simulations

The goal of NASA's Magnetospheric MultiScale (MMS) mission is to understand magnetic reconnection with sensor measurements from four spinning satellites flown in a tight tetrahedron formation. Four of the six electric field sensors on each satellite are located at the end of 60- meter wire booms to increase measurement sensitivity in the spin plane and to minimize motion coupling from perturbations on the main body. A propulsion burn however, might induce boom oscillations that could impact science measurements if oscillations do not damp to values on the order of 0.1 degree in a timely fashion. Large damping time constants could also adversely affect flight dynamics and attitude control performance. In this paper, we will discuss the implementation of a high resolution method for calculating the boom's intrinsic damping, which was used in multi-body dynamics simulations. In summary, experimental data was obtained with a scaled-down boom, which was suspended as a pendulum in vacuum. Optical techniques were designed to accurately measure the natural decay of angular position and subsequently, data processing algorithms resulted in excellent spatial and temporal resolutions. This method was repeated in a parametric study for various lengths, root tensions and vacuum levels. For all data sets, regression models for damping were applied, including: nonlinear viscous, frequency-independent hysteretic, coulomb and some combination of them. Our data analysis and dynamics models have shown that the intrinsic damping for the baseline boom is insufficient, thereby forcing project management to explore mitigation strategies.

Yew, Alvin G.↗

Specifications of FIPD Fission Gas Release Data

All fission gas release data stored in the Fuels Irradiation & Physics Database (FIPD) was originally measured using the Gas Assay, Sample and Recharge (GASR) System in the Hot Fuel Examination Facility (HFEF). It is therefore called GASR data in FIPD. During the measurement of a sample, such as an irradiated EBR-II fuel element/capsule, a pinhole-sized region near the top of the element plenum was melted by a laser. Plenum gas then expanded into a calibrated volume (note: in this document, “sample” and “capsule/element” are used interchangeably consistent with GASR documents in FIPD). The pressure rise in the volume was recorded. Helium backfilling and expansion was then performed to determine the sample (e.g., fuel element plenum) volume using Boyle’s Law and assuming ideal gas behavior at constant temperature. With the plenum volume and the recorded pressure rise, the sample (e.g., fuel element plenum) pressure was derived with assumption of ideal gas law behavior. The plenum volume and pressure as well as the cladding temperature during the measurement were collected (GASR data in FIPD). Other records associated with the fission gas release data include: raw GASR data records including volumes and post-puncture pressures of seal head/sealing head and manifold, calibration data, backfilling gas pressure data, and the data analysis records. A sample(s) of the fission gas released from the plenum was collected by the GASR system into sample bottles. The chemical and isotopic composition of the gas sample was analyzed separately from GASR data, and will be discussed in a separate specification. The plenum volume, pressure, and cladding temperature during the measurement are typically utilized to determine the number of moles of gas in the plenum. This quantity is often compared to the number of moles of gas generated by fission events. However, calculating these values and their associated uncertainty is beyond the scope of this document, as it necessitates additional assumptions. The most important document to understand the FIPD fission gas data is the GASR operational manual (title: Gas Assay, Sample and Recharge System (GASR) operation and maintenance manual, HFEF/N OMM 4381, DOC. NO. W0018-0032-ES-00). This manual provides: (1) description of the GASR and the functions of each component (laser drilling, welding, seal head/sealing head, manifold, vacuum system, sample system, purge and gas tag system, etc.); (2) step-by-step guidance on calibrations, operations, and measurements; and (3) maintenance procedures and other details relating to the structure and operation of the GASR. Note that the original GASR operated until 2020. A new GASR with the same design and measurement methodology was installed in 2021. The specifications of the GASR presented on the HFEF website at this time are consistent with the ones given in the operational manual. The methods to calculate the plenum volume and pressures were not included in the operational manual, but were recorded in the legacy data analysis files. Details of the methods are given in Chapter 3.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Simulating dust grain-radiation coupling on a moving mesh

ABSTRACT We present a model for the interaction between dust and radiation fields in the radiation hydrodynamic code arepo-rt, which solves the moment-based radiative transfer equations on an unstructured moving mesh. Dust is directly treated using live simulation particles, each of which represent a population of grains that are coupled to hydrodynamic motion through a drag force. We introduce methods to calculate radiation pressure on and photon absorption by dust grains. By including a direct treatment of dust, we are able to calculate dust opacities and update radiation fields self-consistently based on the local dust distribution. This hybrid scheme coupling dust particles to an unstructured mesh for radiation is validated using several test problems with known analytic solutions, including dust driven via spherically symmetric flux from a constant luminosity source and photon absorption from radiation incident on a thin layer of dust. Our methods are compatible with the multifrequency scheme in arepo-rt, which treats UV, optical photons as single scattered and IR photons as multi scattered. At IR wavelengths, we model heating of and thermal emission from dust. Dust and gas are not assumed to be in local thermodynamic equilibrium but transfer energy through collisional exchange. We estimate dust temperatures by balancing these dust-radiation and dust-gas energy exchange rates. This framework for coupling dust and radiation can be applied in future radiation hydrodynamic simulations of galaxy formation.

79 ASTRONOMY AND ASTROPHYSICS↗

Open momentum space method for the Hofstadter butterfly and the quantized Lorentz susceptibility

Here we develop a generic k ∙ p open momentum space method for calculating the Hofstadter butterfly of both continuum (moiré) models and tight-binding models, where the quasimomentum is directly substituted by the Landau level (LL) operators. By taking a LL cutoff (and a reciprocal lattice cutoff for continuum models), one obtains the Hofstadter butterfly with in-gap spectral flows. For continuum models such as the moiré model for twisted bilayer graphene, our method gives a sparse Hamiltonian, making it much more efficient than existing methods. The spectral flows in the Hofstadter gaps can be understood as edge states on a momentum space boundary, from which one can determine the two integers (t ν , s ν ) of a gap ν satisfying the Diophantine equation. The spectral flows can also be removed to obtain a clear Hofstadter butterfly. While t ν is known as the Chern number, our theory identifies s ν as a dual Chern number for the momentum space, which corresponds to a quantized Lorentz susceptibility γ xy = eBs ν .

2-dimensional systems↗

Model of Operation-and-Maintenance Costs for Photovoltaic Systems

This article presents a method for calculating costs associated with operation and maintenance (O&M) of photovoltaic (PV) systems. It compiles details regarding the cost and frequency of multiple O&M services to estimate annual O&M costs ($\$$/year) for each year of an analysis period, the net present value ($\$$) of life cycle costs accumulated over the analysis period, and the reserve account amount ($\$$). Here we show that this method is an improvement over the previous averaged or levelized per-unit ($\$$/kW/year) valuations for estimating PV O&M costs, because it allows a detailed selection of services to perform based on system size, market served (e.g., residential, commercial, or utility), type and configuration of system components (e.g., micro-, string or central inverter) and site and environmental conditions (e.g., snow, pollen, bird populations). This model also distinguishes costs that vary from year to year and increase at different rates over time because of heuristic failure distributions (e.g. Weibull or Lognormal distribution) based on actuarial data for many of the services. This cost model was created by the PV O&M Working Group of researchers and industry, sponsored by DOE Solar Energy Technologies Office, and has been published in an on-line version hosted by SunSpec Alliance at apsuite.sunspec.org. A spreadsheet version is included with this paper as supplementary material.

14 SOLAR ENERGY↗

Numerical simulation of airfoil ice accretion

A mathematical model of glaze and rime ice accretion on a two-dimensional airfoil is presented. The model employs standard methods for calculating the flow field and cloud droplet trajectories. This data is used as input to a thermodynamic analysis of the ice accretion process which includes liquid runback. The structure of the model allows for dynamic updating of the droplet collection efficiency and heat and mass transfer processes as the accreted ice shape changes. This results in the improved shape prediction over previously described methods. Simulations made with the model are compared to experimental results obtained for a NACA 0012 airfoil in the NASA/Lewis Icing Research Tunnel. The agreement with experiment is found to be generally satisfactory for some cases, but a need for improvements is indicated by others. Suggestions for improvements are made.

Macarthur, C. D.↗

A Prognostic Launch Vehicle Probability of Failure Assessment Methodology for Conceptual Systems Predicated on Human Causal Factors

Create an improved method to calculate reliability of a conceptual launch vehicle system prior to fabrication by using historic data of actual root causes of failures. While failures have unique "proximate causes", there are typically a finite amount of common "root causes". Heretofore launch vehicle reliability evaluation typically hardware-centric statistical analyses, while most root causes of failures are been shown to be human-centric. A method based on human-centric root causes can be used to quantify reliability assessments and focus proposed actions to mitigate problems. Existing methods have been optimistic in their projections of launch vehicle reliability compared to actuals. Hypothesis: reliability of a conceptual launch vehicle can be more accurately evaluated based on a rational, probabilistic approach using past failure assessment teams' findings predicated on human-centric causes."Human Reliability Analysis Methods Selection Guidance for NASA"Chandler F.T., et al., NASA HQ/OSMA study group, July 2006. Outside HRA experts from academia, other federal labs, and the private sector. 50 system reliability methods considered, fourteen selected for further study, four finally selected as best suited for human spaceflight. Probabilistic Risk Analysis (PRA) + Human Reliability Analysis (HRA) enabled incorporating effects and probabilities of human errors. While four down-selected methods deemed appropriate for failure assessment, it did not appear that these methods could be concisely applied to perform major system-wide assessment of probability of failure of a conceptual design without becoming unwieldy."Engineering a Safer World", Detailed, comprehensive study external to NASA Leveson N. G., MIT, 2011.Systems-Theoretic Accident Model and Processes (STAMP). All-encompassing accident model based on systems theory analyzed accidents after they occurred and created approaches to prevent occurrence in developing systems not focused on failure prevention per se, but rather reducing hazards by influencing human behavior through use of constraints, hierarchical control structures, and process models to improve system safetySystem Theoretic Process Analysis (STPA) addresses predictive part of problem (a "hazard analysis"). Includes all causal factors identified in STAMP: "...design errors, software flaws, component interaction accidents, cognitively complex human decision-making errors, and social organizational and management factors contributing to accidents" can guide design process rather than require it to exist before-hand did not appear capable of concise application for system-wide assessment of probability of failure of a conceptual design without becoming unwieldy.

Williams, Craig H.↗

Capillary Pressure to Relative Permeability: Task 7.3.4 Relative Permeability Analysis Final Report

Surface tension affects all aspects of fluid flow in porous media. Through measurements of surface tension interaction under multiphase conditions, a relative permeability relationship can be determined. Relative permeability is a numerical description of the interplay between two or more fluids and the porous media they flow through. It is a critical parameter for various tools used to characterized subsurface multiphase flow systems, such as numerical simulation for oil and gas development, carbon sequestration, and groundwater contamination remediation. Therefore, it is critical to get a good statistic distribution of relative permeability in the porous media under study. Empirical relationships for determining relative permeability from capillary pressure are already well established but do not provide the needed flexibility required to match laboratory derive relative permeability relationships. By expanding the existing methods for calculating relative permeability from capillary pressure data, it is possible to create both two and three-phase relative permeability relationships. Existing laboratory-measured relative permeability data along with mercury intrusion capillary (MICP) data coupled with interfacial tension and contact angle measurements were used to determine the efficacy of this approach to relative permeability curve creation. The relative permeability relationships determined with this method were fit to existing laboratory data to elucidate common fitting parameters that were then used to create relative permeability relationships from MICP data that does not have an associated laboratory-measured relative permeability relationship.

58 GEOSCIENCES↗

Time-Dependent Parabolic Finite Difference Formulation for Harmonic Sound Propagation in a Two-Dimensional Duct with Flow

An explicit finite difference real time iteration scheme is developed to study harmonic sound propagation in aircraft engine nacelles. To reduce storage requirements for future large 3D problems, the time dependent potential form of the acoustic wave equation is used. To insure that the finite difference scheme is both explicit and stable for a harmonic monochromatic sound field, a parabolic (in time) approximation is introduced to reduce the order of the governing equation. The analysis begins with a harmonic sound source radiating into a quiescent duct. This fully explicit iteration method then calculates stepwise in time to obtain the 'steady state' harmonic solutions of the acoustic field. For stability, applications of conventional impedance boundary conditions requires coupling to explicit hyperbolic difference equations at the boundary. The introduction of the time parameter eliminates the large matrix storage requirements normally associated with frequency domain solutions, and time marching attains the steady-state quickly enough to make the method favorable when compared to frequency domain methods. For validation, this transient-frequency domain method is applied to sound propagation in a 2D hard wall duct with plug flow.

Kreider, Kevin L.↗