Development of a three-dimensional meshfree detonation tracking method
Explore the source record for details and available documents.
Mathematics and Computing↗
SEARCH · Engineering Papers
Results for “Mathematical 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.
Quantifying Uncertainties Through Advanced Theoretical Analysis for Image and Signal Reconstruction, NLV-02-15, Year 2 of 2
Presentation describes project to (1) adapt accelerated matrix-free linear system solvers for applying new statistical approaches to very large–scale problems, (2) characterize the degree to which new statistical methods are superior to classical methods in respect to the sensitivity to “tuning parameters,” and (3) develop a new mathematical theory and associated numerical methods for computing the point spread function of an x-ray imaging system.
97 MATHEMATICS AND COMPUTING↗
Vortex methods for separated flows
The numerical solution of the Euler or Navier-Stokes equations by Lagrangian vortex methods is discussed. The mathematical background is presented and includes the relationship with traditional point-vortex studies, convergence to smooth solutions of the Euler equations, and the essential differences between two and three-dimensional cases. The difficulties in extending the method to viscous or compressible flows are explained. Two-dimensional flows around bluff bodies are emphasized. Robustness of the method and the assessment of accuracy, vortex-core profiles, time-marching schemes, numerical dissipation, and efficient programming are treated. Operation counts for unbounded and periodic flows are given, and two algorithms designed to speed up the calculations are described.
Spalart, Philippe R.↗
Hermite analogs of the lowest order Raviart–Thomas mixed method for convection–diffusion equations
The Raviart–Thomas mixed finite-element method of the lowest order (Raviart and Thomas in mixed finite-element methods for second-order elliptic problems, Lecture Notes in mathematics, Springer, New York, 1977) commonly known as the RT{sub 0} method, is a well-established and popular numerical tool to solve diffusion-like problems providing flux continuity across inter-element boundaries. Douglas and Roberts extended the method to the case of more general second-order boundary-value problems including the convection–diffusion equations (cf. this journal Douglas in Comput Appl Math 1:91–103; 1982). The main drawback of these methods, however, is the poor representation of the primal variable by piecewise constant functions. The Hermite analog of the RT{sub 0} method for treating pure diffusion phenomena proposed in Ruas (J Comput Appl Math, 246:234–242; 2013) proved to be a valid alternative to attain higher order approximation of the primal variable while keeping intact the matrix structure and the quality of the discrete flux variable of the original RT{sub 0} method. Non-trivial extensions of this method are studied here that can be viewed as Hermite analogs of the two Douglas and Roberts’ versions of the RT{sub 0} method, to solve convection–diffusion equations. A detailed convergence study is carried out for one of the Hermite methods, and numerical results illustrate the performance of both of them, as compared to each other and to the corresponding mixed methods.
97 MATHEMATICS AND COMPUTING↗
Particle size distributions and the vertical distribution of suspended matter in the upwelling region off Oregon
Various methods of presenting and mathematically describing particle size distribution are explained and evaluated. The hyperbolic distribution is found to be the most practical but the more complex characteristic vector analysis is the most sensitive to changes in the shape of the particle size distributions. A method for determining onshore-offshore flow patterns from the distribution of particulates was presented. A numerical model of the vertical structure of two size classes of particles was developed. The results show a close similarity to the observed distributions but overestimate the particle concentration by forty percent. This was attributed to ignoring grazing by zooplankton. Sensivity analyses showed the size preference was most responsive to the maximum specific growth rates and nutrient half saturation constants. The verical structure was highly dependent on the eddy diffusivity followed closely by the growth terms.
Kitchen, J. C.↗
Improved Estimation Of Vibrations In Engines And Pumps
Method of estimating vibrational loads upon, and vibrational power dissipated in, rocket engines, pumps, turbomachines, and similar flow-through-type machinery derived from older vibration-scaling concept of estimating vibrations of components of rocket engine in design stage, without complicated mathematical modeling and calculations. Methods involve scaling vibrations via global variables like mass-flow rates and exit velocities of engine exhausts. Present method called "Improved Composite Loads Spectra-Component Scaling Criteria".
Chow, David C.↗
The NASA/industry Design Analysis Methods for Vibrations (DAMVIBS) program: Boeing Helicopters airframe finite element modeling
Mathematical models based on the finite element method of structural analysis, as embodied in the NASTRAN computer code, are routinely used by the helicopter industry to calculate airframe static internal loads used for sizing structural members. Historically, less reliance has been placed on the vibration predictions based on these models. Beginning in the early 1980's NASA's Langley Research Center initiated an industry wide program with the objective of engendering the needed trust in vibration predictions using these models and establishing a body of modeling guides which would enable confident future prediction of airframe vibration as part of the regular design process. Emphasis in this paper is placed on the successful modeling of the Army/Boeing CH-47D which showed reasonable correlation with test data. A principal finding indicates that improved dynamic analysis requires greater attention to detail and perhaps a finer mesh, especially the mass distribution, than the usual stress model. Post program modeling efforts show improved correlation placing key modal frequencies in the b/rev range with 4 percent of the test frequencies.
Gabel, R.↗
Implicit Multirate GARK Methods
This article considers multirate generalized-structure additively partitioned Runge–Kutta methods for solving stiff systems of ordinary differential equations with multiple time scales. These methods treat different partitions of the system with different timesteps for a more targeted and efficient solution compared to monolithic single rate approaches. With implicit methods used across all partitions, methods must find a balance between stability and the cost of solving nonlinear equations for the stages. In order to characterize this important trade-off, we explore multirate coupling strategies, problems for assessing linear stability, and techniques to efficiently implement Newton iterations for stage equations. Unlike much of the existing multirate stability analysis which is limited in scope to particular methods, we present general statements on stability and describe fundamental limitations for certain types of multirate schemes. New implicit multirate methods up to fourth order are derived, and their accuracy and efficiency properties are verified with numerical tests.
97 MATHEMATICS AND COMPUTING↗
A formally certified end-to-end implementation of Shor’s factorization algorithm
Quantum computing technology may soon deliver revolutionary improvements in algorithmic performance, but it is useful only if computed answers are correct. While hardware-level decoherence errors have garnered significant attention, a less recognized obstacle to correctness is that of human programming errors—“bugs.” Techniques familiar to most programmers from the classical domain for avoiding, discovering, and diagnosing bugs do not easily transfer, at scale, to the quantum domain because of its unique characteristics. To address this problem, we have been working to adapt formal methods to quantum programming. With such methods, a programmer writes a mathematical specification alongside the program and semiautomatically proves the program correct with respect to it. The proof’s validity is automatically confirmed—certified—by a “proof assistant.” Formal methods have successfully yielded high-assurance classical software artifacts, and the underlying technology has produced certified proofs of major mathematical theorems. As a demonstration of the feasibility of applying formal methods to quantum programming, we present a formally certified end-to-end implementation of Shor’s prime factorization algorithm, developed as part of a framework for applying the certified approach to general applications. By leveraging our framework, one can significantly reduce the effects of human errors and obtain a high-assurance implementation of large-scale quantum applications in a principled way.
Science & Technology - Other Topics↗
Denoising and Extension of Response Functions in the Time Domain
Response functions of quantum systems, such as electron Green’s functions, magnetic, or charge susceptibilities, describe the response of a system to an external perturbation. They are the central objects of interest in field theories and quantum computing and measured directly in experiment. Further, response functions are intrinsically causal. In equilibrium and steady-state systems, they correspond to a positive spectral function in the frequency domain. Since response functions define an inner product on a Hilbert space and thereby induce a positive definite function, the properties of this function can be used to reduce noise in measured data and, in equilibrium and steady state, to construct positive definite extensions for data known on finite time intervals, which are then guaranteed to correspond to positive spectra.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
The fully nonconforming virtual element method for biharmonic problems
Here in this paper, we address the numerical approximation of linear fourth-order elliptic problems on polygonal meshes. In particular, we present a novel nonconforming virtual element discretization of arbitrary order of accuracy for biharmonic problems. The approximation space is made of possibly discontinuous functions, thus giving rise to the fully nonconforming virtual element method. We derive optimal error estimates in a suitable (broken) energy norm and present numerical results to assess the validity of the theoretical estimates.
97 MATHEMATICS AND COMPUTING↗
Cause and Cure - Deterioration in Accuracy of CFD Simulations With Use of High-Aspect-Ratio Triangular Tetrahedral Grids
Traditionally high-aspect ratio triangular/tetrahedral meshes are avoided by CFD re-searchers in the vicinity of a solid wall, as it is known to reduce the accuracy of gradient computations in those regions and also cause numerical instability. Although for certain complex geometries, the use of high-aspect ratio triangular/tetrahedral elements in the vicinity of a solid wall can be replaced by quadrilateral/prismatic elements, ability to use triangular/tetrahedral elements in such regions without any degradation in accuracy can be beneficial from a mesh generation point of view. The benefits also carry over to numerical frameworks such as the space-time conservation element and solution element (CESE), where triangular/tetrahedral elements are the mandatory building blocks. With the requirement of the CESE method in mind, a rigorous mathematical framework that clearly identities the reason behind the difficulties in use of such high-aspect ratio triangular/tetrahedral elements is presented here. As will be shown, it turns out that the degree of accuracy deterioration of gradient computation involving a triangular element is hinged on the value of its shape factor Gamma def = sq sin Alpha1 + sq sin Alpha2 + sq sin Alpha3, where Alpha1; Alpha2 and Alpha3 are the internal angles of the element. In fact, it is shown that the degree of accuracy deterioration increases monotonically as the value of Gamma decreases monotonically from its maximal value 9/4 (attained by an equilateral triangle only) to a value much less than 1 (associated with a highly obtuse triangle). By taking advantage of the fact that a high-aspect ratio triangle is not necessarily highly obtuse, and in fact it can have a shape factor whose value is close to the maximal value 9/4, a potential solution to avoid accuracy deterioration of gradient computation associated with a high-aspect ratio triangular grid is given. Also a brief discussion on the extension of the current mathematical framework to the tetrahedral-grid case along with some of the practical results of this extension is also provided. Furthermore, through the use of numerical simulations of practical viscous problems involving high-Reynolds number flows, the effectiveness of the gradient evaluation procedures within the CESE framework (that have their basis on the analysis presented here) to produce accurate and stable results on such high-aspect ratio meshes is also showcased.
CESE Method↗
Discontinuous Skeletal Gradient Discretisation methods on polytopal meshes
Here, in this work we develop arbitrary-order Discontinuous Skeletal Gradient Discretisations (DSGD) on general polytopal meshes. Discontinuous Skeletal refers to the fact that the globally coupled unknowns are broken polynomials on the mesh skeleton. The key ingredient is a high-order gradient reconstruction composed of two terms: (i) a consistent contribution obtained mimicking an integration by parts formula inside each element and (ii) a stabilising term for which sufficient design conditions are provided. An example of stabilisation that satisfies the design conditions is proposed based on a local lifting of high-order residuals on a Raviart–Thomas–Nédélec subspace. We prove that the novel DSGDs satisfy coercivity, consistency, limit-conformity, and compactness requirements that ensure convergence for a variety of elliptic and parabolic problems. Lastly, links with Hybrid High-Order, non-conforming Mimetic Finite Difference and non-conforming Virtual Element methods are also studied. Numerical examples complete the exposition.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Designing Shading Schemes for Microphone Phased Arrays
A technique is proposed for designing sensor weightings, or array shading functions, for microphone phased arrays. The presented method attempts to maintain mathematical simplicity while maximizing array gain for a desired mainlobe beamwidth. It does so using the two-parameter optimization of a family of functions. The shading design process is described and applied to an array used recently in the NASA Langley Research Center 14- by 22-Foot Subsonic Tunnel. Several variants of the method are compared to an existing shading scheme for the array, along with a more general optimization formulation. Results are presented for both an idealized point source and experimental data from the wind tunnel test. The shading schemes are shown to improve the behavior of array operations, but have a stronger impact on visual interpretation of maps than on source level calculations. The study results show that the two-parameter optimization scheme provides favorable characteristics while maintaining implementation simplicity.
Christopher J Bahr↗
The Gaussian CL s method for searches of new physics
Here we describe a method based on the CL s approach to present results in searches of new physics, under the condition that the relevant parameter space is continuous. Our method relies on a class of test statistics developed for non-nested hypotheses testing problems, denoted by ΔT, which has a Gaussian approximation to its parent distribution when the sample size is large. This leads to a simple procedure of forming exclusion sets for the parameters of interest, which we call the Gaussian CL s method. Our work provides a self-contained mathematical proof for the Gaussian CL s method, that explicitly outlines the required conditions. These conditions are milder than that required by the Wilks' theorem to set confidence intervals (CIs). We illustrate the Gaussian CL s method in an example of searching for a sterile neutrino, where the CL s approach was rarely used before. We also compare data analysis results produced by the Gaussian CL s method and various CI methods to showcase their differences.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Frequency-Domain Modeling Of Dynamics Of Helicopters
Method for identification, in frequency domain, of parameters of mathematical models of dynamical behaviors of helicopters combines and extends several existing methods. Focuses on linear-state-space models. Comprehensive frequency-domain method for identification of parameters of mathematical model helicopter involves synthesis of several prior methods.
Tischler, Mark B.↗
Energy Saving by Insulating a Cooking Pot
Sandia National Laboratories is one of the most technologically advanced and resourceful national laboratories in the United States. Their mission is to use science and technology, innovative research, and global engagement to counter threats, reduce dangers, and respond to disasters. To achieve this mission, Sandia is divided into ten divisions that are overlooked by Dr. Stephen Younger, the Laboratories Director. Amongst the ten divisions, the internship will take place in Division 8000, which is Integrated Security Solutions. Inside this division, the internship will be a Structural and Thermal Analysis R&D. Apart from being strictly focused on security, the division performs research for small companies that do not count with the resources to perform it themselves. During the summer internship, one of these cases is being attended. A small company invented an insulating cover that is supposed to be placed over the top of a cooking pot lid. They claim that the insulating cover helps in the saving of energy. As an intern and researcher, the task is to determine how much energy is effectively saved by placing the insulating cover over the cooking pot. To achieve the goal, mathematical analyses, experiments and simulations will be done. Then, if appropriate, a better design will be made to increase the energy saved by the insulating cover. For the mathematical analysis, the thermal resistance method was used, resulting in energy saved at a rate of 14.2795 W. After the mathematical analysis was done, a test plan was developed and presented to the supervisor for approval. After some recommendations, the test plan will be modified accordingly.
32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
Computing Misalignments From Measured Contact Resistances
Misalignments of contacts in integrated circuits computed, and measured contact resistances corrected for effects of misalignments by use of method based on relatively simple mathematical model of measured contact resistance. Flange resistance of aligned contact multiplied by semiempirical factor. Method particularly useful in estimating effects of misalignments on four- and six-terminal (Kelvin-tap) contact-resistance measurements taken with L-shaped test structures.
Lieneweg, Udo↗