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,693 records · Page 94
All sourcesNASA NTRSDOE OSTI
NASA NTRS1986

Uniformly high order accurate essentially non-oscillatory schemes 3

In this paper (a third in a series) the construction and the analysis of essentially non-oscillatory shock capturing methods for the approximation of hyperbolic conservation laws are presented. Also presented is a hierarchy of high order accurate schemes which generalizes Godunov's scheme and its second order accurate MUSCL extension to arbitrary order of accuracy. The design involves an essentially non-oscillatory piecewise polynomial reconstruction of the solution from its cell averages, time evolution through an approximate solution of the resulting initial value problem, and averaging of this approximate solution over each cell. The reconstruction algorithm is derived from a new interpolation technique that when applied to piecewise smooth data gives high-order accuracy whenever the function is smooth but avoids a Gibbs phenomenon at discontinuities. Unlike standard finite difference methods this procedure uses an adaptive stencil of grid points and consequently the resulting schemes are highly nonlinear.

Harten, A.↗
NASA NTRS1995

Feasibility of Coupling Between a Single-Mode Elliptical-Core Fiber and a Single Mode Rib Waveguide Over Temperature

To determine the feasibility of coupling the output of an optical fiber to a rib waveguide in a temperature environment ranging from 20 C to 300 C, a theoretical calculation of the coupling efficiency between the two was investigated. This is a significant problem which needs to be addressed to determine whether an integrated optic device can function in a harsh temperature environment. Because the behavior of the integrated-optic device is polarization sensitive, a polarization-preserving optic fiber, via its elliptical core, was used to couple light with a known polarization into the device. To couple light energy efficiently from an optical fiber into a channel waveguide, the design of both components should provide for well-matched electric field profiles. The rib waveguide analyzed was the light input channel of an integrated-optic pressure sensor. Due to the complex geometry of the rib waveguide, there is no analytical solution to the wave equation for the guided modes. Approximation or numerical techniques must be utilized to determine the propagation constants and field patterns of the guide. In this study, three solution methods were used to determine the field profiles of both the fiber and guide: the effective-index method (EIM), Marcatili's approximation, and a Fourier method. These methods were utilized independently to calculate the electric field profile of a rib channel waveguide and elliptical fiber at two temperatures, 20 C and 300 C. These temperatures were chosen to represent a nominal and a high temperature that the device would experience. Using the electric field profile calculated from each method, the theoretical coupling efficiency between the single-mode optical fiber and rib waveguide was calculated using the overlap integral and results of the techniques compared. Initially, perfect alignment was assumed and the coupling efficiency calculated. Then, the coupling efficiency calculation was repeated for a range of transverse offsets at both temperatures. Results of the calculation indicate a high coupling efficiency can be achieved when the two components were properly aligned. The coupling efficiency was more sensitive to alignment offsets in the y direction than the x, due to the elliptical modal profile of both components. Changes in the coupling efficiency over temperature were found to be minimal.

Tuma, Margaret L.↗
DOE OSTI2015

On the theory of the galvanomagnetic properties of composite materials: Lattice model

The problem of the galvanomagnetic properties of composite materials is formulated for a lattice model. The effective galvanomagnetic characteristics of a weakly heterogeneous lattice are determined in the quadratic approximation in the deviation of local conductivity tensor σ̂ (r) from average value <σ̂>. In the case of a low concentration (c ≪ 1) of “defect” bonds, effective conductivity tensor σ̂{sub e} of a binary lattice model is calculated in the c-linear approximation. Effective medium method equations are derived for the formulated lattice problem, and the results are compared with the results obtained in a continuous medium model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Conduction of the three-dimensional model of a composite with structural anisotropy

The electric conduction of the model of a structurally anisotropic composite that represents an isotropic matrix with disk-shaped inclusions (oblate spheroids) is analytically studied. The disks have the same orientation and their centers are randomly distributed in the composite volume. The cases of dielectric and ideally conducting inclusions are considered. The conduction of the model in a broad concentration range is described by the approximate effective medium method. Conduction in the critical region (in the vicinity of the metal–dielectric phase transition point) is studied using the similarity hypothesis.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1978

