Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical approximations”

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 577 records · Page 32

Essentially nonoscillatory postprocessing filtering methods

High order accurate centered flux approximations used in the computation of numerical solutions to nonlinear partial differential equations produce large oscillations in regions of sharp transitions. Here, we present a new class of filtering methods denoted by Essentially Nonoscillatory Least Squares (ENOLS), which constructs an upgraded filtered solution that is close to the physically correct weak solution of the original evolution equation. Our method relies on the evaluation of a least squares polynomial approximation to oscillatory data using a set of points which is determined via the ENO network. Numerical results are given in one and two space dimensions for both scalar and systems of hyperbolic conservation laws. Computational running time, efficiency, and robustness of method are illustrated in various examples such as Riemann initial data for both Burgers' and Euler's equations of gas dynamics. In all standard cases, the filtered solution appears to converge numerically to the correct solution of the original problem. Some interesting results based on nonstandard central difference schemes, which exactly preserve entropy, and have been recently shown generally not to be weakly convergent to a solution of the conservation law, are also obtained using our filters.

Lafon, F.↗

Numerical studies of identification in nonlinear distributed parameter systems

An abstract approximation framework and convergence theory for the identification of first and second order nonlinear distributed parameter systems developed previously by the authors and reported on in detail elsewhere are summarized and discussed. The theory is based upon results for systems whose dynamics can be described by monotone operators in Hilbert space and an abstract approximation theorem for the resulting nonlinear evolution system. The application of the theory together with numerical evidence demonstrating the feasibility of the general approach are discussed in the context of the identification of a first order quasi-linear parabolic model for one dimensional heat conduction/mass transport and the identification of a nonlinear dissipation mechanism (i.e., damping) in a second order one dimensional wave equation. Computational and implementational considerations, in particular, with regard to supercomputing, are addressed.

Banks, H. T.↗

Application of Newton's method to the postbuckling of rings under pressure loadings

The postbuckling response of circular rings (or long cylinders) is examined. The rings are subjected to four types of external pressure loadings; each type of pressure is defined by its magnitude and direction at points on the buckled ring. Newton's method is applied to the nonlinear differential equations of the exact inextensional theory for the ring problem. A zeroth approximation for the solution of the nonlinear equations, based on the mode shape corresponding to the first buckling pressure, is derived in closed form for each of the four types of pressure. The zeroth approximation is used to start the iteration cycle in Newton's method to compute numerical solutions of the nonlinear equations. The zeroth approximations for the postbuckling pressure-deflection curves are compared with the converged solutions from Newton's method and with similar results reported in the literature.

Thurston, Gaylen A.↗

Discrete observability and numerical quadrature

The authors consider the problem of approximate observability of a one-dimensional diffusion equation on a finite spatial domain with spatial point measurements. The problem of the optimal selection of the measurement points is considered under three conditions: (1) no preassigned measurement nodes; (2) one preassigned node and; (3) two preassigned nodes. The main observation is that the optimal choice is related to three classical procedures in numerical analysis: (1) Gaussian quadrature; (2) Radau quadrature and; (3) Lobatto quadrature. It is shown that the existence of the Radau and Lobatto quadrature is closely related to classical root locus theory.

Martin, Clyde F.↗

A comparison of transport algorithms for premixed, laminar steady state flames

The effects of different methods of approximating multispecies transport phenomena in models of premixed, laminar, steady state flames were studied. Five approximation methods that span a wide range of computational complexity were developed. Identical data for individual species properties were used for each method. Each approximation method is employed in the numerical solution of a set of five H2-02-N2 flames. For each flame the computed species and temperature profiles, as well as the computed flame speeds, are found to be very nearly independent of the approximation method used. This does not indicate that transport phenomena are unimportant, but rather that the selection of the input values for the individual species transport properties is more important than the selection of the method used to approximate the multispecies transport. Based on these results, a sixth approximation method was developed that is computationally efficient and provides results extremely close to the most sophisticated and precise method used.

Coffee, T. P.↗

Connections between conventional and singular-value-based multi-variable flight control system design techniques

