Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical approximation & analysis”

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

AutoBayes Program Synthesis System Users Manual

Program synthesis is the systematic, automatic construction of efficient executable code from high-level declarative specifications. AutoBayes is a fully automatic program synthesis system for the statistical data analysis domain; in particular, it solves parameter estimation problems. It has seen many successful applications at NASA and is currently being used, for example, to analyze simulation results for Orion. The input to AutoBayes is a concise description of a data analysis problem composed of a parameterized statistical model and a goal that is a probability term involving parameters and input data. The output is optimized and fully documented C/C++ code computing the values for those parameters that maximize the probability term. AutoBayes can solve many subproblems symbolically rather than having to rely on numeric approximation algorithms, thus yielding effective, efficient, and compact code. Statistical analysis is faster and more reliable, because effort can be focused on model development and validation rather than manual development of solution algorithms and code.

Schumann, Johann↗

Measuring sub-surface spatially varying thermal conductivity of silicon implanted with krypton

The thermal properties of semiconductors following exposure to ion irradiation are of great interest for the cooling of electronic devices; however, gradients in composition and structure due to irradiation often make the measurement difficult. Furthermore, the nature of spatial variations in thermal resistances due to spatially varying ion irradiation damage is not well understood. In this work, we develop an advancement in the analysis of time-domain thermoreflectance to account for spatially varying thermal conductivity in a material resulting from a spatial distribution of defects. We then use this method to measure the near-surface (≲1 μm) thermal conductivity of silicon wafers irradiated with Kr + ions, which has an approximate Gaussian distribution centered 260 nm into the sample. Our numerical analysis presented here allows for the spatial gradient of thermal conductivity to be extracted via what is fundamentally a volumetric measurement technique. We validate our findings via transmission electron microscopy, which is able to confirm the spatial variation of the sub-surface silicon structure, and provide additional insight into the local structure resulting from the effects of ion bombardment. Thermal measurements found the ion stopping region to have a nearly 50x reduction in thermal conductivity as compared to pristine silicon, while TEM showed the region was not fully amorphized. Our results suggest this drastic reduction in silicon thermal conductivity is primarily driven by structural defects in crystalline regions along with boundary scattering between amorphous and crystalline regions, with a negligible contribution being due to implanted krypton ions themselves.

36 MATERIALS SCIENCE↗

Exponential Runge-Kutta Parareal for non-diffusive equations

Parareal is a well-known parallel-in-time algorithm that combines a coarse and fine propagator within a parallel iteration. It allows for large-scale parallelism that leads to significantly reduced computational time compared to serial time-stepping methods. However, like many parallel-in-time methods it can fail to converge when applied to non-diffusive equations such as hyperbolic systems or dispersive nonlinear wave equations. Here, this paper explores the use of exponential integrators within the Parareal iteration. Exponential integrators are particularly interesting candidates for Parareal because of their ability to resolve fast-moving waves, even at the large stepsizes used by coarse propagators. This work begins with an introduction to exponential Parareal integrators followed by several motivating numerical experiments involving the nonlinear Schrödinger equation. These experiments are then analyzed using linear analysis that approximates the stability and convergence properties of the exponential Parareal iteration on nonlinear problems. The paper concludes with two additional numerical experiments involving the dispersive Kadomtsev-Petviashvili equation and the hyperbolic Vlasov-Poisson equation. These experiments demonstrate that exponential Parareal methods offer improved time-to-solution compared to serial exponential integrators when solving certain non-diffusive equations.

97 MATHEMATICS AND COMPUTING↗

Multilevel Monte Carlo Estimation of Unbiased Expectation via Sample Reuse and the Low Variance Estimation of Asymptotic Rates

