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 109 records · Page 6

An analysis method for two-dimensional transonic viscous flow

A method for the approximate calculation of transonic flow over airfoils, including shock waves and viscous effects, is described. Numerical solutions are obtained by use of a computer program which is discussed in the appendix. The importance of including the boundary layer in the analysis is clearly demonstrated, as well as the need to improve on existing procedures near the trailing edge. Comparisons between calculations and experimental data are presented for both conventional and supercritical airfoils, emphasis being on the surface pressure distribution, and good agreement is indicated.

Bavitz, P. C.↗

On an integro-differential equation model for the study of the response of an acoustically coupled panel

The response of a clamped panel to supersonically convected turbulence is considered. A theoretical model in the form of an integro-differential equation is employed that takes into account the coupling between the panel motion and the surrounding acoustic medium. The kernels of the integrals, which represent induced pressures due to the panel motion, are Green's functions for sound radiations under various moving and stationary sources. An approximate analysis is made by following a finite-element Ritz-Galerkin procedure. Preliminary numerical results, in agreement with experimental findings, indicate that the acoustic damping is the controlling mechanism of the response.

Yen, D. H. Y.↗

The effect of torsional flexibility on the rolling characteristics at supersonic speeds of tapered unswept wings

An analysis is presented of the effect of torsional flexibility on the rolling characteristics at supersonic speeds of tapered unswept wings with partial-span constant-percent-chord ailerons extending inboard from the wing tip. The geometric variables considered are aspect ratio, taper ratio, aileron span, and aileron chord. The shape of the wing-torsional-stiffness curve is assumed and the twisting moment is considered to result solely from the pressure distribution caused by aileron deflection, so that the necessity of using a successive-approximation method is avoided. Because of the complexity of the equations resulting from the analysis, numerical calculations from the equations are presented in a series of figures. A computational form is provided to be used in conjunction with these figures so that calculations can be made without reference to the analysis.

Tucker, Warren A↗

RandONets: Shallow networks with random projections for learning linear and nonlinear operators

Deep neural networks have been extensively used for the solution of both the forward and the inverse problem for dynamical systems. However, their implementation necessitates optimizing a high-dimensional space of parameters and hyperparameters. This fact, along with the requirement of substantial computational resources, pose a barrier to achieving high numerical accuracy, but also interpretability. Here, to address the above challenges, we present Random Projection-based Operator Networks (RandONets): shallow networks with random projections and tailor-made numerical analysis methods that learn accurately and fast linear and nonlinear operators. Building on previous works, we prove that RandOnets are universal approximators of linear and nonlinear operators. Due to their simplicity, RandONets provide a one-step transformation of the input space, facilitating interpretability. For the evaluation of their performance, we focus on operators of PDEs. We show, that RandONets outperform by several orders of magnitude, both in terms of numerical approximation accuracy and computational cost, the “vanilla” DeepONets. Hence, we believe that our method will trigger further developments in the field of scientific machine learning, for the development of new ‘’light”schemes that will provide high accuracy while reducing dramatically the computational cost. A MATLAB toolbox for RandONets, including demos, is available on GitHub at https://github.com/GianlucaFabiani/RandONets.

Interpretable machine learning↗

Conservative numerical schemes with optimal dispersive wave relations: Part I. Derivation and analysis

An energy-conserving and an energy-and-enstrophy conserving numerical schemes are derived by approximating the Hamiltonian formulation of the inviscid shallow water flows based on the vorticity-divergence variables. These schemes also conserve the first-order moments such as mass and vorticity, as usual. The conservative properties of the schemes stem from the skew-symmetry and singularities of the Poisson brackets, which are carefully retained in the discrete approximations. Here, the schemes operate on unstructured orthogonal dual meshes, over bounded or unbounded domains, and they are also shown to possess the same optimal dispersive wave relations as those of the Z-grid scheme, which is a consequence of the use of the vorticity and divergence variables.

54 ENVIRONMENTAL SCIENCES↗

Error analysis of analytic solutions for self-excited near-symmetric rigid bodies - A numerical study

Analytic error bounds are presented for the solutions of approximate models for self-excited near-symmetric rigid bodies. The error bounds are developed for analytic solutions to Euler's equations of motion. The results are applied to obtain a simplified analytic solution for Eulerian rates and angles. The results of a sample application of the range and error bound expressions for the case of the Galileo spacecraft experiencing transverse torques demonstrate the use of the bounds in analyses of rigid body spin change maneuvers.