It is proposed that frequency-domain multivariable robustness techniques, when combined with classical multivariable procedures, can offer an additional means of evaluating FCS designs. A lateral-directional FCS for an advanced fighter is used as an example. Robustness to unstructured aircraft-input uncertainties is assessed using purely numerical singular-value procedures. Literal approximations for the singular values of the open-loop plant and controller and for the inverse return difference are shown to provide a means of decomposing and diagnosing robustness problems that are insoluble via purely numerical methods.

Mcruer, D. T.↗

OH PLIF Visualization of a Premixed Ethylene-fueled Dual-Mode Scramjet Combustor

Hydroxyl radical (OH) planar induced laser fluorescence (PLIF) measurements have been performed in a small-scale scramjet combustor at the University of Virginia Aerospace Research Laboratory at nominal simulated Mach 5 enthalpy. OH lines were carefully chosen to have fluorescent signal that is independent of pressure and temperature but linear with mole fraction. The OH PLIF signal was imaged in planes orthogonal to and parallel to the freestream flow at different equivalence ratios. Flameout limits were tested and identified. Instantaneous planar images were recorded and analyzed to compare the results with width increased dual-pump enhanced coherent anti-Stokes Raman spectroscopy (WIDECARS) measurements in the same facility and large eddy simulation/Reynolds average Navier-Stokes (LES/RANS) numerical simulation. The flame angle was found to be approximately 10 degrees for several different conditions, which is in agreement with numerical predictions and measurements using WIDECARS. Finally, a comparison between NO PLIF non-combustion cases and OH PLIF combustion cases is provided: the comparison reveals that the dominant effect of flame propagation is freestream turbulence rather than heat release and concentration gradients.

Cantu, Luca M. L.↗

Vibration-translation energy transfer in anharmonic diatomic molecules. 2: The vibrational quantum number dependence

A semiclassical model of the inelastic collision between a vibrationally excited anharmonic oscillator and a structureless atom was used to predict the variation of thermally averaged vibration-translation rate coefficients with temperature and initial-state quantum number. Multiple oscillator states were included in a numerical solution for collinear encounters. The results are compared with CO-He experimental values for both ground and excited initial states using several simplified forms of the interaction potential. The numerical model was also used as a basis for evaluating several less complete but analytic models. Two computationally simple analytic approximations were found that successfully reproduced the numerical rate coefficients for a wide range of molecular properties and collision partners. Their limitations were also identified. The relative rates of multiple-quantum transitions from excited states were evaluated for several molecular types.

Mckenzie, R. L.↗

Vibration-translation energy transfer in anharmonic diatomic molecules. II - The vibrational quantum-number dependence

A semiclassical model of the inelastic collision between a vibrationally excited anharmonic oscillator and a structureless atom is used to predict the variation of thermally averaged vibrational-translational rate coefficients with temperature and initial-state quantum number. Multiple oscillator states are included in a numerical solution for collinear encounters. The results are compared with CO-He experimental values for both ground and excited initial states using several simplified forms of the interaction potential. The numerical model is also used as a basis for evaluating several less complete, but analytic, models. Two computationally simple analytic approximations are found that successfully reproduce the numerical rate coefficients for a wide range of molecular properties and collision partners. Their limitations are identified, and the relative rates of multiple-quantum transitions from excited states are evaluated for several molecular types.

Mckenzie, R. L.↗

Numerical Estimation of the Curvature of Biological Surfaces

Many biological systems may profitably be studied as surface phenomena. A model consisting of isotropic growth of a curved surface from a flat sheet is assumed. With such a model, the Gaussian curvature of the final surface determines whether growth rate of the surface is subharmonic or superharmonic. These properties correspond to notions of convexity and concavity, and thus to local excess growth and local deficiency of growth. In biological models where the major factors controlling surface growth are intrinsic to the surface, researchers thus gained from geometrical study information on the differential growth undergone by the surface. These ideas were applied to an analysis of the folding of the cerebral cortex, a geometrically rather complex surface growth. A numerical surface curvature technique based on an approximation to the Dupin indicatrix of the surface was developed. A metric for comparing curvature estimates is introduced, and considerable numerical testing indicated the reliability of this technique.

Todd, P. H.↗

Data Structure and Parallel Decomposition Considerations on a Fibonacci Grid

