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 397 records · Page 22

Predicting the mineral composition of dust aerosols – Part 2: Model evaluation and identification of key processes with observations

A global compilation of nearly sixty measurement studies is used to evaluate two methods of simulating the mineral composition of dust aerosols in an Earth system model. Both methods are based upon a Mean Mineralogical Table (MMT) that relates the soil mineral fractions to a global atlas of arid soil type. The Soil Mineral Fraction (SMF) method assumes that the aerosol mineral fractions match the fractions of the soil. The MMT is based upon soil measurements after wet sieving, a process that destroys aggregates of soil particles that would have been emitted from the original, undisturbed soil. The second method approximately reconstructs the emitted aggregates. This model is referred to as the Aerosol Mineral Fraction (AMF) method because the mineral fractions of the aerosols differ from those of the wet-sieved parent soil, partly due to reaggregation. The AMF method remedies some of the deficiencies of the SMF method in comparison to observations. Only the AMF method exhibits phyllosilicate mass at silt sizes, where they are abundant according to observations. In addition, the AMF quartz fraction of silt particles is in better agreement with measured values, in contrast to the overestimated SMF fraction. Measurements at distinct clay and silt particle sizes are shown to be more useful for evaluation of the models, in contrast to the sum over all particles sizes that is susceptible to compensating errors, as illustrated by the SMF experiment. Model errors suggest that allocation of the emitted silt fraction of each mineral into the corresponding transported size categories is an important remaining source of uncertainty. Evaluation of both models and the MMT is hindered by the limited number of size-resolved measurements of mineral content that sparsely sample aerosols from the major dust sources. In conclusion, the importance of climate processes dependent upon aerosol mineral composition shows the need for global and routine mineral measurements.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Aerodynamic sensitivities from subsonic, sonic and supersonic unsteady, nonplanar lifting-surface theory

The technique of implicit differentiation has been used in combination with linearized lifting-surface theory to derive analytical expressions for aerodynamic sensitivities (i.e., rates of change of lifting pressures with respect to general changes in aircraft geometry, including planform variations) for steady or oscillating planar or nonplanar lifting surfaces in subsonic, sonic, or supersonic flow. The geometric perturbation is defined in terms of a single variable, and the user need only provide simple expressions or similar means for defining the continuous or discontinuous global or local perturbation of interest. Example expressions are given for perturbations of the sweep, taper, and aspect ratio of a wing with trapezoidal semispan planform. In addition to direct computational use, the analytical method presented here should provide benchmark criteria for assessing the accuracy of aerodynamic sensitivities obtained by approximate methods such as finite geometry perturbation and differencing. The present process appears to be readily adaptable to more general surface-panel methods.

Yates, E. Carson, Jr.↗

The conformal transformation of an airfoil into a straight line and its application to the inverse problem of airfoil theory

A method of conformal transformation is developed that maps an airfoil into a straight line, the line being chosen as the extended chord line of the airfoil. The mapping is accomplished by operating directly with the airfoil ordinates. The absence of any preliminary transformation is found to shorten the work substantially over that of previous methods. Use is made of the superposition of solutions to obtain a rigorous counterpart of the approximate methods of thin-airfoils theory. The method is applied to the solution of the direct and inverse problems for arbitrary airfoils and pressure distributions. Numerical examples are given. Applications to more general types of regions, in particular to biplanes and to cascades of airfoils, are indicated. (author)

Mutterperl, William↗

Application of radial-equilibrium condition to axial-flow compressor and turbine design

Basic general equations governing the three-dimensional compressible flow of gas through a compressor or turbine are given in terms of total enthalpy, entropy, and velocity components of the gas. Two methods of solution are obtained for the simplified, steady axially symmetric flow; one involves the use of a number of successive planes normal to the axis of the machine and short distances apart, and the other involves only three stations for a stage in which an appropriate radial-flow path is used. Methods of calculation for the limiting cases of zero and infinite blade aspect ratios and an approximate method of calculation for finite blade aspect ratio are also given. In these methods, the blade loading and the shape of the annular passage wall may be arbitrarily specified.