A new variant of the multilevel Monte Carlo estimator [5, 3, 9, 12] is presented for the estimation of expectation statistics that utilizes sample reuse in specified levels, explicitly removes approximation error bias associated with numerically computed output quantities of interest that have an asymptotic limit behavior, and permits a low variance estimate of the asymptotic rate of convergence to that limit. In addition, it is shown that this new multilevel Monte Carlo variant can yield a computational cost savings. A review of Monte Carlo and multilevel Monte Carlo estimators is presented that includes analysis of expected value, expected mean squared error, and the calculation of optimized multilevel sample size parameters. The multilevel Monte Carlo estimator produces estimates of expectation for numerically approximated output quantities of interest that are biased by approximation error. When the quantity of interest can be modeled as the asymptotic limit of numerically approximated output quantities of interest, it is theoretically possible to remove this approximation error bias in the multilevel Monte Carlo estimator. In actual implementations, however, this procedure is unreliable due to statistical variability and inaccuracy in estimating the needed asymptotic limit. Analysis and numerical experiment show that the proposed variant of the multilevel Monte Carlo method greatly reduces (in some cases eliminates) the statistical variability in this limit estimation.

Barth, Timothy↗

Eigenvalue error analysis of viscously damped structures using a Ritz reduction method

The efficient solution of the eigenvalue problem that results from inserting passive dampers with variable stiffness and damping coefficients into a structure is addressed. Eigenanalysis of reduced models obtained by retaining a number of normal modes augmented with Ritz vectors corresponding to the static solutions resulting from the load patterns introduced by the dampers has been empirically shown to yield excellent approximations to the full eigenvalue problem. An analysis of this technique in the case of a single damper is presented. A priori and a posteriori error estimates are generated and tested on numerical examples. Comparison theorems with modally truncated models and a Markov parameter matching reduced-order model are derived. These theorems corroborate the heuristic that residual flexibility methods improve low-frequency approximation of the system. The analysis leads to other techniques for eigenvalue approximation. Approximate closed-form solutions are derived that include a refinement to eigenvalue derivative methods for approximation. An efficient Newton scheme is also developed. A numerical example is presented demonstrating the effectiveness of each of these methods.

Chu, Cheng-Chih↗

Impulsive phase of solar flares. I - Characteristics of high energy electrons

The models describing the operative physical processes occurring during the impulsive phase of a solar flare can be broadly classified as thermal or nonthermal models. The current investigation is part of a series of studies which analyze the general characteristics of the 'nonthermal' models. The results are presented of an analysis regarding the steady state distribution of energetic electrons within the flare plasma. The time independent Fokker-Planck equation is used to determine the steady state distribution of the electrons in the ambient plasma. The model investigated encompasses thick-target, thin-target, and trap aspects of nonthermal hard X-ray models. High energy electrons with some initial energy spectrum and pitch angle distribution are injected at a height above the photosphere either into a closed magnetic loop structure or into an open field structure. Some approximate analytic solutions are discussed and the results of a numerical analysis are presented.

Leach, J.↗

Learning continuous models for continuous physics

Abstract Dynamical systems that evolve continuously over time are ubiquitous throughout science and engineering. Machine learning (ML) provides data-driven approaches to model and predict the dynamics of such systems. A core issue with this approach is that ML models are typically trained on discrete data, using ML methodologies that are not aware of underlying continuity properties. This results in models that often do not capture any underlying continuous dynamics—either of the system of interest, or indeed of any related system. To address this challenge, we develop a convergence test based on numerical analysis theory. Our test verifies whether a model has learned a function that accurately approximates an underlying continuous dynamics. Models that fail this test fail to capture relevant dynamics, rendering them of limited utility for many scientific prediction tasks; while models that pass this test enable both better interpolation and better extrapolation in multiple ways. Our results illustrate how principled numerical analysis methods can be coupled with existing ML training/testing methodologies to validate models for science and engineering applications.

97 MATHEMATICS AND COMPUTING↗

Composite box beam analysis - Theory and experiments

Beam theory is widely used as a first approximation in numerous structural applications. When applied to composite beams, the accuracy of beam theory becomes questionable because (1) the shearing and warping deformations become significant, as the shearing stiffness of composite laminates is often very low, and (2) several elastic couplings can occur that strongly influence the behavior of composite beams. The torsional behavior of thin-walled composite beams has important implications for aeronautical structures and is deeply modified by the above nonclassical effects. This paper presents two comprehensive analysis methodologies for composite beams and describes experimental results obtained from a thin-walled, rectangular cross-sectional beam. The theoretical predictions are found in good agreement with the observed twist and strain distributions. Out-of-plane torsional warping of the cross-section is found to be the key factor for an accurate modeling of the torsional behavior of such structures.

