SEARCH · Engineering Papers

Results for “numerical computation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 289 records · Page 16

All sources NASA NTRS DOE OSTI

NASA NTRS1989

Frequency response modeling and control of flexible structures: Computational methods

The dynamics of vibrations in flexible structures can be conventiently modeled in terms of frequency response models. For structural control such models capture the distributed parameter dynamics of the elastic structural response as an irrational transfer function. For most flexible structures arising in aerospace applications the irrational transfer functions which arise are of a special class of pseudo-meromorphic functions which have only a finite number of right half place poles. Computational algorithms are demonstrated for design of multiloop control laws for such models based on optimal Wiener-Hopf control of the frequency responses. The algorithms employ a sampled-data representation of irrational transfer functions which is particularly attractive for numerical computation. One key algorithm for the solution of the optimal control problem is the spectral factorization of an irrational transfer function. The basis for the spectral factorization algorithm is highlighted together with associated computational issues arising in optimal regulator design. Options for implementation of wide band vibration control for flexible structures based on the sampled-data frequency response models is also highlighted. A simple flexible structure control example is considered to demonstrate the combined frequency response modeling and control algorithms.

Bennett, William H.↗

NASA NTRS1999

Portable Parallel Programming for the Dynamic Load Balancing of Unstructured Grid Applications

The ability to dynamically adapt an unstructured -rid (or mesh) is a powerful tool for solving computational problems with evolving physical features; however, an efficient parallel implementation is rather difficult, particularly from the view point of portability on various multiprocessor platforms We address this problem by developing PLUM, tin automatic anti architecture-independent framework for adaptive numerical computations in a message-passing environment. Portability is demonstrated by comparing performance on an SP2, an Origin2000, and a T3E, without any code modifications. We also present a general-purpose load balancer that utilizes symmetric broadcast networks (SBN) as the underlying communication pattern, with a goal to providing a global view of system loads across processors. Experiments on, an SP2 and an Origin2000 demonstrate the portability of our approach which achieves superb load balance at the cost of minimal extra overhead.

Biswas, Rupak↗

DOE OSTI2021

JSPEC – An Open Source Program for IBS and Electron Cooling Simulation

The intrabeam scattering can affect the accumulation, the lifetime, and the property of a high-intensity beam. Electron cooling is a method to mitigate the intrabeam scattering effect. JSPEC (JLab Simulation Package on Electron Cooling) is an open-source program developed at Jefferson Lab, which includes various numerical models and friction force formulas for intrabeam scattering and electron cooling simulations. JSPEC has been benchmarked with BETACOOL and experimental data. In this report, we will introduce what features JSPEC provides to the users and how it carries out the computations. Numerical examples are presented to demonstrate the performance and the validity of JSPEC.

Zhang, He↗

DOE OSTI2023

Computation of high-frequency magnetoelastic waves in layered materials

Here, the direct calculation of magnetoelastic wave dispersion in layered media is presented using an efficient, accurate computational technique. The governing, coupled equations for elasticity and magnetism, the Navier and Landau-Lifshitz equations, respectively, are linearized to form a quadratic eigenvalue problem that determines a complex web of wave-number–frequency dispersion branches and their corresponding mode profiles. Numerical discretization of the eigenvalue problem via a spectral collocation method (SCM) is employed to determine the complete dispersion maps for both a single, finite-thickness magnetic layer and a finite magnetic-nonmagnetic double-layer. The SCM, previously used to study elastic waves in nonmagnetic media, is fast, accurate, and adaptable to a variety of sample configurations and geometries. Emphasis is placed on the extremely high-frequency regimes being accessed in ultrafast magnetism experiments. The dispersion maps and modes provide insight into how energy propagates through the coupled system, including how energy can be transferred between elastic- and magnetic-dominated waves as well as between different layers. The numerical computations for a single layer are further understood by a simplified analytical calculation in the high-frequency, exchange-dominated regime where the resonance condition required for energy exchange (an anticrossing) between quasi-elastic and quasi-magnetic dispersion branches is determined. Nonresonant interactions are shown to be well approximated by the dispersion of uncoupled elastic and magnetic waves. The methods and results provide fundamental theoretical tools to model and understand current and future magnetic devices powering spintronic innovation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS2019

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↗

NASA NTRS1980

Experimental unsteady aerodynamics of conventional and supercritical airfoils

Experimental data on the unsteady aerodynamics of oscillating airfoils in transonic flow are presented. Two 0.5 m-chord airfoil models - an NACA 64A010 and an NLR 7301 - were tested in the NASA-Ames 11 by 11 foot Transonic Wind Tunnel at Mach numbers to 0.85, at chord Reynolds numbers to 12 million and at mean angles of attack to 4 deg. The airfoils were subjected to both pitching and plunging motions at reduced frequencies to 0.3 (physical frequencies to 53 Hz). The new hardware and the extensive use of computer-experiment integration developed for this test are described. The geometrical configuration of the model and the test arrangement are described in detail. Mean and first harmonic data are presented in both tabular and graphical form to aid in comparisons with other data and with numerical computations.

