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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Late-time small body disruptions for planetary defense

Here, diverting hazardous small bodies on impact trajectories with the Earth can in some circumstances be impossible without risking disrupting them. Disruption is a much more difficult planetary defense scenario to assess, being linked both to the response of the body to shock loading and the much more complicated gravitational dynamics of the fragments in the solar system relative to pure deflection scenarios. In this work we present a new simulation suite built on N-body gravitational methods that solves fragment orbits in the full gravitational system without recourse to more approximate methods. We assess the accuracy of our simulations and the simplifying assumptions we adopt to make the system tractable, and then discuss in more detail several specific, plausible planetary defense scenarios based on real close encounters. We find that disruption can be a very effective planetary defense strategy even for very late (sub-year) interventions, and should be considered an effective backup strategy should preferred methods, which require long warning times, fail.

42 ENGINEERING↗
DOE OSTI2015

Ion energy-angle distribution functions at the plasma-material interface in oblique magnetic fields

The ion energy-angle distribution (IEAD) at the wall of a magnetized plasma is of fundamental importance for the determination of the material processes occurring at the plasma-material interface, comprising secondary emissions and material sputtering. Here, we present a numerical characterization of the IEAD at the wall of a weakly collisional magnetized plasma with the magnetic field inclined at an arbitrary angle with respect to the wall. The analysis has been done using two different techniques: (1) a fluid-Monte Carlo method, and (2) particle-in-cell simulations, the former offering a fast but approximate method for the determination of the IEADs, the latter giving a computationally intensive but self-consistent treatment of the plasma behavior from the quasi-neutral region to the material boundary. The two models predict similar IEADs, whose similarities and differences are discussed. Data are presented for magnetic fields inclined at angles from normal to grazing incidence (0°–85°). We show the scaling factors of the average and peak ion energy and trends of the pitch angle at the wall as a function of the magnetic angle, for use in the correlation of fluid plasma models to material models.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1974

A recurrence matrix method for the analysis of longitudinal and torsional vibrations in non-uniform multibranch beams with variable boundary conditions

An approximate method for calculating the longitudinal and torsional natural frequencies and associated modal data of a beamlike, variable cross section multibranch structure is presented. The procedure described is the numerical integration of the first order differential equations that characterize the beam element in longitudinal motion and that satisfy the appropriate boundary conditions.

Davis, R. B.↗
NASA NTRS1979

Comparison of three-dimensional panel methods with strip boundary-layer simulations to experiment

Several three dimensional surface panel methods coupled on a streamwise strip basis to two dimensional boundary layer methods were compared to experimental data. Several different boundary layer simulations were investigated. Comparable boundary layer effects which converge to the expected two dimensional result wind increasing aspect ratio were indicated by either a transpiration or camber line displacement simulation. The method was applied to several varied configurations. The lift predictions agreed well with an approximate method based on correcting the inviscid lift using the two dimensional section viscous characteristics. The wing pressures and the effect of the boundary layer on the wing pressures were predicted accurately.

Kjelgaard, S. O.↗
NASA NTRS1981

Helicopter rotor loads using a matched asymptotic expansion technique

The theoretical basis and computational feasibility of the Van Holten method, and its performance and range of validity by comparison with experiment and other approximate methods was examined. It is found that within the restrictions of incompressible, potential flow and the assumption of small disturbances, the method does lead to a valid description of the flow. However, the method begins to break down under conditions favoring nonlinear effects such as wake distortion and blade/rotor interaction.

Pierce, G. A.↗
NASA NTRS1993

Nongray gas analyses for reflecting walls utilizing a flux technique

A flux formulation for a planar slab of molecular gas radiation bounded by diffuse reflecting walls is developed. While this formulation is limited to the planar geometry, it is useful for studying approximations necessary in modeling nongray radiative heat transfer. The governing equations are derived by considering the history of multiple reflections between the walls. Accurate solutions are obtained by explicitly accounting for a finite number of reflections and approximating the spectral effects of the remaining reflections. Four approximate methods are presented and compared using a single absorption band of H2O. All four methods reduce to an identical zeroth-order formulation, which accounts for all reflections approximately but does handle nonreflected radiation correctly. A single absorption band of CO2 is also considered using the best-behaved approximation for higher orders. A zeroth-order formulation is sufficient to predict the radiative transfer accurately for many cases considered. For highly reflecting walls, higher order solutions are necessary for better accuracy. Including all the important bands of H2O, the radiative source distributions are also obtained for two different temperature and concentration profiles.

Menart, J. A.↗
DOE OSTI2016

Application of Variance Reduction Techniques in RMC Burnup and Thermal-Hydraulic Coupled Calculation