Wu, Chung-Hua↗

Predicting the Mineral Composition of Dust Aerosols: Model Evaluation and Identification of Key Processes with Observations - Part 2

A global compilation of nearly sixty measurement studies is used to evaluate two methods of simulating the mineral composition of dust aerosols in an Earth system model. Both methods are based upon a Mean Mineralogical Table (MMT) that relates the soil mineral fractions to a global atlas of arid soil type. The Soil Mineral Fraction (SMF) method assumes that the aerosol mineral fractions match the fractions of the soil. The MMT is based upon soil measurements after wet sieving, a process that destroys aggregates of soil particles that would have been emitted from the original, undisturbed soil. The second method approximately reconstructs the emitted aggregates. This model is referred to as the Aerosol Mineral Fraction (AMF) method because the mineral fractions of the aerosols differ from those of the wet-sieved parent soil, partly due to reaggregation. The AMF method remedies some of the deficiencies of the SMF method in comparison to observations. Only the AMF method exhibits phyllosilicate mass at silt sizes, where they are abundant according to observations. In addition, the AMF quartz fraction of silt particles is in better agreement with measured values, in contrast to the overestimated SMF fraction. Measurements at distinct clay and silt particle sizes are shown to be more useful for evaluation of the models, in contrast to the sum over all particles sizes that is susceptible to compensating errors, as illustrated by the SMF experiment. Model errors suggest that allocation of the emitted silt fraction of each mineral into the corresponding transported size categories is an important remaining source of uncertainty. Evaluation of both models and the MMT is hindered by the limited number of size-resolved measurements of mineral content that sparsely sample aerosols from the major dust sources. The importance of climate processes dependent upon aerosol mineral composition shows the need for global and routine mineral measurements.

aerosols↗

A method for the calculation of external lift, moment, and pressure drag of slender open-nose bodies of revolution at supersonic speeds

An approximate method is presented for the calculation of the external lift, moment, and pressure drag of slender open-nose bodies of revolution of supersonic speeds. The lift, moment, and pressure drag of a typical ram-jet body shape are calculated at Mach numbers 1.45, 1.60, 1.75, and 3.00; and the lift and moment results are compared with available experimental data. The agreement of the calculated lift and moment data with the experimental data is excellent. The pressure-drag comparison was not presented because of the uncertainty of the amount of skin-friction drag present in the experimental results.

Brown, Clinton E↗

Review of radiative transfer methods in scattering atmospheres

The radiative transfer in a scattering plane-parallel atmosphere is discussed, considering the exact analytical, the computational and the approximate methods. Some results of numerical comparisons are given. Finally, the difficulties of realistic atmospheric models are emphasized.

Lenoble, J.↗

Method of successive approximations for the solution of certain problems in aerodynamics

A method of successive approximations for the solution of problems in the fields of diffusion, boundary-layer flow, and heat-transfer is illustrated by solving problems in each of these fields. In most of the examples, the approximate solutions are compared with known accurate solutions and the agreement is shown to be good.

FLOW, TURBULENT - DIFFUSION↗

Solving high-dimensional partial integral differential equations: The finite expression method

Partial integro-differential equations (PIDEs) have broad applications in the sciences, from electro-magnetism to options pricing. Here, in this paper, we introduce a new finite expression method (FEX) to solve PIDEs. This approach builds upon the original FEX and its inherent advantages with new advances: 1) A novel method of parameter grouping is proposed to reduce the number of coefficients in high-dimensional function approximation; 2) A Taylor series approximation method is implemented to significantly improve the computational efficiency and accuracy of the evaluation of the integral terms of PIDEs. The new FEX based method, denoted FEX-PG to indicate the addition of the parameter grouping (PG) step to the algorithm, provides both high accuracy and interpretable numerical solutions, with the outcome being an explicit equation that facilitates intuitive understanding of the underlying solution structures. These features are often absent in traditional methods, such as finite element methods (FEM) and finite difference methods, as well as in deep learning-based approaches. To benchmark our method against recent advances, we apply the new FEX-PG to solve benchmark PIDEs in the literature. In high-dimensional settings, FEX-PG exhibits strong and robust performance, achieving relative errors on the order of single precision machine epsilon, significantly outperforming existing approaches based on neural networks.