On least squares approximations to indefinite problems of the mixed type

A least squares method is presented for computing approximate solutions of indefinite partial differential equations of the mixed type such as those that arise in connection with transonic flutter analysis. The method retains the advantages of finite difference schemes namely simplicity and sparsity of the resulting matrix system. However, it offers some great advantages over finite difference schemes. First, the method is insensitive to the value of the forcing frequency, i.e., the resulting matrix system is always symmetric and positive definite. As a result, iterative methods may be successfully employed to solve the matrix system, thus taking full advantage of the sparsity. Furthermore, the method is insensitive to the type of the partial differential equation, i.e., the computational algorithm is the same in elliptic and hyperbolic regions. In this work the method is formulated and numerical results for model problems are presented. Some theoretical aspects of least squares approximations are also discussed.

Fix, G. J.↗
NASA NTRS1994

Development of advanced Navier-Stokes solver

The objective of research was to develop and validate new computational algorithms for solving the steady and unsteady Euler and Navier-Stokes equations. The end-products are new three-dimensional Euler and Navier-Stokes codes that are faster, more reliable, more accurate, and easier to use. The three-dimensional Euler and full/thin-layer Reynolds-averaged Navier-Stokes equations for compressible/incompressible flows are solved on structured hexahedral grids. The Baldwin-Lomax algebraic turbulence model is used for closure. The space discretization is based on a cell-centered finite-volume method augmented by a variety of numerical dissipation models with optional total variation diminishing limiters. The governing equations are integrated in time by an implicit method based on lower-upper factorization and symmetric Gauss-Seidel relaxation. The algorithm is vectorized on diagonal planes of sweep using two-dimensional indices in three dimensions. Convergence rates and the robustness of the codes are enhanced by the use of an implicit full approximation storage multigrid method.

Yoon, Seokkwan↗
NASA NTRS1991

Procedure for computer-controlled milling of accurate surfaces of revolution for millimeter and far-infrared mirrors

A simple method for milling accurate off-axis parabolic mirrors with a computer-controlled milling machine is discussed. For machines with a built-in circle-cutting routine, an exact paraboloid can be milled with few computer commands and without the use of the spherical or linear approximations. The proposed method can be adapted easily to cut off-axis sections of elliptical or spherical mirrors.

Emmons, Louisa↗
DOE OSTI2022

Virtual diffraction simulations using the quasi-coarse-grained dynamics method to understand and interpret plasticity contributions during in situ shock experiments

The dynamic deformation response of metallic materials has contributions from dislocations, deformation twinning, and plastic deformation. The current state-of-art techniques can detail the complex mechanistic history of deformation modes under shock loading in real-time using in situ X-ray diffraction (XRD). However, the capability of these experiments to unravel plasticity contributions is challenging due to limitations in interpreting results and the lack of validation from atomistic simulations. Molecular dynamics (MD) simulations can successfully capture various deformation modes in metals and complement experiments using simulated diffractograms at various stages of evolution. However, the difference in length and time scales of MD simulations and experiments is a substantial obstacle in completing and interpreting in situ diffractograms. Therefore, the existing modeling methods require various approximations to model defect evolution and interaction at the mesoscales. However, while using approximations to correlate the peak broadening behavior to the density of dislocations or the shifts/splitting due to the presence of twins, the interpretations of the plasticity contributions from diffractograms are non-trivial, especially when multiple modes of deformation may be operating. This viewpoint discusses combining a mesoscale modeling method called quasi-coarse-grained dynamics and virtual XRD to characterize the plasticity contributions in BCC metals from slip, twinning, and phase transformation behavior. Furthermore, the combined approach shows promise in bridging the mesoscale gap between the capabilities of atomic-scale simulations and in situ experiments to characterize the dynamic deformation of materials.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Identification of stable adsorption sites and diffusion paths on nanocluster surfaces: an automated scanning algorithm