Davis, S. S.↗

DOE OSTI2021

Computational approaches to Coherent Synchrotron Radiation in two and three dimensions

Coherent Synchrotron Radiation (CSR) is an important and often detrimental effect in particle accelerators. While one-dimensional models have been successfully used to design and explain the behavior of modern machines, questions remain about their domain of validity. In recent years, two- and three-dimensional models have been developed that are amenable to efficient numerical computation. This article gives an overview of CSR computation from its discovery through the present state of the art.

43 PARTICLE ACCELERATORS↗

NASA NTRS1973

The Computerized Anatomical Man (CAM) model

A computerized anatomical man (CAM) model, representing the most detailed and anatomically correct geometrical model of the human body yet prepared, has been developed for use in analyzing radiation dose distribution in man. This model of a 50-percentile standing USAF man comprises some 1100 unique geometric surfaces and some 2450 solid regions. Internal body geometry such as organs, voids, bones, and bone marrow are explicitly modeled. A computer program called CAMERA has also been developed for performing analyses with the model. Such analyses include tracing rays through the CAM geometry, placing results on magnetic tape in various forms, collapsing areal density data from ray tracing information to areal density distributions, preparing cross section views, etc. Numerous computer drawn cross sections through the CAM model are presented.

Billings, M. P.↗

NASA NTRS1993

Improvement of method for experimental determination of flutter speed by parameter identification

The method investigated in the current research work is based on the identification of the equations of motion during test flights, followed by the solution of these identified equations to numerically compute the flutter dynamic pressure. The current research work is aimed at overcoming two main difficulties which arise when: (1) A continuous system is truncated into a system with a finite number of degrees of freedom. (2) The desired range of frequencies is wide; in the latter case, numerical difficulties occur that can only partly be overcome by giving relatively more weight to the high frequency modes. Results pertaining to the above first point will be presented under section entitled 'Truncation Effects', and results pertaining to the above second point will be presented under section entitled 'Frequency Range Effects'. The numerical example used herein consists of a continuous simply supported beam. It is further assumed that the object of the identification procedure is to identify the beam's first seven modes.

Nissim, E.↗

NASA NTRS1995

A well-posed numerical method to track isolated conformal map singularities in Hele-Shaw flow

We present a new numerical method for calculating an evolving 2D Hele-Shaw interface when surface tension effects are neglected. In the case where the flow is directed from the less viscous fluid into the more viscous fluid, the motion of the interface is ill-posed; small deviations in the initial condition will produce significant changes in the ensuing motion. This situation is disastrous for numerical computation, as small round-off errors can quickly lead to large inaccuracies in the computed solution. Our method of computation is most easily formulated using a conformal map from the fluid domain into a unit disk. The method relies on analytically continuing the initial data and equations of motion into the region exterior to the disk, where the evolution problem becomes well-posed. The equations are then numerically solved in the extended domain. The presence of singularities in the conformal map outside of the disk introduces specific structures along the fluid interface. Our method can explicitly track the location of isolated pole and branch point singularities, allowing us to draw connections between the development of interfacial patterns and the motion of singularities as they approach the unit disk. In particular, we are able to relate physical features such as finger shape, side-branch formation, and competition between fingers to the nature and location of the singularities. The usefulness of this method in studying the formation of topological singularities (self-intersections of the interface) is also pointed out.

Baker, Gregory↗

DOE OSTI2024

A simple solution to the Rietveld refinement recipe problem

Rietveld refinements are widely used for many purposes in the physical sciences. Conducting a Rietveld refinement typically requires expert input because correct results may require that parameters be added to the fit in the proper order. This order will depend on the nature of the data and the initial parameter values. A mechanism for computing the next parameter to add to the refinement is shown. The fitting function is evaluated with the current parameter value set and each parameter incremented and decremented by a small offset. This provides the partial derivatives with respect to each parameter, along with information to discriminate meaningful values from numerical computational errors. The implementation of this mechanism in the open-source GSAS-II program is discussed. This new method is discussed as an important step towards the development of automated Rietveld refinement technology.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1994

A high-speed linear algebra library with automatic parallelism

Parallel or distributed processing is key to getting highest performance workstations. However, designing and implementing efficient parallel algorithms is difficult and error-prone. It is even more difficult to write code that is both portable to and efficient on many different computers. Finally, it is harder still to satisfy the above requirements and include the reliability and ease of use required of commercial software intended for use in a production environment. As a result, the application of parallel processing technology to commercial software has been extremely small even though there are numerous computationally demanding programs that would significantly benefit from application of parallel processing. This paper describes DSSLIB, which is a library of subroutines that perform many of the time-consuming computations in engineering and scientific software. DSSLIB combines the high efficiency and speed of parallel computation with a serial programming model that eliminates many undesirable side-effects of typical parallel code. The result is a simple way to incorporate the power of parallel processing into commercial software without compromising maintainability, reliability, or ease of use. This gives significant advantages over less powerful non-parallel entries in the market.

