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,747 records · Page 97
All sourcesNASA NTRSDOE OSTI
NASA NTRS1974

Aerodynamic heating on 3-D bodies including the effects of entropy-layer swallowing

A relatively simple method was developed previously (authors, 1973) for calculating laminar, transitional, and turbulent heating rates on three-dimensional bodies in hypersonic flows. This method was shown to yield reasonably accurate results for laminar heating on blunted circular and elliptical cones and an earlier version of the space shuttle vehicle. As the boundary layer along the surface grows, more and more of the inviscid-flow mass is entrained into the boundary layer, and the streamlines which passed through the nearly normal portion of the bow shock wave are 'swallowed' by the boundary layer. This phenomenon is often referred to as entropy-layer or streamline swallowing, and it can have a significant effect on the calculated heating rates. An approximate, yet simple, method for including the effects of entropy-layer swallowing in the heating-rate calculations is given.

Dejarnette, F. R.↗
NASA NTRS1998

Application of the ENS3DAE Euler/Navier-Stokes Aeroelastic method

This paper discusses recent applications of the ENS3DAE computational aeroelasticity method. In particular, it describes aeroelastic and unsteady aerodynamics calculations performed on wings with trailing edge control surfaces. These simulations include the investigation of control reversal for a structurally flexible wing with a deflected control surface, and a static and dynamic analysis of a rigid wing with an oscillating control surface. The two sets of calculations were performed independently on different wings using different grid topologies. The control reversal simulation represents an inviscid Euler static aeroelastic analysis of a thin wing with a rectangular planform. The geometry of this wing makes it suitable for computations using more approximate, inviscid aerodynamics methods. Thus, the results of the present Euler computations are compared with numerical data generated by a validated computational aeroelasticity code which uses a simpler aerodynamic formulation. The second illustrated case involves the simulation of a significantly more complex flowfield and the static and dynamic analyses of this geometry were performed using the viscous Navier-Stokes equation option in ENS3DAE. Results of both the steady and unsteady calculations on this wing are compared with existing experimental data.

Schuster, David M.↗
DOE OSTI2015

Superfluidity of heated Fermi systems in the static fluctuation approximation

Superfluidity properties of heated finite Fermi systems are studied in the static fluctuation approximation, which is an original method. This method relies on a single and controlled approximation, which permits taking correctly into account quasiparticle correlations and thereby going beyond the independent-quasiparticle model. A closed self-consistent set of equations for calculating correlation functions at finite temperature is obtained for a finite Fermi system described by the Bardeen–Cooper–Schrieffer Hamiltonian. An equation for the energy gap is found with allowance for fluctuation effects. It is shown that the phase transition to the supefluid state is smeared upon the inclusion of fluctuations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1973

On singlet s-wave electron-hydrogen scattering.

Discussion of various zeroth-order approximations to s-wave scattering of electrons by hydrogen atoms below the first excitation threshold. The formalism previously developed by the author (1967, 1968) is applied to Feshbach operators to derive integro-differential equations, with the optical-potential set equal to zero, for the singlet and triplet cases. Phase shifts of s-wave scattering are computed in the zeroth-order approximation of the Feshbach operator method and in the static-exchange approximation. It is found that the convergence of numerical computations is faster in the former approximation than in the latter.

Madan, R. N.↗
NASA NTRS1982

Development of a reduced basis technique for transient thermal analysis

A technique to reduce the degrees of freedom in static and dynamic problems, the reduced basis method, is described. The method combines the classical Rayleigh-Ritz approximation with contemporary finite element methods to retain modeling versatility as the degrees of freedom are reduced. Applications to a nonlinear dynamic response problem are discussed efforts to apply the method to nonlinear transient thermal response problems are summarized. The selection of basis vectors for reducing the system of equations is addressed.

Shore, C. P.↗
NASA NTRS1994

On the wavelet optimized finite difference method

When one considers the effect in the physical space, Daubechies-based wavelet methods are equivalent to finite difference methods with grid refinement in regions of the domain where small scale structure exists. Adding a wavelet basis function at a given scale and location where one has a correspondingly large wavelet coefficient is, essentially, equivalent to adding a grid point, or two, at the same location and at a grid density which corresponds to the wavelet scale. This paper introduces a wavelet optimized finite difference method which is equivalent to a wavelet method in its multiresolution approach but which does not suffer from difficulties with nonlinear terms and boundary conditions, since all calculations are done in the physical space. With this method one can obtain an arbitrarily good approximation to a conservative difference method for solving nonlinear conservation laws.