In reactor analog simulation, the Monte Carlo method has relatively high fidelity since it can reduce calculation errors resulting from using approximation methods by describing complex geometry and neutron energy spectra. The Monte Carlo method is, in essence, a stochastic simulation of the neutron transport process, and statistical uncertainty of all the computed results must be provided. In large-scale reactor computation, this method can accurately estimate integral quantities such as the effective multiplication factor (K{sub eff}), and it can limit the uncertainty within an acceptable range via reasonably selecting simulation particles. Nevertheless, it faces great challenges in estimating local quantities like neutron flux and fission power, with uncertainty of estimators beyond acceptance. To decrease the uncertainty of local quantities by increasing simulation particles may cost much time and reduce efficiency. The main reasons are as follows. Firstly, the scale of reactors is so large that common pressurized water reactors have tens of thousands of lattice cells; secondly, uneven distribution of neutron flux and fission power causes uneven distribution of uncertainty. Smith K. puts forward the 95/95 principle to estimate whether uncertainty calculated with Monte Carlo Code satisfies the requirement. According to the principle, with a confidence level of 95%, statistical errors of power density in more than 95% of the regions should be less than 1%. To meet the requirement of the principle, the uniform fission site (UFS) method is raised currently. UFS aims at redistributing fission sources to uniformize the fission source sampling, which can make statistical variance distribution more uniform. This method, which is easy to use, will not obviously increase computation time. Many studies show that UFS has a good variance reduction effect in Monte Carlo criticality calculation, but less attention is paid to its performance in burnup calculation and thermal-hydraulic coupled calculation. Uncertainty arising from the Monte Carlo method has a large impact on burnup calculation and thermal-hydraulic coupled calculation. Variance propagates between burnup steps and between iterative steps of coupled calculation, which causes asymmetry of results of the physical symmetrical regions and instability of coupled iterative convergence errors. To further verify the effect of UFS, this study uses it in the RMC (the Reactor Monte Carlo code) burnup calculation module and the neutronics/thermohydraulics coupled calculation module, and conducts two-dimensional pressurized water reactor burnup calculation and BEAVRS benchmark thermalhydraulic coupled calculation. In this way, it verifies the improvement effect of UFS on variance propagation. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS1971

Investigation of thermal protection systems effects on viscid and inviscid flow fields for manned entry systems

Procedures and methods for predicting aerothermodynamic heating to delta orbiter shuttle vehicles were reviewed. A number of approximate methods were found to be adequate for large scale parameter studies, but are considered inadequate for final design calculations. It is recommended that final design calculations be based on a computer code which accounts for nonequilibrium chemistry, streamline spreading, entropy swallowing, and turbulence. It is further recommended that this code be developed with the intent that it can be directly coupled with an exact inviscid flow field calculation when the latter becomes available. A nonsimilar, equilibrium chemistry computer code (BLIMP) was used to evaluate the effects of entropy swallowing, turbulence, and various three dimensional approximations. These solutions were compared with available wind tunnel data. It was found study that, for wind tunnel conditions, the effect of entropy swallowing and three dimensionality are small for laminar boundary layers but entropy swallowing causes a significant increase in turbulent heat transfer. However, it is noted that even small effects (say, 10-20%) may be important for the shuttle reusability concept.

Bartlett, E. P.↗
NASA NTRS2008

Accuracy of Gradient Reconstruction on Grids with High Aspect Ratio

Gradient approximation methods commonly used in unstructured-grid finite-volume schemes intended for solutions of high Reynolds number flow equations are studied comprehensively. The accuracy of gradients within cells and within faces is evaluated systematically for both node-centered and cell-centered formulations. Computational and analytical evaluations are made on a series of high-aspect-ratio grids with different primal elements, including quadrilateral, triangular, and mixed element grids, with and without random perturbations to the mesh. Both rectangular and cylindrical geometries are considered; the latter serves to study the effects of geometric curvature. The study shows that the accuracy of gradient reconstruction on high-aspect-ratio grids is determined by a combination of the grid and the solution. The contributors to the error are identified and approaches to reduce errors are given, including the addition of higher-order terms in the direction of larger mesh spacing. A parameter GAMMA characterizing accuracy on curved high-aspect-ratio grids is discussed and an approximate-mapped-least-square method using a commonly-available distance function is presented; the method provides accurate gradient reconstruction on general grids. The study is intended to be a reference guide accompanying the construction of accurate and efficient methods for high Reynolds number applications

Thomas, James↗
DOE OSTI2022

Photoionization of Atomic Systems Using the Random-Phase Approximation Including Relativistic Interactions

Approximation methods are unavoidable in solving a many-electron problem. One of the most successful approximations is the random-phase approximation (RPA). Miron Amusia showed that it can be used successfully to describe atomic photoionization processes of many-electron atomic systems. In this article, the historical reasons behind the term “random-phase approximation” are revisited. A brief introduction to the relativistic RPA (RRPA) developed by Walter Johnson and colleagues is provided and some of its illustrative applications are presented.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1980

Theoretical mean colors for RR Lyrae variables

A hydrodynamically pulsating 0.6 solar mass model of a typical RR Lyrae variable was studied with a radiation transport-hydrodynamic computer program to predict theoretical T sub 3 and colors at many phases and to find the proper methods for getting mean colors and the consequent mean effective temperatures. The variable Eddington radiation approximation method was used with gray and with multifrequency absorption coefficients to represent the radiation flow in the outer optically thin layers. Comparison between observed and computed B-V colors indicate that these low Z population 2 models are reasonably accurate using King 1A composition opacities. The well known Oke, Giver, and Searle relation between B-V and T sub e reproduced. Mean colors were found by four different averaging methods. The method that gives a mean color and the mean T sub e closest to the nonpulsating model was the separate intensity means of B and V.

