Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “linear extrapolation 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 73 records · Page 4

Verification of Data-Driven Models of Physical Phenomena using Interpretable Approximation

Machine-learned models, specifically neural networks, are increasingly used as “closures” or “constitutive models” in engineering simulators to represent fine-scale physical phenomena that are too computationally expensive to resolve explicitly. However, these neural net models of unresolved physical phenomena tend to fail unpredictably and are therefore not used in mission-critical simulations. In this report, we describe new methods to authenticate them, i.e., to determine the (physical) information content of their training datasets, qualify the scenarios where they may be used and to verify that the neural net, as trained, adhere to physics theory. We demonstrate these methods with neural net closure of turbulent phenomena used in Reynolds Averaged Navier-Stokes equations. We show the types of turbulent physics extant in our training datasets, and, using a test flow of an impinging jet, identify the exact locations where the neural network would be extrapolating i.e., where it would be used outside the feature-space where it was trained. Using Generalized Linear Mixed Models, we also generate explanations of the neural net (à la Local Interpretable Model agnostic Explanations) at prototypes placed in the training data and compare them with approximate analytical models from turbulence theory. Finally, we verify our findings by reproducing them using two different methods.

42 ENGINEERING↗

Approximate Green’s Function Coupled Cluster A Method Employing Effective Dimension Reduction

The Green’s function coupled cluster (GFCC) method, originally proposed in the early 1990s, is a powerful many-body tool for computing and analyzing the electronic structures of molecular structures and periodic systems, especially when electrons are strongly correlated. However, in order for the GFCC method to become an approach that can be routinely used in the electronic structure calculations, more robust linear algebraic techniques and approximations need to be employed to reduce its extremely high computational overhead. In our recent studies, we showed that the GFCC equation can be numerically solved by seeking solutions of multiple linear systems in the frequency regime of interest, which can be easily distributed in a massively parallel environment, and broaden the avenue for more state-of-the-art numerical linear algebraic tools associated with multiple linear systems to be kicked in to this field. In the present pursuit of this practice, we show a successful application of model order reduction (MOR) techniques in the GFCC framework. Briefly speaking, for a frequency regime of interest that requires a high resolution of spectral function, instead of solving GFCC linear equation in the original dimension for every single frequency point, an easily-solvable linear system model with reduced dimension can be built upon projecting the original GFCC linear system into a subspace, from which the interpolation and extrapolation can give a reasonable approximation to the results of the original GFCC linear equations. Here, we show that the subspace can be properly constructed in an iterative manner from the auxiliary vectors of the GFCC linear equations at some select frequencies inside the frequency regime. During the iterations, the quality of the subspace, as well as the linear system model, can be systematically improved. The approach is tested in this work in terms of the efficiency and accuracy of computing spectral functions for some typical molecular systems including carbon monoxide, 1,3-butadiene, benzene, and adenine molecules. As a byproduct, the high quality linear system model is also found to be able to provide high quality initial guess for the existing linear solver to improve the convergence performance.

Peng, Bo↗

Cosmological forecasts for future galaxy surveys with the linear point standard ruler: Toward consistent BAO analyses far from a fiducial cosmology

The linear point (LP) standard ruler was identified as the basis of a purely geometric method for exploiting the baryon acoustic oscillations (BAO). The LP exploits the BAO feature imprinted in the galaxy two-point correlation function to measure cosmological distances independent of any specific cosmological model. We forecast the expected precision of future and ongoing spectroscopic galaxy surveys to measure distances leveraging the linear point. We investigate the cosmological implications of our forecasted results. We focus in particular on a relevant working example: the detection of the late-time cosmic acceleration independent of other cosmological probes. Further, our findings show that, even within the Λ⁢CDM standard cosmological paradigm, estimated distances need to be reliable over a very wide parameter range in order to realize their maximum utility. This is particularly relevant if we aim to properly characterize cosmological tensions. The LP is a promising candidate approach to achieve this reliability. In contrast, widely employed procedures in BAO analysis estimate distances keeping fixed cosmological parameters to fiducial values close to cosmic-microwave-background constraints in flat-Λ⁢CDM. It is unclear whether they are purely geometric methods. Moreover, they rely on untested extrapolations to explore the parameter space away from those fiducial flat-Λ⁢CDM values. We recommend that all BAO methodologies be validated across the full range of models and parameters over which their results are quoted, first by means of linear predictions and then N -body simulations.