The diverse coordination environments on the surfaces of discrete, three-dimensional (3D) nanoclusters contribute significantly to their unique catalytic properties. Identifying the numerous adsorption sites and diffusion paths on these clusters is however tedious and time-consuming, especially for large, asymmetric nanoclusters. Here, we present a simple, automated method for constructing approximate 2D potential energy surfaces for the adsorption of atomic species on the surfaces of 3D nanoclusters with minimal human intervention. These potential energy surfaces fully characterize the important adsorption sites and diffusion paths on the nanocluster surfaces with accuracies similar to current approaches and at comparable computational cost. Our method can treat complex nanoclusters, such as alloy nanoclusters, and accounts for cluster relaxation and adsorbate-induced reconstruction, important for obtaining accurate energetics. Moreover, its highly parallelizable nature is ideal for modern supercomputer architectures. We showcase our method using two clusters: Au 18 and Pt 55 . For Au 18 , diffusion of atomic hydrogen between the most stable sites occurs via non-intuitive paths, underlining the necessity of exploring the complete potential energy surface. By enabling the rapid and unbiased assessment of adsorption and diffusion on large, complex nanoclusters, which are particularly difficult to handle manually, our method will help advance materials discovery and the rational design of catalysts.

36 MATERIALS SCIENCE↗
NASA NTRS1987

Error analysis for spectral approximation of the Korteweg-De Vries equation

The conservation and convergence properties of spectral Fourier methods for the numerical approximation of the Korteweg-de Vries equation are analyzed. It is proved that the (aliased) collocation pseudospectral method enjoys the same convergence properties as the spectral Galerkin method, which is less effective from the computational point of view. This result provides a precise mathematical answer to a question raised by several authors in recent years.

Maday, Y.↗
NASA NTRS2000

Cascade Optimization for Aircraft Engines With Regression and Neural Network Analysis - Approximators

The NASA Engine Performance Program (NEPP) can configure and analyze almost any type of gas turbine engine that can be generated through the interconnection of a set of standard physical components. In addition, the code can optimize engine performance by changing adjustable variables under a set of constraints. However, for engine cycle problems at certain operating points, the NEPP code can encounter difficulties: nonconvergence in the currently implemented Powell's optimization algorithm and deficiencies in the Newton-Raphson solver during engine balancing. A project was undertaken to correct these deficiencies. Nonconvergence was avoided through a cascade optimization strategy, and deficiencies associated with engine balancing were eliminated through neural network and linear regression methods. An approximation-interspersed cascade strategy was used to optimize the engine's operation over its flight envelope. Replacement of Powell's algorithm by the cascade strategy improved the optimization segment of the NEPP code. The performance of the linear regression and neural network methods as alternative engine analyzers was found to be satisfactory. This report considers two examples-a supersonic mixed-flow turbofan engine and a subsonic waverotor-topped engine-to illustrate the results, and it discusses insights gained from the improved version of the NEPP code.

Patnaik, Surya N.↗
NASA NTRS1971

The closure approximation in the hierarchy equations.

The expectation of the solution process in a stochastic operator equation can be obtained from averaged equations only under very special circumstances. Conditions for validity are given and the significance and validity of the approximation in widely used hierarchy methods and the ?self-consistent field' approximation in nonequilibrium statistical mechanics are clarified. The error at any level of the hierarchy can be given and can be avoided by the use of the iterative method.

Adomian, G.↗
NASA NTRS1995

Electronic Compensation For Distortion Of Antenna Reflector

Proposed method of obtaining approximately desired radiation or reception pattern from antenna that includes reflector based on concept of superposition of electromagnetic fields generated by multiple feedhorns or feed antenna elements arrayed at various positions near reflector and excited at electronically adjustable magnitudes and phases. In intended application, reflector nominally paraboloidal, feed elements N feedhorns in hexagonal array, and method used to compensate for deviations of real reflector surface from nominal paraboloidal shape. Method and concept also applicable to electronic beam steering and electronic antenna compensation in other situations.

Bailey, M. C.↗
DOE OSTI2025