Combinatorial optimization↗

Combining global and local approximations

A method based on a linear approximation to a scaling factor, designated the 'global-local approximation' (GLA) method, is presented and shown capable of extending the range of usefulness of derivative-based approximations to a more refined model. The GLA approach refines the conventional scaling factor by means of a linearly varying, rather than constant, scaling factor. The capabilities of the method are demonstrated for a simple beam example with a crude and more refined FEM model.

Haftka, Raphael T.↗

Depletion capabilities in the OpenMC Monte Carlo particle transport code

A depletion solver has been implemented in OpenMC and is described herein. The depletion solver is implemented in Python and interfaces with OpenMC’s transport solver through a C++ application programming interface, which enables an in-memory transport-depletion coupling. Multiple integration methods for advancing in time have been implemented and exhibit tradeoffs in cost, accuracy, and memory use. For all time integration methods, evaluation of the matrix exponential is performed by using the incomplete partial fraction form of the Chebyshev rational approximation method. Simulations of a pressurized water reactor (PWR) pincell and a sodium-cooled fast reactor (SFR) assembly were carried out with OpenMC and Serpent. For both problems, the use of a high-fidelity depletion chain results in predictions of k eff that agree within 20–30 pcm between OpenMC and Serpent. Predicted actinide concentrations were found to agree to a fraction of a percent, and most fission product concentrations were found to agree within 1%. Here, the few cases where larger differences were observed can be attributed either to differences in how the energy dependence of fission product yields is handled or deficiencies in the nuclear data used. OpenMC simulations of the PWR and SFR problems using a simplified 228-nuclide depletion chain demonstrate that it achieves accuracy close to that of the full, high-fidelity depletion chain with respect to the studied responses.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Zero-point energy, tunneling, and vibrational adiabaticity in the Mu + H2 reaction

Abstract: Isotopic substitution of muonium for hydrogen provides an unparalleled opportunity to deepen our understanding of quantum mass effects on chemical reactions. A recent topical review [Aldegunde et al., Mol. Phys. 111, 3169 (2013)] of the thermal and vibrationally-stateselected reaction of Mu with H2 raises a number of issues that are addressed here. We show that some earlier quantum mechanical calculations of the Mu + H2 reaction, which are highlighted in this review and which have been used to benchmark approximate methods, are in error by as much as 19% in the low-temperature limit. We demonstrate that an approximate treatment of the Born–Oppenheimer diagonal correction that was used in some recent studies is not valid for treating the vibrationally-state-selected reaction. We also discuss why vibrationally adiabatic potentials that neglect bend zero-point energy are not a useful analytical tool for understanding reaction rates and why vibrationally nonadiabatic transitions cannot be understood by considering tunneling through vibrationally adiabatic potentials. Finally, we present calculations on a hierarchical family of potential energy surfaces to assess the sensitivity of rate constants to the quality of the potential surface.

Born–Oppenheimer diagonal correction, kinetic isot↗

Application of axisymmetric analogue for calculating heating in three-dimensional flows

A rapid, approximate method has been developed for calculating the heating rates on three-dimensional vehicles such as the Space Shuttle Orbiter and other advanced reentry configurations. The method is based on the axisymmetric analogue for three-dimensional boundary layers. It uses information obtained from a three-dimensional inviscid flowfield solution, such as HALIS, to calculate inviscid surface streamlines along which approximate heating rates are calculated independent of what happens along other streamlines. Three-dimensional effects are included through the metric coefficient that describes the divergence or convergence of streamlines. Boundary-layer edge properties are obtained from the inviscid flowfield solution by interpolating in the inviscid flowfield at a distance equal to the boundary-layer thickness away from the wall. This accounts, approximately, for the variable boundary-layer edge entropy. Using this method, heating calculations can be made along a typical streamline in a few seconds. This method has been used to accurately predict heating rates for simple shapes such as a spherically blunted cone and more complex shapes such as the Shuttle Orbiter for a variety of wind-tunnel and flight conditions. A unique feature of the method is its ability to accurately predict heating rates on the Shuttle Orbiter wing.