79 ASTRONOMY AND ASTROPHYSICS↗

Investigation of the mechanical and microstructural evolution of a Cu based bulk metallic glass during ion irradiation

Ion irradiation and annealing experiments were performed on Cu 60 Zr 20 Hf 10 Ti 10 bulk metallic glass (BMG) specimens to investigate their irradiation- and temperature-induced microstructural and mechanical property evolution. For the ion irradiations, samples were exposed to 9 MeV Ni 3+ ions to a midrange (~1.2 μm depth) dose of 10 displacements per atom (dpa) at temperatures ranging from room temperature to 360 °C (the corresponding peak dose at ~2.8 μm depth was ~25 dpa). Bulk X-ray diffraction (XRD) and transmission electron microscopy (TEM) characterization revealed that the alloy did not crystallize during irradiation up to 290 °C but did partially crystallize at 360 °C. XRD analysis revealed that the crystallization which occurred in the sample irradiated at 360 °C was caused by thermal effects instead of irradiation displacement damage. Subsequent Rietveld refinement analysis of the XRD measurements revealed the presence of two distinct crystal phases, namely a CuTiZr hexagonal structure belonging to the P63/mmc space group and a CuTi tetragonal structure belonging to the P4/mmm space group. Nanoindentation experiments revealed that no pronounced hardness changes occurred in the specimens irradiated at room temperature and 290°C, although significant hardening was observed in the sample irradiated at 360 °C. The significant increase in the hardness at 360°C was ascribed to thermally induced partial crystallization of the alloy instead of the ion irradiation. In general, the results of the nanoindentation experiments and XRD characterization suggest that although the Cu BMG exhibits good stability during irradiation at temperatures up to 290 °C it is not suitable for irradiation environments where the temperature is 360 °C for extended periods of time. The Lam and Chong extrapolation method, which has been used to study the indentation size effect (ISE) in amorphous alloys, was employed to quantify how irradiation and temperature affect this type of behavior in the BMG. However, the poor linear fitting of the indentation hardness data by this model indicate that a new ISE model is likely needed to quantify indentation hardening in BMGs.

36 MATERIALS SCIENCE↗

The radiation environment of OSO missions from 1974 to 1978

Trapped particle radiation levels on several OSO missions were calculated for nominal trajectories using improved computational methods and new electron environment models. Temporal variations of the electron fluxes were considered and partially accounted for. Magnetic field calculations were performed with a current field model and extrapolated to a later epoch with linear time terms. Orbital flux integration results, which are presented in graphical and tabular form, are analyzed, explained, and discussed.

Stassinopoulos, E. G.↗

Third-order multi-dimensional Euler/Navier-Stokes solver

A nonoscillatory multidimensional advection code is presented which is applicable to the study of very high Reynolds number steady flow, and results are obtained for the oblique advection of a step profile. The scheme is based on an exponential upwinding or linear extrapolation refinement of a previous QUICK scheme, allowing sharp monotonic resolution of the step while retaining third-order accuracy. The resulting simple flux-limiting criteria have been applied to a higher order method and a formally third-order supercompressive algorithm. Results are presented for first-order, second-order, and third-order upwinding.

Leonard, B. P.↗

Prediction of Undsteady Flows in Turbomachinery Using the Linearized Euler Equations on Deforming Grids