Bauchau, O. A.↗

Gas Evolution Dynamics in Godunov-Type Schemes and Analysis of Numerical Shock Instability

In this paper we are going to study the gas evolution dynamics of the exact and approximate Riemann solvers, e.g., the Flux Vector Splitting (FVS) and the Flux Difference Splitting (FDS) schemes. Since the FVS scheme and the Kinetic Flux Vector Splitting (KFVS) scheme have the same physical mechanism and similar flux function, based on the analysis of the discretized KFVS scheme the weakness and advantage of the FVS scheme are closely observed. The subtle dissipative mechanism of the Godunov method in the 2D case is also analyzed, and the physical reason for shock instability, i.e., carbuncle phenomena and odd-even decoupling, is presented.

Xu, Kun↗

Comment on high resolution simulations of cosmic strings. 1: Network evoloution

Comments are made on recent claims (Albrecht and Turok, 1989) regarding simulations of cosmic string evolution. Specially, it was claimed that results were dominated by a numerical artifact which rounds out kinks on a scale of the order of the correlation length on the network. This claim was based on an approximate analysis of an interpolation equation which is solved herein. The typical rounding scale is actually less than one fifth of the correlation length, and comparable with other numerical cutoffs. Results confirm previous estimates of numerical uncertainties, and show that the approximations poorly represent the real solutions to the interpolation equation.

Turok, Neil↗

Optical Time Transfer for Bistatic SAR Spacecraft

A spacecraft-to-spacecraft optical time-transfer simulation has been developed as a tool for informing NASA’s Surface Deformation and Change (SDC) mission architecture. The SDC mission will combine radar images from multiple spacecraft to improve understanding of the Earth’s sea-level and landscape changes. Spacecraft must be precisely synchronized in order to create sharp radar images. Simulation of multiple spacecraft time-synchronizing via laser communication can inform technology choices of a mission by providing a picosecond-precision level estimate of clock error. This timing and ranging simulation has been combined with a radar system performance analysis pipeline. The simulated timing errors are used to predict performance of bistatic SAR systems in the presence of oscillator noise and time synchronization in accuracy. This analysis includes both analytic approximation equations from existing literature, and a numerical radar simulation to extract key system performance parameters like phase error and signal-to-noise ratio (SNR)degradation. Precision time-transfer techniques facilitate the accurate synchronization of clocks between any combination of terminals. Most time-transfer technology for comparing two clocks at different terminals use radio frequencies (RF) to measure the time delay between the sending and receiving of signals. Laser technology offers the capability to transmit high data rates with systems that are of smaller size and lower power than comparable RF systems. The clocks on independent spacecraft will have some phase and frequency errors between them that result in clock drift. The two clock models that are included in this bi-directional MATLAB simulation are a cesium-based Chip-Scale Atomic Clock (CSAC) and a rubidium-based Miniature Atomic Clock (MAC). The CSAC has flown as hardware for small satellite missions such as the University of Florida’s CHOMPTT mission. A study of example orbits, including that of NASA NASA-ISRO Synthetic Aperture Radar Mission (NISAR) mission, and lasing rates demonstrate the impact of flight configuration parameters on the synchronization error between two spacecraft. The MATLAB timing simulation uses a Runge-Kutta 4th-order method to propagate spacecraft orbits and computes the light-travel time estimate between them. The simulation outputs the estimated range and estimated clock error based on a user-defined spacecraft cluster configuration. The radar simulation and analytic approximations are applied to evaluate a potential future NASA bistatic SAR constellation architecture. In the proposed architecture, satellites follow each other in the same orbit at 800 km altitude, with a 210 km baseline. We also baseline the CSAC as an ultra stable oscillator, and use NASA’s NISAR for baseline radar system parameters to compute a clock-system introduced phase error of 5.6 degrees without synchronization by frequent time transfer. We build on this base case with a sensitivity analysis of radar performance over a proposed range of constellation and radar system parameters. With this analysis pipeline, we comment on which radar parameters should or should not be changed to minimize synchronization requirements. This analysis technique could be extended or modified to evaluate the timing requirements of other geometries for other future multistatic SAR missions, or other interferometric satellite missions.