Hamilton, H. H.↗

The reduction of airplane flight test data to standard atmosphere conditions

This report was prepared for the National Advisory Committee for Aeronautics in order to supply the need of practical methods of reducing observed performance to standard conditions with a minimum of labor. The first part gives a very simple approximate method of reducing performance in climb, and is particularly adapted to work not requiring extreme accuracy. The second part gives a somewhat more elaborate and more accurate method which is well suited to general flight test reduction. The third part gives the conventional method of calibrating air-speed indicators and reducing the indicated speeds to true air speeds. An appendix gives working tables and charts for the standard atmosphere. (author)

Diehl, Walter S↗

Polarimetric thermal emission from rough surfaces

Recent theoretical works have suggested the potential of passive polarimetry in the remote sensing of geophysical media. It was shown that the third Stokes parameter U of the thermal emission may become larger for azimuthally asymmetric fields of observation. In order to investigate the potential applicability of passive polarimetry to the remote sensing of ocean surface, measurements of the polarimetric thermal emission from a sinusoidal water surface and a numerical study of the polarimetric thermal emission from randomly rough ocean surfaces were performed. Measurements of sinusoidal water surface thermal emission were performed using a sinusoidal water surface which was created by placing a thin sheet of fiberglass with a sinusoidal profile in two dimensions extended infinitely in the third dimension onto a water surface. The theory of thermal emission from a 'two-layer' periodic surface is derived and the exact solution is performed using both the extended boundary condition method (EBC) and the method of moments (MOM). The theoretical predictions are found to be in good agreement with the experimental results once the effects of the radiometer antenna pattern are included and the contribution of background noise to the measurements is modeled. The experimental results show that the U parameter indicates the direction of periodicity of the water surface and can approach values of up to 30 K for the surface observed. Next, a numerical study of polarimetric thermal emission from randomly rough surfaces was performed. A Monte Carlo technique utilizing an exact method for calculating thermal emission was chosen for the study to avoid any of the limitations of the commonly used approximate methods in rough surface scattering. In this Monte Carlo technique, a set of finite rough surface profiles in two dimensions with desired statistics was generated and extended periodically. The polarimetric thermal emission from each surface of the set was then calculated using both the EBC and the MOM and the results were averaged. The surface statistics chosen were intended to model a wind perturbed ocean surface in the X to Ku band microwave region. The results indicate that the U parameter is sensitive to the azimuthal angle between the surface periodicity and the looking angle and to the rms height of the surface, and that the U parameter is fairly insensitive to variations in polar angle, permittivity, surface power law spectrum, and surface spectrum high frequency cutoff. These properties give further strength to the idea of using the U parameter to detect wind direction over the ocean.

Johnson, J. T.↗

Numerical Results for Polarized Light Scattering in a Spherical Atmosphere

We report numerical results for polarized light reflection from the top of a Rayleigh scattering spherical atmosphere with height-dependent single scattering albedo over a dark surface. Michael Mishchenko considered this scenario back in the 1990’s, for a plane-parallel atmosphere of unit optical thickness (OT), for which radiance errors arising from neglecting polarization reaches their highest values. To further extend Mishchenko’s results, we consider a value of OT = 0.25, for which the effect of atmospheric curvature is pronounced. New results are generated using three state-of-the art radiative transfer (RT) codes. These are: the MYSTIC and MCSSA models, which simulate light scattering in a true- spherical atmosphere using Monte Carlo methods; and the discrete ordinate code VLIDORT, operating with a new multiple-scatter spherical correction designed to deliver reasonable approximations to spherical-medium scattering. In this work, we report results for both single and multiple scattering; this will help to support the validation of existing and future polarized spherical RT codes, especially those using approximative methods to deal with sphericity.

S Korkin↗