A linearized Euler solver for calculating unsteady flows in turbomachinery blade rows due to both incident gusts and blade motion is presented. The model accounts for blade loading, blade geometry, shock motion, and wake motion. Assuming that the unsteadiness in the flow is small relative to the nonlinear mean solution, the unsteady Euler equations can be linearized about the mean flow. This yields a set of linear variable coefficient equations that describe the small amplitude harmonic motion of the fluid. These linear equations are then discretized on a computational grid and solved using standard numerical techniques. For transonic flows, however, one must use a linear discretization which is a conservative linearization of the non-linear discretized Euler equations to ensure that shock impulse loads are accurately captured. Other important features of this analysis include a continuously deforming grid which eliminates extrapolation errors and hence, increases accuracy, and a new numerically exact, nonreflecting far-field boundary condition treatment based on an eigenanalysis of the discretized equations. Computational results are presented which demonstrate the computational accuracy and efficiency of the method and demonstrate the effectiveness of the deforming grid, far-field nonreflecting boundary conditions, and shock capturing techniques. A comparison of the present unsteady flow predictions to other numerical, semi-analytical, and experimental methods shows excellent agreement. In addition, the linearized Euler method presented requires one or two orders-of-magnitude less computational time than traditional time marching techniques making the present method a viable design tool for aeroelastic analyses.

Clark, William S.↗

Deep material network via a quilting strategy: visualization for explainability and recursive training for improved accuracy

Recent developments integrating micromechanics and neural networks offer promising paths for rapid predictions of the response of heterogeneous materials with similar accuracy as direct numerical simulations. The deep material network is one such approaches, featuring a multi-layer network and micromechanics building blocks trained on anisotropic linear elastic properties. Once trained, the network acts as a reduced-order model, which can extrapolate the material’s behavior to more general constitutive laws, including nonlinear behaviors, without the need to be retrained. However, current training methods initialize network parameters randomly, incurring inevitable training and calibration errors. Here, we introduce a way to visualize the network parameters as an analogous unit cell and use this visualization to “quilt” patches of shallower networks to initialize deeper networks for a recursive training strategy. The result is an improvement in the accuracy and calibration performance of the network and an intuitive visual representation of the network for better explainability.

97 MATHEMATICS AND COMPUTING↗

Development of a Supercharged Octane Number and a Supercharged Octane Index

Gasoline knock resistance is characterized by the Research and Motor Octane Number (RON and MON), which are rated on the CFR octane rating engine at naturally aspirated conditions. However, modern automotive downsized boosted spark ignition (SI) engines generally operate at higher cylinder pressures and lower temperatures relative to the RON and MON tests. Using the naturally aspirated RON and MON ratings, the octane index (OI) characterizes the knock resistance of gasolines under boosted operation by linearly extrapolating into boosted “beyond RON” conditions via RON, MON, and a linear regression K factor. Using OI solely based on naturally aspirated RON and MON tests to extrapolate into boosted conditions can lead to significant errors in predicting boosted knock resistance between gasolines due to non-linear changes in autoignition and knocking characteristics with increasing pressure conditions. Here, a new “Supercharged Octane Number” (SON) method was developed on the CFR engine at increased intake pressures, which improved the correlation to boosted knock-limited automotive SI engine data over RON for several surrogate fuels and gasolines, including five “Co-Optima” RON 98 fuels and an E10 regular grade gasoline. Furthermore, the conventional OI was extended to a newly introduced Supercharged Octane Index (OI S ) based on SON and RON, which significantly improved the correlation to fuel knock resistance measurements from modern boosted SI engine knock-limited spark advance tests. This demonstrated the first proof of concept of a SON and OI S to better characterize a fuel’s knock resistance in modern boosted SI engines.

42 ENGINEERING↗

Evaluation of algorithms for geological thermal-inertia mapping

The errors incurred in producing a thermal inertia map are of three general types: measurement, analysis, and model simplification. To emphasize the geophysical relevance of these errors, they were expressed in terms of uncertainty in thermal inertia and compared with the thermal inertia values of geologic materials. Thus the applications and practical limitations of the technique were illustrated. All errors were calculated using the parameter values appropriate to a site at the Raft River, Id. Although these error values serve to illustrate the magnitudes that can be expected from the three general types of errors, extrapolation to other sites should be done using parameter values particular to the area. Three surface temperature algorithms were evaluated: linear Fourier series, finite difference, and Laplace transform. In terms of resulting errors in thermal inertia, the Laplace transform method is the most accurate (260 TIU), the forward finite difference method is intermediate (300 TIU), and the linear Fourier series method the least accurate (460 TIU).

Miller, S. H.↗

Sonic Boom Prediction and Minimization of the Douglas Reference OPT5 Configuration