Surface Deformation and Change↗

A symmetric modal formulation of fluid-structure interaction, including a static approximation to higher order fluid modes

A method of fluid-structure coupling which provides symmetrical matrix equations of standard form solved by existing finite element computer programs is presented. The method postulates that the uncoupled vibration modes of the fluid or the structure be calculated before the coupled analysis. A numerical solution of vibration modes in an axisymmetric container demonstrated that a static approximation to higher order fluid modes can improve the accuracy of dynamic response computations using modal methods.

Macneal, R. H.↗

Interpolation Method Needed for Numerical Uncertainty Analysis of Computational Fluid Dynamics

Using Computational Fluid Dynamics (CFD) to predict a flow field is an approximation to the exact problem and uncertainties exist. There is a method to approximate the errors in CFD via Richardson's Extrapolation. This method is based off of progressive grid refinement. To estimate the errors in an unstructured grid, the analyst must interpolate between at least three grids. This paper describes a study to find an appropriate interpolation scheme that can be used in Richardson's extrapolation or other uncertainty method to approximate errors. Nomenclature

CFD↗

Elimination of Gibbs' phenomena from error analysis of finite element results

This paper is one of a series on error analysis and correction of finite element solutions for plates and shells. The error analysis in the earlier papers used half-range double Fourier sine series for numerical harmonic analysis. The half-range formulas are simple to apply, but they can be inaccurate near the ends of the ranges of the independent variables. The Gibbs' phenomenon exhibited by half-range sine series in one independent variable has a two-dimensional analog; a classic example is the Navier solution in a double half-range sine series for the simply supported plate under a uniform load. A simple change of variables is introduced in the paper to improve the accuracy of the double Fourier sine series without adding complexity to the numerical analysis. The change of variables is applied to the problem of approximating a transverse load that is tabulated on a rectangular grid. A solution based on the change of variables is compared with results from the Navier solution for the simply supported plate problem and finite element results for the same problem.

Thurston, Gaylen A.↗

Loop current eddy shedding estimated using Geosat altimeter data

A method for extracting periodic signals from altimeter data is applied in the Gulf of Mexico using data from the first two years of the Geosat Exact Repeat Mission. The least squares fitting of selected frequencies of sine and cosine functions to the time series of sea surface height at each point of the Geosat ground track is used to estimate the dominant frequency in the gulf for this time period. Optimal interpolation of the sine and cosine coefficients at the dominant frequency to a uniformly spaced grid allows reconstruction of the time series in the gulf and shows eddies shedding from the Loop Current with a period of approximately 10.5 months, consistent with shedding periods determined from numerical models and analysis of satellite infrared data. This technique is a novel method for investigation of periodic signals in Geosat altimetry and is not severely affected by the data outages present in the time series.

Jacobs, Gregg A.↗

Recent developments in equivalent plate modeling for wing shape optimization

A new technique for structural modeling of airplane wings is presented taking transverse shear effects into account. The kinematic assumptions of first order shear deformation plate theory in combination with numerical analysis based on simple polynomials which define geometry, construction and displacement approximations lead to analytical expressions for elements of the stiffness and mass matrices and load vector. Contributions from the cover skins, spar and rib caps and spar and rib webs are included as well as concentrated springs and concentrated masses. Limitations of current equivalent plate wing modeling techniques based on classical plate theory are discussed, and the improved accuracy of the new equivalent plate technique is demonstrated through comparison to finite element analysis and test results. Analytical derivatives of stiffness, mass and load terms with respect to wing shape lead to analytic sensitivities of displacements, stresses and natural modes with respect to planform shape and depth distribution. This makes the new capability an effective structural tool for wing shape optimization.

Livne, Eli↗