Jameson, Leland↗
NASA NTRS1976

Evaluation of two swept-infinite-wing potential/viscous-flow computer programs

Two computer programs capable of predicting the potential and viscous interacting flow around wings of infinite aspect ratio was evaluated. The programs are compared in terms of their capabilities, the approximations and the methods of solution used, and the input requirements. Six airfoils, each representative of a class of airfoils, are used as test airfoils. The results predicted by the programs are presented for each airfoil at sweep angles of 0, 20, and 40 degrees over a range of angles of attack. The results show that at zero sweep both programs predicted the aerodynamic coefficients well and generally in good agreement with measurements. At 20 and 40 degrees of sweep, as there are no experimental data available, definitive conclusions cannot be drawn about the accuracy of the predictions although the results are presented and discussed. The execution times are approximately the same for the two programs.

Arimilli, R. V.↗
NASA NTRS1985

A fast, time-accurate unsteady full potential scheme

The unsteady form of the full potential equation is solved in conservation form by an implicit method based on approximate factorization. At each time level, internal Newton iterations are performed to achieve time accuracy and computational efficiency. A local time linearization procedure is introduced to provide a good initial guess for the Newton iteration. A novel flux-biasing technique is applied to generate proper forms of the artificial viscosity to treat hyperbolic regions with shocks and sonic lines present. The wake is properly modeled by accounting not only for jumps in phi, but also for jumps in higher derivatives of phi, obtained by imposing the density to be continuous across the wake. The far field is modeled using the Riemann invariants to simulate nonreflecting boundary conditions. The resulting unsteady method performs well which, even at low reduced frequency levels of 0.1 or less, requires fewer than 100 time steps per cycle at transonic Mach numbers. The code is fully vectorized for the CRAY-XMP and the VPS-32 computers.

Shankar, V.↗
DOE OSTI2014

Quantum confinement in Si and Ge nanostructures: Theory and experiment

The role of quantum confinement (QC) in Si and Ge nanostructures (NSs) including quantum dots, quantum wires, and quantum wells is assessed under a wide variety of fabrication methods in terms of both their structural and optical properties. Structural properties include interface states, defect states in a matrix material, and stress, all of which alter the electronic states and hence the measured optical properties. We demonstrate how variations in the fabrication method lead to differences in the NS properties, where the most relevant parameters for each type of fabrication method are highlighted. Si embedded in, or layered between, SiO{sub 2}, and the role of the sub-oxide interface states embodies much of the discussion. Other matrix materials include Si{sub 3}N{sub 4} and Al{sub 2}O{sub 3}. Si NSs exhibit a complicated optical spectrum, because the coupling between the interface states and the confined carriers manifests with varying magnitude depending on the dimension of confinement. Ge NSs do not produce well-defined luminescence due to confined carriers, because of the strong influence from oxygen vacancy defect states. Variations in Si and Ge NS properties are considered in terms of different theoretical models of QC (effective mass approximation, tight binding method, and pseudopotential method). For each theoretical model, we discuss the treatment of the relevant experimental parameters.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Compact video content representation for video coding using low multi-linear tensor rank approximation with dynamic core tensor order

Tensor is a multi-way representation of an entity that is widely used in many applications of mathematics and physics. The wide range of tensor algebra operations attract the researchers of graphical and vision field to identify unique methods for various applications such as segmentation, compression, indexing and so on. In this paper, low multi-linear rank approximation of a tensor is used for decomposition so as to perform compression on a video that is represented as a 4D tensor. Thus, the encoder has blocks for (1) multilinear rank approximation method so as to aid decomposition, (2) sparsity removal of residual data to reduce the required memory storage to a small level, (3) quantization on the sparse less residual data and (4) LZ77 dictionary coding for the video. The decoder uses a tensor reconstruction block along with residual error correction block that reduces loss/error in the reconstructed video. The size of core tensor is an important factor in rank approximation. It is identified at run time depending on the video contents. This size is determined adaptively using Tikhonov’s regularization method. The best value of core tensor size is identified with the corner of L-curve which preserves the video contents from heavy loss/error. Experimental results of the proposed work on comparison with H.264/AVC compression method indicate that the proposed method achieves up to 10% improvement in compression ratio with the least amount of loss and a maximum signal to noise ratio of 0.14. The proposed system gives better compression ratio than H.265 for videos with more number of frames without compromising quality.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1993