Conventional CFD methods and grids do not yield adequate resolution of the complex shock flow pattern generated by a real aircraft geometry. As a result, a unique grid topology and supersonic flow solver was developed at Northrop Grumman based on the characteristic behavior of supersonic wave patterns emanating from the aircraft. Using this approach, it was possible to compute flow fields with adequate resolution several body lengths below the aircraft. In this region, three-dimensional effects are diminished and conventional two-dimensional modified linear theory (MLT) can be applied to estimate ground pressure signatures or sonic booms. To accommodate real aircraft geometries and alleviate the burdensome grid generation task, an implicit marching multi-block, multi-grid finite-volume Euler code was developed as the basis for the sonic boom prediction methodology. The Thomas two-dimensional extrapolation method is built into the Euler code so that ground signatures can be obtained quickly and efficiently with minimum computational effort suitable to the aircraft design environment. The loudness levels of these signatures can then be determined using a NASA generated noise code. Since the Euler code is a three-dimensional flow field solver, the complete circumferential region below the aircraft is computed. The extrapolation of all this field data from a cylinder of constant radius leads to the definition of the entire boom corridor occurring directly below and off to the side of the aircraft's flight path yielding an estimate for the entire noise "annoyance" corridor in miles as well as its magnitude. An automated multidisciplinary sonic boom design optimization software system was developed during the latter part of HSR Phase 1. Using this system, it was found that sonic boom signatures could be reduced through optimization of a variety of geometric aircraft parameters. This system uses a gradient based nonlinear optimizer as the driver in conjunction with a computationally efficient Euler CFD solver (NIIM3DSB) for computing the three-dimensional near-field characteristics of the aircraft. The intent of the design system is to identify and optimize geometric design variables that have a beneficial impact on the ground sonic boom. The system uses a simple wave drag data format to specify the aircraft geometry. The geometry is internally enhanced and analytic methods are used to generate marching grids suitable for the multi-block Euler solver. The Thomas extrapolation method is integrated into this system, and hence, the aircraft's centerline ground sonic boom signature is also automatically computed for a specified cruise altitude and yields the parameters necessary to evaluate the design function. The entire design system has been automated since the gradient based optimization software requires many flow analyses in order to obtain the required sensitivity derivatives for each design variable in order to converge on an optimal solution. Hence, once the problem is defined which includes defining the objective function and geometric and aerodynamic constraints, the system will automatically regenerate the perturbed geometry, the necessary grids, the Euler solution, and finally the ground sonic boom signature at the request of the optimizer.

Siclari, Michael J.↗

Microbial community structure at the U.S.-Joint Global Ocean Flux Study Station ALOHA: Inverse methods for estimating biochemical indicator ratios

Modeling biogeochemical fluxes in the marine plankton requires the application of factors for extrapolation of biomass indicators measured in the field (chlorophyll a, adenosine triphosphate, bacterial counts) to biomass carbon or nitrogen. These are often inferred from culture studies and are poorly constrained for natural populations. At least squares inverse method with a simple linear model constrains the values of several common indicator ratios, giving self-consistent solutions that provide useful information about the structure of the microbial community at our North Pacific Ocean study site (Station ALOHA (A Long-term Oligotrophic Habitat Assessment)). These results indicate that the fraction of the microbial biomass that is autotrophic (pigmented) is greater in the mixed layer than at the deep chlorophyll maximum layer and that heterotrophic bacteria are a significant but not necessarily predominant component of the microbial community in the euphotic zone.

Christian, James R.↗

Predictive aging of polymers

A method of predicting aging of polymers operates by heating a polymer in the outdoors to an elevated temperature until a change of property is induced. The test is conducted at a plurality of temperatures to establish a linear Arrhenius plot which is extrapolated to predict the induction period for failure of the polymer at ambient temperature. An Outdoor Photo Thermal Aging Reactor (OPTAR) is also described including a heatable platen for receiving a sheet of polymer, means to heat the platen, and switching means such as a photoelectric switch for turning off the heater during dark periods.

Cuddihy, Edward F.↗

Predictive aging of polymers