The Fibonacci grid, proposed by Swinbank and Purser (see companion abstract), provides attractive properties for global numerical atmospheric prediction by offering an optimally homogeneous, geometrically regular, and approximately isotropic discretization, with only the polar regions requiring special numerical treatment. It is a mathematical idealization, applied to the sphere, of the multi-spiral patterns often found in botanical structures, such as in pine cones and sunflower heads. Computationally, it is natural to organize the domain, into zones, in each of which the same pair, or triple, of "Fibonacci spirals" dominate. But the further subdivision of such zones into "tiles" of a shape and size suitable for distribution to the processors of a massively parallel computer requires very careful consideration if the subsequent spatial computations along the respective spirals, especially those computations (such as compact differencing schemes) that involve recursion, can be implemented in an efficient "load-balanced "manner without requiring excessive amounts of inter-processor communications. In this paper we show how certain "number theoretic" properties of the Fibonacci sequence (whose numbers prescribe the multiplicity of successive spirals) may be exploited in the decomposition of grid zones into tidy arrangements of triangular grid tiles, each tile possessing one side approximately parallel to the constant-latitude zone boundary. We also describe how the spatially recursive processes may be decomposed across such a tiling, and the directionality of the recursions reversed on alternate grid lines, to ensure a very high degree of load balancing throughout the execution of the computations required for one time step of a global model.

Michalakes, John↗

Structural optimization with approximate sensitivities

Computational efficiency in structural optimization can be enhanced if the intensive computations associated with the calculation of the sensitivities, that is, gradients of the behavior constraints, are reduced. Approximation to gradients of the behavior constraints that can be generated with small amount of numerical calculations is proposed. Structural optimization with these approximate sensitivities produced correct optimum solution. Approximate gradients performed well for different nonlinear programming methods, such as the sequence of unconstrained minimization technique, method of feasible directions, sequence of quadratic programming, and sequence of linear programming. Structural optimization with approximate gradients can reduce by one third the CPU time that would otherwise be required to solve the problem with explicit closed-form gradients. The proposed gradient approximation shows potential to reduce intensive computation that has been associated with traditional structural optimization.

Patnaik, S. N.↗

Comparison of Aircraft Models and Integration Schemes for Interval Management in the TRACON

Reusable models of common elements for communication, computation, decision and control in air traffic management are necessary in order to enable simulation, analysis and assurance of emergent properties, such as safety and stability, for a given operational concept. Uncertainties due to faults, such as dropped messages, along with non-linearities and sensor noise are an integral part of these models, and impact emergent system behavior. Flight control algorithms designed using a linearized version of the flight mechanics will exhibit error due to model uncertainty, and may not be stable outside a neighborhood of the given point of linearization. Moreover, the communication mechanism by which the sensed state of an aircraft is fed back to a flight control system (such as an ADS-B message) impacts the overall system behavior; both due to sensor noise as well as dropped messages (vacant samples). Additionally simulation of the flight controller system can exhibit further numerical instability, due to selection of the integration scheme and approximations made in the flight dynamics. We examine the theoretical and numerical stability of a speed controller under the Euler and Runge-Kutta schemes of integration, for the Maintain phase for a Mid-Term (2035-2045) Interval Management (IM) Operational Concept for descent and landing operations. We model uncertainties in communication due to missed ADS-B messages by vacant samples in the integration schemes, and compare the emergent behavior of the system, in terms of stability, via the boundedness of the final system state. Any bound on the errors incurred by these uncertainties will play an essential part in a composable assurance argument required for real-time, flight-deck guidance and control systems,. Thus, we believe that the creation of reusable models, which possess property guarantees, such as safety and stability, is an innovative and essential requirement to assessing the emergent properties of novel airspace concepts of operation.

Neogi, Natasha↗

Reflection and transmission of acoustical waves from a layer with space-dependent velocity.

The refraction of acoustical waves by a moving medium layer is theoretically treated and the reflection and transmission coefficients are determined. The moving-medium-layer velocity is uniform but with a space dependence in one direction. A partitioning of the moving medium layer into constant-velocity sublayers is introduced and numerical results for a three-sublayer approximation of Poiseuille flow are presented. The degenerate case of a single constant-velocity layer is also treated theoretically and numerically. The numerical results show the reflection and transmission coefficients as functions of the peak moving-medium-layer normalized velocity for several angles of incidence.