Kia, T.↗

A charged particle transport approximation for thick and thin plasmas

In this work, we will construct a simple method to be utilized in the interpretation of the properties of a thin plasma from a study of the transport of charged particles through it. We do so by first demonstrating how charged particle fluxes encode local plasma information in their spectra for sufficiently thick plasmas and then propose how to extend that analysis to the thin plasma limit where the extensive geometry cannot be neglected through arguments from scale separation. We provide a numerical treatment of the transport problem and demonstrate that it is in good agreement with the approximation we present. Finally, we utilize both the numerical solution and our novel approximation to study the impact of the extensive scale of the plasma on the yield and flux normalized high-energy neutron spectra resulting from the upscattering of charged fuel ions. Using this analysis, we show that (after controlling for fusion yield) for the same uniform densities and temperatures, larger plasmas have higher magnitude but softer reaction-in-flight (RIF) neutron spectra relative to smaller plasmas. This is because larger plasmas retain more knocked-on suprathermal ions within their bulk and can downscatter them to lower average energies, while smaller plasmas allow a larger fraction of high energy particles to “range out” of the system prior to substantial downscattering or inducing an RIF reaction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Stable boundary approximations for a class of implicit schemes for the one-dimensional inviscid equations of gas dynamics

The applicability to practical calculations of recent theoretical developments in the stability analysis of difference approximations for initial-boundary-value problems of the hyperbolic type. For the numerical experiments, select the one-dimensional inviscid gas-dynamic equations in conservation-law form is selected. A class of implicit schemes based on linear multistep methods for ordinary differential equations is chosen and the use of space or space-time extrapolations as implicit or explicit boundary schemes is emphasized. Some numerical examples with various inflow-outflow conditions highlight the commonly discussed issues: explicit versus implicit boundary schemes, unconditionally stable schemes, and underspecification or overspecification of boundary conditions.

Yee, H. C.↗

The Analysis of Turbulent Flow by Hot Wire Signals

When measuring velocities in turbulent gas flow, approximation signal analysis with hot wire anemometers having one and two wire probes are used. A numeric test of standard analyses shows the resulting systemmatic error increases quickly with increasing turbulent intensity. Since it also depends on the turbulence structure, it cannot be corrected. The use of such probes is thus restricted to low turbulence. By means of three wire probes (in two dimensional flows with X wire probes) in principle, instantaneous values of velocity can be determined, and an asymmetric arrangement of wires has a theoretical advantage.

Bartenwerfer, M.↗

Flow Effects on the Flammability Diagrams of Solid Fuels

A research program is currently underway with the final objective of developing a fundamental understanding of the controlling mechanisms underlying the flammability diagrams of solid combustible materials and their derived fire properties. Given that there is a high possibility of an accidental fire occurring in a space-based facility, understanding the fire properties of materials that will be used in such facilities is of critical importance. With this purpose, the flammability diagrams of the materials, as those produced by the Lateral Ignition and Flame Spread Test (LIFT) apparatus and by a new forced flow device, the Forced Flow Ignition and Flame Spread Test (FIST) apparatus, will be obtained. The specific objective of the program is to apply the new flammability apparatus, which will more accurately reflect the potential ambient conditions of space-based environments, to the characterization of the materials for space applications. This paper presents a parametric study of oxidizer flow effects on the ignition curve of the flammability diagrams of PMMA. The dependence of the ignition delay time on the external radiant flux and either the sample width (LIFT) or the flow velocity (FIST) has been studied. Although preliminary, the results indicate that natural and forced convection flow changes, affect the characteristics of the ignition curves of the flammability diagrams. The major effect on the ignition time appears to be due to convective transfer variations at the fuel surface. At high radiant fluxes or high flow velocities, however, it appears that gas phase processes become increasingly important, affecting the overall ignition delay time. A numerical analysis of the solid fuel heating and pyrolysis has also been developed. The theoretical predictions approximate the experiments well for conditions in which the gas phase induction time is negligible.

Cordova, J. L.↗

Error analysis of finite element solutions for postbuckled cylinders

A general method of error analysis and correction is investigated for the discrete finite-element results for cylindrical shell structures. The method for error analysis is an adaptation of the method of successive approximation. When applied to the equilibrium equations of shell theory, successive approximations derive an approximate continuous solution from the discrete finite-element results. The advantage of this continuous solution is that it contains continuous partial derivatives of an order higher than the basis functions of the finite-element solution. Preliminary numerical results are presented in this paper for the error analysis of finite-element results for a postbuckled stiffened cylindrical panel modeled by a general purpose shell code. Numerical results from the method have previously been reported for postbuckled stiffened plates. A procedure for correcting the continuous approximate solution by Newton's method is outlined.