A method of predicting aging of polymers operates by heating a polymer in the outdoors to an elevated temperature until a change of property is induced. The test is conducted at a plurality of temperatures to establish a linear Arrhenius plot which is extrapolated to predict the induction period for failure of the polymer at ambient temperature. An Outdoor Photo Thermal Aging Reactor (OPTAR) is also described including a heatable platen for receiving a sheet of polymer, means to heat the platen and switching means such as a photoelectric switch for turning off the heater during dark periods.

Cuddihy, Edward F.↗

Machine learning a molecular Hamiltonian for predicting electron dynamics

We develop a computational method to learn a molecular Hamiltonian matrix from matrix-valued time series of the electron density. As we demonstrate for three small molecules, the resulting Hamiltonians can be used for electron density evolution, producing highly accurate results even when propagating 1,000 time steps beyond the training data. As a more rigorous test, we use the learned Hamiltonians to simulate electron dynamics in the presence of an applied electric field, extrapolating to a problem that is beyond the field-free training data. We find that the resulting electron dynamics predicted by our learned Hamiltonian are in close quantitative agreement with the ground truth. Our method relies on combining a reduced-dimensional, linear statistical model of the Hamiltonian with a time-discretization of the quantum Liouville equation within time-dependent Hartree Fock theory. In conclusion, we train the model using a least-squares solver, avoiding numerous, CPU-intensive optimization steps. For both field-free and field-on problems, we quantify training and propagation errors, highlighting areas for future development.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Strain-Gage Loads Calibration Parametric Study

This paper documents a parametric study of various aircraft wing-load test features that affect the quality of the resultant derived shear, bending-moment, and torque strain-gage load equations. The effect of the following on derived strain-gage equation accuracy are compared: single-point loading compared with distributed loading, variation in applied test load magnitude, number of applied load cases, and wing-box-only compared with control-surface loading. The subject of this study is an extensive wing-load calibration test of the Active Aeroelastic Wing F/A-18 airplane. Selected subsets of the available test data were used to derive load equations using the linear regression method. Results show the benefit of distributed loading and the diminishing-return benefits of test load magnitudes and number of load cases. The use of independent check cases as a quality metric for the derived load equations is shown to overcome blind extrapolating beyond the load data used to derive the load equations.

Lokos, William A.↗

Calculation of unsteady linearized Euler flows in cascades using harmonically deforming grids

A method for calculating unsteady, inviscid, compressible flows in cascades is presented. Using the linearized Euler technique, the flow is decomposed into a steady or mean flow plus a harmonically varying small disturbance flow. The equations that describe the small disturbance flow are linear variable coefficient equations, and are solved using a pseudo-time time marching Lax-Wendroff technique. Unlike previous linearized methods, however, the solution is computed on a harmonically deforming computational grid that conforms to the motion of the vibrating airfoils. The mean flow and perturbation flow solutions are defined in the deforming coordinate system rather than in a coordinate system fixed in space. Hence, no extrapolation terms are required to implement the upwash boundary conditions at the airfoil surfaces, significantly improving the accuracy of the method. For transonic flow calculations, unsteady shock motions are modelled using shock capturing. The unsteady loads due to the shock motion are then seen as pressure impulses. Representative computational results are presented for transonic channel flows and subsonic and transonic cascade flows.

Hall, Kenneth C.↗

Some observations on the Houbolt-Rainey and peak-hold methods of flutter onset prediction

A subcritical response method for flutter onset prediction developed by Houbolt and Rainey in 1958 is compared with the Peak-Hold Method which was apparently first applied to flutter onset prediction by Sandford, Abel, and Gray in the early 1970's. The rational argument presented shows that the two methods are not different, but are actually the same. So, because there is an analytical foundation for the Houbolt-Rainey Method, then there is the same analytical foundation for the Peak-Hold Method. Further, it is suggested that, in applying Peak-Hold Method in cases where turbulence is used as the excitation force, the variation of the reciprocal of the response amplitude with the reciprocal of the dynamic pressure to be used to extrapolate to flutter onset rather than the variation with dynamic pressure which is the current practice because the linear trend which is predicted to occur for the former is easier to extrapolate to the flutter condition than the nonlinear trend predicted to occur for the latter.

Doggett, Robert V., Jr.↗