Steinmetz, G. G.↗

Twisted bilayer graphene. II. Stable symmetry anomaly

We show that the entire continuous model of twisted bilayer graphene (TBG) (and not just the two active bands) with particle-hole symmetry is anomalous and hence incompatible with lattice models. Previous works [e.g., Z. Song, Z. Wang, W. Shi, G. Li, C. Fang, and B. A. Bernevig, Phys. Rev. Lett. 123, 036401 (2019); J. Ahn, S. Park, and B.-J. Yang, Phys. Rev. X 9, 021013 (2019); H. C. Po, L. Zou, T. Senthil, and A. Vishwanath, Phys. Rev. B 99, 195455 (2019); J. Kang and O. Vafek, Phys. Rev. X 8, 031088 (2018); M. Koshino, N. F. Q. Yuan, T. Koretsune, M. Ochi, K. Kuroki, and L. Fu, Phys. Rev. X 8, 031087 (2018); J. Liu, J. Liu, and X. Dai, Phys. Rev. B 99, 155415 (2019); L. Zou, H. C. Po, A. Vishwanath, and T. Senthil, Phys. Rev. B 98, 085435 (2018)] found that the two flatbands in TBG possess a fragile topology protected by the C 2z T symmetry. Song et al. [Z. Song, Z. Wang, W. Shi, G. Li, C. Fang, and B. A. Bernevig, Phys. Rev. Lett. 123, 036401 (2019)] also pointed out an approximate particle-hole symmetry (P) in the continuous model of TBG. In this paper, we numerically confirm that P is indeed a good approximation for TBG and show that the fragile topology of the two flatbands is enhanced to a P-protected stable topology. This stable topology implies 4l + 2 (l ϵ N) Dirac points between the middle two bands. The P-protected stable topology is robust against arbitrary gap closings between the middle two bands and the other bands. We also show that, remarkably, this P-protected stable topology, as well as the corresponding 4l + 2 Dirac points, cannot be realized in lattice models that preserve both C 2z T and P symmetries. In other words, the continuous model of TBG is anomalous and cannot be realized on lattices. Two other topology related topics, with consequences for the interacting TBG problem, i.e., the choice of Chern band basis in the two flatbands and the perfect metal phase of TBG in the so-called second chiral limit, are also discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A rotor optimization using regression analysis

The design and development of helicopter rotors is subject to the many design variables and their interactions that effect rotor operation. Until recently, selection of rotor design variables to achieve specified rotor operational qualities has been a costly, time consuming, repetitive task. For the past several years, Kaman Aerospace Corporation has successfully applied multiple linear regression analysis, coupled with optimization and sensitivity procedures, in the analytical design of rotor systems. It is concluded that approximating equations can be developed rapidly for a multiplicity of objective and constraint functions and optimizations can be performed in a rapid and cost effective manner; the number and/or range of design variables can be increased by expanding the data base and developing approximating functions to reflect the expanded design space; the order of the approximating equations can be expanded easily to improve correlation between analyzer results and the approximating equations; gradients of the approximating equations can be calculated easily and these gradients are smooth functions reducing the risk of numerical problems in the optimization; the use of approximating functions allows the problem to be started easily and rapidly from various initial designs to enhance the probability of finding a global optimum; and the approximating equations are independent of the analysis or optimization codes used.

Giansante, N.↗

The numerical calculation of laminar boundary-layer separation

Iterative finite-difference techniques are developed for integrating the boundary-layer equations, without approximation, through a region of reversed flow. The numerical procedures are used to calculate incompressible laminar separated flows and to investigate the conditions for regular behavior at the point of separation. Regular flows are shown to be characterized by an integrable saddle-type singularity that makes it difficult to obtain numerical solutions which pass continuously into the separated region. The singularity is removed and continuous solutions ensured by specifying the wall shear distribution and computing the pressure gradient as part of the solution. Calculated results are presented for several separated flows and the accuracy of the method is verified. A computer program listing and complete solution case are included.

Klineberg, J. M.↗

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.↗