Boucher, Michael L.↗

NASA NTRS1983

Pair production in superstrong magnetic fields

The production of electron-positron pairs by single photons in magnetic fields 10 to the twelth power G was investigated in detail for photon energies near threshold as well as for the asymptotic limit of high photon energy. The exact attenuation coefficient, which is derived and then evaluated numerically, is strongly influenced by the discrete energy states of the electron and positron. Near threshold, it exhibits a sawtooth pattern as a function of photon energy, and its value is significantly below that predicted by the asymptotic expression for the attenuation coefficient. The energy distributions of the created pair are computed numerically near threshold and analytic expressions are derived in the asymptotic limit. These results indicate that as field strength and photon energy increase, it becomes increasingly probable for the pair to divide the photon energy unequally. This effect, as well as the threshold behavior of the attenuation coefficient, could have important consequences for pulsar models.

Daugherty, J. K.↗

NASA NTRS1983

Pair production in superstrong magnetic fields

The production of electron-positron pairs by single photons in magnetic fields of not less than about 10 to the 12th gauss has been investigated in detail for photon energy. The exact attenuation coefficient, which is derived and then evaluated numerically, is strongly influenced by the discrete energy states of the electron and positron. Near threshold, it exhibits a 'sawtooth' pattern as a function of photon energy and its value is significantly below that predicted by the asymptotic expression for the attenuation coefficient. The energy distributions of the created pair are computed numerically near threshold, and analytic expressions are derived in the asymptotic limit. These results indicate that as field strength and photon energy increase, it becomes increasingly probable for the pair to divide strength energy unequally. This effect, as well as the threshold behavior of the attenuation coefficient, could have important consequences to pulsar models.

Daugherty, J. K.↗

NASA NTRS2023

Implementing a VLBI Time Delay Model for Earth-orbiting Satellites: Partial derivatives and Verification

This document describes the partial derivatives for commonly estimated parameters including antenna positions, satellite position, and satellite velocity from a near-field VLBI delay model for Earth-orbiting satellites. This model was presented in the Journal of Geodesy by Jaron & Nothnagel (2019). In this context, we present a streamlined version of the near-field VLBI delay. From this simplified model, we deduce a delay rate expression and calculate partial derivatives, maintaining only the terms with significant impact on the computed derivatives’ magnitude. To verify the accuracy of the simplified model and the partial derivatives computed from it, we have created a simple simulation in Matlab of an Earth-orbiting satellite at the altitude of a typical Global Navigation Satellite Systems (GNSS) satellite. From this simulation, we compare the simplified and original delays, and we verify the magnitude and direction of the partial derivatives against the numerically computed derivatives from the original delay model. The partial derivatives and simplified VLBI delay model detailed here are implemented in Fortran in the open-source library VTD.

Geodesy↗

NASA NTRS1987

An efficient method to compute spurious end point contributions in PO solutions

A method is given to compute the spurious endpoint contributions in the physical optics solution for electromagnetic scattering from conducting bodies. The method is applicable to general three-dimensional structures. The only information required to use the method is the radius of curvature of the body at the shadow boundary. Thus, the method is very efficient for numerical computations. As an illustration, the method is applied to several bodies of revolution to compute the endpoint contributions for backscattering in the case of axial incidence. It is shown that in high-frequency situations, the endpoint contributions obtained using the method are equal to the true endpoint contributions.

Gupta, Inder J.↗

NASA NTRS1991

The unsteady laminar boundary layer on an axisymmetric body subject to small-amplitude fluctuations in the free-stream velocity

The effect of small-amplitude, time-periodic, free-stream disturbances on an otherwise steady axisymmetric boundary layer on a circular cylinder is considered. Numerical solutions to the problem are presented, and an asymptotic solution to the flow, valid far downstream along the axis of the cylinder is detailed. Particular emphasis is placed on the unsteady eigensolutions that occur far downstream, which turn out to be very different from the analogous planar eigensolutions. These axisymmetric eigensolutions are computed numerically and also are described by asymptotic analyses valid for low and high frequencies of oscillation.

Duck, Peter W.↗

NASA NTRS1990

Research in progress in applied mathematics, numerical analysis, and computer science

Research conducted at the Institute in Science and Engineering in applied mathematics, numerical analysis, and computer science is summarized. The Institute conducts unclassified basic research in applied mathematics in order to extend and improve problem solving capabilities in science and engineering, particularly in aeronautics and space.

Source record↗