Davis, C. G.↗
NASA NTRS1994

A rapid radiative transfer model for reflection of solar radiation

A rapid analytical radiative transfer model for reflection of solar radiation in plane-parallel atmospheres is developed based on the Sobolev approach and the delta function transformation technique. A distinct advantage of this model over alternative two-stream solutions is that in addition to yielding the irradiance components, which turn out to be mathematically equivalent to the delta-Eddington approximation, the radiance field can also be expanded in a mathematically consistent fashion. Tests with the model against a more precise multistream discrete ordinate model over a wide range of input parameters demonstrate that the new approximate method typically produces average radiance differences of less than 5%, with worst average differences of approximately 10%-15%. By the same token, the computational speed of the new model is some tens to thousands times faster than that of the more precise model when its stream resolution is set to generate precise calculations.

Xiang, X.↗
NASA NTRS1972

Adjoint variational methods in nonconservative stability problems.

A general nonself-adjoint eigenvalue problem is examined and it is shown that the commonly employed approximate methods, such as the Galerkin procedure, the method of weighted residuals and the least square technique lack variational descriptions. When used in their previously known forms they do not yield stationary eigenvalues and eigenfunctions. With the help of an adjoint system, however, several analogous variational descriptions may be developed and it is shown in the present study that by properly restating the method of least squares, stationary eigenvalues may be obtained. Several properties of the adjoint eigenvalue problem, known only for a restricted group, are shown to exist for the more general class selected for study.

Prasad, S. N.↗
DOE OSTI2020

Connected three-body terms in single-reference unitary many-body theories: Iterative and perturbative approximations

This work introduces various approaches to include connected three-body terms in unitary many-body theories, focusing on the driven similarity renormalization group (DSRG). Starting from the least approximate method—the linearized DSRG truncated to one-, two-, and three-body operators [LDSRG(3)]—we develop several approximate LDSRG(3) models with reduced computational cost. Through a perturbative analysis, we motivate a family of iterative LDSRG(3)-n and -n' (n = 1, 2, 3, 4) methods that contain a subset of the LDSRG(3) diagrams. Among these variants, the LDSRG(3)-2 scheme has the same computational complexity of coupled cluster theory with singles, doubles, and triples (CCSDT), but it outperforms CCSDT in the accuracy of the predicted correlation energies. We also propose and implement two perturbative triples corrections based on the linearized DSRG truncated to one- and two-body operators augmented with recursive semi-quadratic commutators [qDSRG(2)]. Overall, the resulting qDSRG(2)+(T) approach matches the accuracy of the “gold-standard” coupled cluster theory with singles, doubles, and perturbative triples model on the energetics of twenty-eight closed-shell atoms and small molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

A method for including external feed in depletion calculations with CRAM and implementation into ORIGEN

A method for including external feed with polynomial time dependence in depletion calculations with the Chebyshev Rational Approximation Method (CRAM) is presented and the implementation of CRAM to the ORIGEN module of the SCALE suite is described. In addition to being able to handle time-dependent feed rates, the new solver also adds the capability to perform adjoint calculations. Results obtained with the new CRAM solver and the original depletion solver of ORIGEN are compared to high precision reference calculations, which shows the new solver to be orders of magnitude more accurate. Lastly, in most cases, the new solver is up to several times faster due to not requiring similar substepping as the original one.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

The regulation mechanism of yitJ and metF riboswitches

Riboswitches which function at the transcriptional level are sensitive to cotranscriptional folding. Based on the recently proposed theory of cotranscriptional folding, we developed a transition node approximation method to effectively decrease the conformation space of long RNA chains. Our results indicate that this approximation is reliable for calculating the cotranscriptional folding kinetics of long mRNA chains. We theoretically studied the cotranscriptional folding behavior of the yitJ and metF riboswitches in the absence/presence of S-adenosylmethionine. Although the two S-box riboswitches have similar OFF-state structures and share common features of riboswitches operated at the transcriptional level, their regulation mechanisms are different. The yitJ riboswitch is regulated by a combination of thermodynamic and kinetic mechanisms, while the metF riboswitch is solely kinetically controlled. For the yitJ riboswitch, transcriptional pausing at the U-stretch directly following the terminator decreases the amount of ligand required to trigger the switch. The different regulation mechanisms and binding affinities of the two riboswitches result from the different lengths of the anti-terminator helix, which in yitJ is short and only disrupts helix P1 of the riboswitch aptamer, but in metF is long and breaks both the helices P1 and P4.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

Majorant recursions to determine eigenstate bounds of a symmetric exponential quantum anharmonic oscillator

The determination of the energy states is highly studied issue in the quantum mechanics. Based on expectation values dynamics, energy states can be observed. But conditions and calculations vary depending on the created system. In this work, a symmetric exponential anharmonic oscillator is considered and development of a recursive approximation method is studied to find its ground energy state. The use of majorant values facilitates the approximate calculation of expectation values.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