Automated Finite Element Modeling of Wing Structures for Shape Optimization

The displacement formulation of the finite element method is the most general and most widely used technique for structural analysis of airplane configurations. Modem structural synthesis techniques based on the finite element method have reached a certain maturity in recent years, and large airplane structures can now be optimized with respect to sizing type design variables for many load cases subject to a rich variety of constraints including stress, buckling, frequency, stiffness and aeroelastic constraints (Refs. 1-3). These structural synthesis capabilities use gradient based nonlinear programming techniques to search for improved designs. For these techniques to be practical a major improvement was required in computational cost of finite element analyses (needed repeatedly in the optimization process). Thus, associated with the progress in structural optimization, a new perspective of structural analysis has emerged, namely, structural analysis specialized for design optimization application, or.what is known as "design oriented structural analysis" (Ref. 4). This discipline includes approximation concepts and methods for obtaining behavior sensitivity information (Ref. 1), all needed to make the optimization of large structural systems (modeled by thousands of degrees of freedom and thousands of design variables) practical and cost effective.

Harvey, Michael Stephen↗
DOE OSTI2013

Rapid Radiochemical Analyses in Support of Fukushima Nuclear Accident - 13196

There is an increasing need to develop faster analytical methods for emergency response, including emergency soil and air filter samples [1, 2]. The Savannah River National Laboratory (SRNL) performed analyses on samples received from Japan in April, 2011 as part of a U.S. Department of Energy effort to provide assistance to the government of Japan, following the nuclear event at Fukushima Daiichi, resulting from the earthquake and tsunami on March 11, 2011. Of particular concern was whether it was safe to plant rice in certain areas (prefectures) near Fukushima. The primary objectives of the sample collection, sample analysis, and data assessment teams were to evaluate personnel exposure hazards, identify the nuclear power plant radiological source term and plume deposition, and assist the government of Japan in assessing any environmental and agricultural impacts associated with the nuclear event. SRNL analyzed approximately 250 samples and reported approximately 500 analytical method determinations. Samples included soil from farmland surrounding the Fukushima reactors and air monitoring samples of national interest, including those collected at the U.S. Embassy and American military bases. Samples were analyzed for a wide range of radionuclides, including strontium-89, strontium-90, gamma-emitting radionuclides, and plutonium, uranium, americium and curium isotopes. Technical aspects of the rapid soil and air filter analyses will be described. The extent of radiostrontium contamination was a significant concern. For {sup 89,90}Sr analyses on soil samples, a rapid fusion technique using 1.5 gram soil aliquots to enable a Minimum Detectable Activity (MDA) of <1 pCi {sup 89,90}Sr /g of soil was employed. This sequential technique has been published recently by this laboratory for actinides and radiostrontium in soil and vegetation [3, 4]. It consists of a rapid sodium hydroxide fusion, pre-concentration steps using iron hydroxide and calcium fluoride precipitations, followed by Sr-Resin separation and gas flow proportional counting. To achieve a lower detection limit for analysis of some of the Japanese soil samples, a 10 gram aliquot of soil was taken, acid-leached and processed with similar preconcentration chemistry. The MDA using this approach was ∼0.03 pCi/g (1.1 mBq/g)/, which is less than the 0.05-0.10 pCi/g {sup 90}Sr levels found in soil as a result of global fallout. The chemical yields observed for the Japanese soil samples was typically 75-80% and the laboratory control sample (LCS) and matrix spike (MS) results looked very good for this work Individual QC results were well within the ± 25% acceptable range and the average of these results does not show significant bias. Additional data for a radiostrontium in soil method for 50 gram samples will also be presented, which appears to be a significant step forward based on looking at the current literature, with higher chemical yields for even larger sample aliquots and lower MDA [5, 6, 7] Hou et al surveyed a wide range of separation methods for Pu in waters and environmental solid samples [8]. While there are many actinide methods in the scientific literature, few would be considered rapid due to the tedious and time-consuming steps involved. For actinide analyses in soil, a new rapid method for the determination of actinide isotopes in soil samples using both alpha spectrometry and inductively-coupled plasma mass spectrometry was employed. The new rapid soil method utilizes an acid leaching method, iron/titanium hydroxide precipitation, a lanthanum fluoride soil matrix removal step, and a rapid column separation process with TEVA Resin. The large soil matrix is removed easily and rapidly using these two simple precipitations with high chemical recoveries and effective removal of interferences. [9, 10] Vacuum box technology and rapid flow rates were used to reduce analytical time. Challenges associated with the mineral content in the volcanic soil will be discussed. Air filter samples were reported within twenty-four (24) hours of receipt using rapid techniques published previously. [11] The rapid reporting of high quality analytical data arranged through the U.S. Department of Energy Consequence Management Home Team was critical to allow the government of Japan to readily evaluate radiological impacts from the nuclear reactor incident to both personnel and the environment. SRNL employed unique rapid methods capability for radionuclides to support Japan that can also be applied to environmental, bioassay and waste management samples. New rapid radiochemical techniques for radionuclides in soil and other environmental matrices as well as some of the unique challenges associated with this work will be presented that can be used for application to environmental monitoring, environmental remediation, decommissioning and decontamination activities. (authors)

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