Sistla, Rajaram↗

Benchmark solutions for the galactic ion transport equations: Energy and spatially dependent problems

Nontrivial benchmark solutions are developed for the galactic ion transport (GIT) equations in the straight-ahead approximation. These equations are used to predict potential radiation hazards in the upper atmosphere and in space. Two levels of difficulty are considered: (1) energy independent, and (2) spatially independent. The analysis emphasizes analytical methods never before applied to the GIT equations. Most of the representations derived have been numerically implemented and compared to more approximate calculations. Accurate ion fluxes are obtained (3 to 5 digits) for nontrivial sources. For monoenergetic beams, both accurate doses and fluxes are found. The benchmarks presented are useful in assessing the accuracy of transport algorithms designed to accommodate more complex radiation protection problems. In addition, these solutions can provide fast and accurate assessments of relatively simple shield configurations.

Ganapol, Barry D.↗

Low Reynolds Number Nozzle Flow Study

An experimental study of low Reynolds number nozzle flow was performed. A brief comparison was made between some of the experimental performance data and performance predicted by a viscous flow code. The performance of 15, 20, and 25 deg conical nozzles, bell nozzles, and trumpet nozzles was evaluated with unheated nitrogen and hydrogen. The numerical analysis was applied to the conical nozzles only, using an existing viscous flow code that was based on a slender-channel approximation. Although the trumpet and 25 deg conical nozzles had slightly better performance at lower Reynolds numbers, it is unclear which nozzle is superior as all fell within the experimental error band. The numerical rssults were found to agree with experimental results for nitrogen and for some of the hydrogen data. Some code modification is recommended to improve confidence in the performance prediction.

Whalen, Margaret V.↗

A Taylor weak-statement algorithm for hyperbolic conservation laws

Finite element analysis, applied to computational fluid dynamics (CFD) problem classes, presents a formal procedure for establishing the ingredients of a discrete approximation numerical solution algorithm. A classical Galerkin weak-statement formulation, formed on a Taylor series extension of the conservation law system, is developed herein that embeds a set of parameters eligible for constraint according to specification of suitable norms. The derived family of Taylor weak statements is shown to contain, as special cases, over one dozen independently derived CFD algorithms published over the past several decades for the high speed flow problem class. A theoretical analysis is completed that facilitates direct qualitative comparisons. Numerical results for definitive linear and nonlinear test problems permit direct quantitative performance comparisons.

Baker, A. J.↗

Method for compression testing of composite materials at high strain rates

A method is presented for testing composite materials in compression at strain rates up to approximately 500 per s. The method uses a thin ring specimen (4 in. in diameter, 1 in. wide, six-eight plies thick) loaded dynamically by an external pressure pulse applied explosively through a liquid. Strains in the specimen and in a steel calibration ring are recoorded with a digital processing oscilloscope. Results are plotted by an x-y plotter in the form of a dynamic stress-strain curve. Data analysis is based on a numerical solution of the equation of motion. A computer program is used which involves smoothing and approximation of the strain magnitude, strain rate, and strain acceleration. Dynamic stress-strain curves obtained for 0-deg and 90-deg specimens of two graphite/epoxy composites are presented.

Daniel, I. M.↗

Improving the five-point bootstrap

We present a new algorithm for the numerical evaluation of five-point conformal blocks in d-dimensions, greatly improving the efficiency of their computation. To do this we use an appropriate ansatz for the blocks as a series expansion in radial coordinates, derive a set of recursion relations for the unknown coefficients in the ansatz, and evaluate the series using a Padé approximant to accelerate its convergence. We then study the 〈σσϵσσ〉 correlator in the 3d critical Ising model by truncating the operator product expansion (OPE) and only including operators with conformal dimension below a cutoff ∆ ⩽ ∆cutoff. We approximate the contributions of the operators above the cutoff by the corresponding contributions in a suitable disconnected five-point correlator. Using this approach, we compute a number of OPE coefficients with greater accuracy than previous methods.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Systems analysis research

Shadow equation for satellite, consistent ephemeris of major planets in solar system, correction process for lunar orbit, and numerical approximation of Hermite expansion

Source record↗