A comparative study of calibration techniques for finite strain elastoplasticity: Numerically-exact sensitivities for FEMU and VFM

Accurate identification of material parameters is crucial for predictive modeling in computational mechanics. Here, the two primary approaches in the experimental mechanics community for calibration from full-field digital image correlation data are known as finite element model updating (FEMU) and the virtual fields method (VFM). In VFM, the objective function is a squared mismatch between internal and external virtual work or power. In FEMU, the objective function quantifies the weighted mismatch between model predictions and corresponding experimentally measured quantities of interest. It is minimized by iteratively updating the parameters of an FE model. While FEMU is seen as more flexible, VFM is commonly used instead of FEMU due to its considerably greater computational expense. However, comparisons between the two methods usually involve approximations of gradients or sensitivities with finite difference schemes, thereby making direct assessments difficult. Hence, in this study, we compare VFM and FEMU in the context of numerically-exact sensitivities obtained through local sensitivity analyses and the application of automatic differentiation software. To this end, we conduct a series of test cases to assess both methods under practical challenges using a finite strain elastoplasticity model.

Automatic differentiation↗
DOE OSTI2018

Estimating the joint distribution of rate parameters across multiple reactions in the absence of experimental data

Here, a procedure for determining the joint uncertainty of Arrhenius parameters across multiple combustion reactions of interest is demonstrated. This approach is capable of constructing the joint distribution of the Arrhenius parameters arising from the uncertain measurements performed in specific target experiments without having direct access to the underlying experimental data. The method involves constructing an ensemble of hypothetical data sets with summary statistics consistent with the available information reported by the experimentalists, followed by a fitting procedure that learns the structure of the joint parameter density across reactions using this consistent hypothetical data as evidence. The procedure is formalized in a Bayesian statistical framework, employing maximum-entropy and approximate Bayesian computation methods and utilizing efficient Markov chain Monte Carlo techniques to explore data and parameter spaces in a nested algorithm. We demonstrate the application of the method in the context of experiments designed to measure the rates of selected chain reactions in the H 2 -O 2 system and highlight the utility of this approach for revealing the critical correlations between the parameters within a single reaction and across reactions, as well as for maximizing consistency when utilizing rate parameter information in predictive combustion modeling of systems of interest.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Electronic structure and x-ray magnetic circular dichroism in YbAgCu4 and YbInCu4

X-ray absorption and magnetic circular dichroism spectra at the L-3 edge of Yb have been studied theoretically in the mixed-valent compounds YbAgCu4 and YbInCu4 using the fully relativistic Dirac linear muffin-tin orbital band structure method. The effect of the spin-orbit (SO) interaction and Coulomb repulsion U on the electronic structure of YbAgCu4 and YbInCu4 is examined in the frame of the local spin-density approximation (LSDA)+SO+U method. The core-hole effect in the final states as well as the effect of the electric quadrupole E-2 and magnetic dipole M-1 transitions have been investigated. The calculated results are compared with available experimental data.

Antonov, VN↗
DOE OSTI2014

A method for embedding circular force-free flux ropes in potential magnetic fields

We propose a method for constructing approximate force-free equilibria in pre-eruptive configurations in which a thin force-free flux rope is embedded into a locally bipolar-type potential magnetic field. The flux rope is assumed to have a circular-arc axis, a circular cross-section, and electric current that is either concentrated in a thin layer at the boundary of the rope or smoothly distributed across it with a maximum of the current density at the center. The entire solution is described in terms of the magnetic vector potential in order to facilitate the implementation of the method in numerical magnetohydrodynamic (MHD) codes that evolve the vector potential rather than the magnetic field itself. The parameters of the flux rope can be chosen so that its subsequent MHD relaxation under photospheric line-tied boundary conditions leads to nearly exact numerical equilibria. To show the capabilities of our method, we apply it to several cases with different ambient magnetic fields and internal flux-rope structures. These examples demonstrate that the proposed method is a useful tool for initializing data-driven simulations of solar eruptions.

79 ASTRONOMY AND ASTROPHYSICS↗