Multiscale formulation of frictional contact mechanics at the pore scale

Direct numerical simulation (DNS) yields the highest fidelity predictions of mechanical deformation at the pore scale, but is prohibitively expensive for analyzing large or many samples. Discrete element methods (DEM) are an efficient alternative, but are limited to granular media and incapable of estimating or controlling prediction errors. We present a pore-level multiscale method (PLMM) that approximates DNS efficiently and with controllable accuracy. We focus on the linear elastic response of a consolidated geologic porous medium with arbitrary microstructure, heterogeneous mineralogy, containing cracks or defects. PLMM decomposes the solid phase into non-overlapping subdomains, on which local basis functions are constructed. The bases are then coupled with a global interface problem that accounts for slip or stick contact conditions between the subdomains. PLMM produces an initial, but accurate, approximation to DNS that can be iteratively improved. It is amenable to parallelism and allows for different mesh, models, and physics in each subdomain. An algebraic interpretation of PLMM as a preconditioner is also presented to allow non-intrusive implementation into existing solvers. Lastly, this work extends previous developments of PLMM in fluid dynamics to solid mechanics and enables future extensions towards modeling coupled flow and mechanics problems.

42 ENGINEERING↗
NASA NTRS1986

Solving Nonlinear Coupled Differential Equations

Harmonic balance method developed to obtain approximate steady-state solutions for nonlinear coupled ordinary differential equations. Method usable with transfer matrices commonly used to analyze shaft systems. Solution to nonlinear equation, with periodic forcing function represented as sum of series similar to Fourier series but with form of terms suggested by equation itself.

Mitchell, L.↗
NASA NTRS1989

Nonlinear filters for efficient shock computation

A new type of methods for the numerical approximation of hyperbolic conservation laws with discontinuous solution is introduced. The methods are based on standard finite difference schemes. The difference solution is processed with a nonlinear conservation form filter at every time level to eliminate spurious oscillations near shocks. It is proved that the filter can control the total variation of the solution and also produce sharp discrete shocks. The method is simpler and faster than many other high resolution schemes for shock calculations. Numerical examples in one and two space dimensions are presented.

Engquist, Bjorn↗
DOE OSTI2013

Subcell resolution in simplex stochastic collocation for spatial discontinuities

Subcell resolution has been used in the Finite Volume Method (FVM) to obtain accurate approximations of discontinuities in the physical space. Stochastic methods are usually based on local adaptivity for resolving discontinuities in the stochastic dimensions. However, the adaptive refinement in the probability space is ineffective in the non-intrusive uncertainty quantification framework, if the stochastic discontinuity is caused by a discontinuity in the physical space with a random location. The dependence of the discontinuity location in the probability space on the spatial coordinates then results in a staircase approximation of the statistics, which leads to first-order error convergence and an underprediction of the maximum standard deviation. To avoid these problems, we introduce subcell resolution into the Simplex Stochastic Collocation (SSC) method for obtaining a truly discontinuous representation of random spatial discontinuities in the interior of the cells discretizing the probability space. The presented SSC–SR method is based on resolving the discontinuity location in the probability space explicitly as function of the spatial coordinates and extending the stochastic response surface approximations up to the predicted discontinuity location. The applications to a linear advection problem, the inviscid Burgers’ equation, a shock tube problem, and the transonic flow over the RAE 2822 airfoil show that SSC–SR resolves random spatial discontinuities with multiple stochastic and spatial dimensions accurately using a minimal number of samples.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