Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “method of multipliers”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS2011

In-Service Monitoring of Steam Pipe Systems at High Temperatures

An effective, in-service health monitoring system is needed to track water condensation in real time through the walls of steam pipes. The system is required to measure the height of the condensed water from outside the pipe, while operating at temperatures that are as high as 250 C. The system needs to account for the effects of water flow and cavitation. In addition, it is desired that the system does not require perforating the pipes and thereby reducing the structural integrity. Generally, steam pipes are used as part of the district heating system carrying steam from central power stations under the streets to heat, cool, or supply power to high-rise buildings and businesses. This system uses ultrasonic waves in pulse-echo and acquires reflected signal data. Via autocorrelation, it determines the water height while eliminating the effect of noise and multiple reflections from the wall of the pipe. The system performs nondestructive monitoring through the walls of steam pipes, and automatically measures the height of condensed water while operating at the high-temperature conditions of 250 C. For this purpose, the ultrasonic pulse-echo method is used where the time-of-flight of the wave reflections inside the water are measured, and it is multiplied by the wave velocity to determine the height. The pulse-echo test consists of emitting ultrasonic wave pulses from a piezoelectric transducer and receiving the reflections from the top and bottom of the condensed water. A single transducer is used as a transmitter as well as the receiver of the ultrasonic waves. To obtain high resolution, a broadband transducer is used and the frequency can be in the range of 2.25 to 10 MHz, providing sharp pulses in the time domain allowing for higher resolution in identifying the individual reflections.

Bar-Cohen, Yoseph↗
DOE OSTI2021

Derivation of the Future Time Equation for Analog, Non-Multiplying Monte Carlo Simulation

The expected computational time required to simulate a particle from a point in phase space through a Monte Carlo history, termed the expected future time, is found by solving the Future Time Equation (FTE). The expected future time may be useful when generating variance reduction parameters for a Monte Carlo simulation with a method such as Consistent Adjoint Driven Importance Sampling (CADIS). This report presents a detailed derivation of the Future Time Probability Density Function (FTPDF) and FTE for neutral particle Monte Carlo transport to aid future researchers. For simplicity, this derivation only considers analog transport in non-multiplying media.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1972

The radiative transfer equation and environmental effects in the upper atmosphere.

A highly effective approach to the solution of the radiative transfer equation, including the effects of emission, absorption, and multiply scattered radiation, is described. The transfer equation is cast in its unsteady hyperbolic form which is amenable to several different numerical methods of solution. Comparisons of solutions with emission, absorption, and scattering of radiation, are presented to demonstrate the validity of the present approach.

Callis, L. B.↗
NASA NTRS1994

Three-dimensional hydrodynamic Bondi-Hoyle accretion. 1: Code validation and stationary accretors

We investigate the hydrodynamics of three-dimensional classical Bondi-Hoyle accretion. Totally absorbing stationary spheres of varying sizes (from 10.0 down to 0.02 Bondi radii) accrete matter from a homogeneous and slightly perturbed medium, which is taken to be an ideal gas (gamma = 5/3 or 1.2). To accommodate the long-range gravitational forces, the extent of the computational volume is typically a factor of 100 larger than the radius of the accretor. We compare the numerical mass accretion rates with the theoretical predictions of Bondi, to assess the validity of the code. The hydrodynamics is modeled by the piecewise parabolic method. No energy sources (nuclear burning) or sinks (radiation, conduction) are included. The resolution in the vicinity of the accretor is increased by multiply nesting several (6-8) grids around the stationary sphere, each finer grid being a factor of 2 smaller spatially than the next coarser grid. This allows us to include a coarse model for the surface of the accretor (vacuum sphere) on the finest grid while at the same time evolving the gas on the coarser grids. The accretion rates derived numerically are in in very good agreement (to about 10% over several orders of magnitude) with the values given by Bondi for a stationary accretor within a hydrodynamic medium. However, the equations have to be changed in order to include the finite size of the accretor (in some cases very large compared to the sonic point or even to the Bondi radius).

Ruffert, Maximilian↗
NASA NTRS1994

Three-dimensional hydrodynamic Bondi-Hoyle accretion. 2: Homogeneous medium at Mach 3 with gamma = 5/3

We investigate the hydrodynamics of three-dimensional classical Bondi-Hoyle accretion. Totally absorbing spheres of varying sizes (from 10 down to 0.01 accretion radii) move at Mach 3 relative to a homogeneous and slightly perturbed medium, which is taken to be an ideal gas (gamma = 5/3). To accommodate the long-range gravitational forces, the extent of the computational volume is 32(exp 3) accretion radii. We examine the influence of numerical procedure on physical behavior. The hydrodynamics is modeled by the 'piecewise parabolic method.' No energy sources (nuclear burning) or sinks (radiation, conduction) are included. The resolution in the vicinity of the accretor is increased by multiply nesting several (5-10) grids around the sphere, each finer grid being a factor of 2 smaller in zone dimension that the next coarser grid. The largest dynamic range (ratio of size of the largest grid to size of the finest zone) is 16,384. This allows us to include a coarse model for the surface of the accretor (vacuum sphere) on the finest grid, while at the same time evolving the gas on the coarser grids. Initially (at time t = 0-10), a shock front is set up, a Mach cone develops, and the accretion column is observable. Eventually the flow becomes unstable, destroying axisymmetry. This happens approximately when the mass accretion rate reaches the values (+/- 10%) predicted by the Bondi-Hoyle accretion formula (factor of 2 included). However, our three-dimensional models do not show the highly dynamic flip-flop flow so prominent in two-dimensional calculations performed by other authors. The flow, and thus the accretion rate of all quantities, shows quasi-periodic (P approximately equals 5) cycles between quiescent and active states. The interpolation formula proposed in an accompanying paper is found to follow the collected numerical data to within approximately 30%. The specific angular momentum accreted is of the same order of magnitude as the values previously found for two-dimensional flows.

Ruffert, Maximilian↗
NASA NTRS1970

COPTRAN - A method of optimum communication systems design

Single set of mathematical expressions describes system cost and probability of error of data transmission in terms of four basic parameters in the link equation. A Lagrange multiplier sets up equations whose solutions yield the optimum values for system design considerations and weight and cost values.

Brinkman, K. L.↗
NASA NTRS1988

Nonlinear programming extensions to rational function approximation methods for unsteady aerodynamic forces

The approximation of unsteady generalized aerodynamic forces in the equations of motion of a flexible aircraft are discussed. Two methods of formulating these approximations are extended to include the same flexibility in constraining the approximations and the same methodology in optimizing nonlinear parameters as another currently used extended least-squares method. Optimal selection of nonlinear parameters is made in each of the three methods by use of the same nonlinear, nongradient optimizer. The objective of the nonlinear optimization is to obtain rational approximations to the unsteady aerodynamics whose state-space realization is lower order than that required when no optimization of the nonlinear terms is performed. The free linear parameters are determined using the least-squares matrix techniques of a Lagrange multiplier formulation of an objective function which incorporates selected linear equality constraints. State-space mathematical models resulting from different approaches are described and results are presented that show comparative evaluations from application of each of the extended methods to a numerical example.

Tiffany, Sherwood H.↗
DOE OSTI2023

Discovering strongly lensed quasar candidates with catalogue-based methods from DESI Legacy Surveys

The Hubble tension, revealed by a ~5σ discrepancy between measurements of the Hubble-Lemaitre constant among observations of the early and local Universe, is one of the most significant problems in modern cosmology. In order to better understand the origin of this mismatch, independent techniques to measure H 0 , such as strong lensing time delays, are required. Notably, the sample size of such systems is key to minimising the statistical uncertainties and cosmic variance, which can be improved by exploring the datasets of large-scale sky surveys such as Dark Energy Spectroscopic Instrument (DESI). We identify possible strong lensing time-delay systems within DESI by selecting candidate multiply imaged lensed quasars from a catalogue of 24 440 816 candidate QSOs contained in the ninth data release of the DESI Legacy Imaging Surveys (DESI-LS). Using a friend-of-friends-like algorithm on spatial co-ordinates, our method generates an initial list of compact quasar groups. This list is subsequently filtered using a measure of the similarity of colours among a group’s members and the likelihood that they are quasars. A visual inspection finally selects candidate strong lensing systems based on the spatial configuration of the group members. We identified 620 new candidate multiply imaged lensed quasars (101 grade-A, 214 grade-B, 305 grade-C). This number excludes 53 known spectroscopically confirmed systems and existing candidate systems identified in other similar catalogues. When available, these new candidates will be further checked by combining the spectroscopic and photometric data from DESI.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

A priori error analysis of high-order LL* (FOSLL*) finite element methods

A number of non-standard finite element methods have been proposed in recent years, each of which derives from a specific class of PDE-constrained norm minimization problems. The most notable examples are LL* methods. In this work, we argue that all high-order methods in this class should be expected to deliver substandard uniform h-refinement convergence rates. In fact, one may not even see rates proportional to the polynomial order p > 1 when the exact solution is a constant function. Here, we show that the convergence rate is limited by the regularity of an extraneous Lagrange multiplier variable which naturally appears via a saddle-point analysis. In turn, limited convergence rates appear because the regularity of this Lagrange multiplier is determined, in part, by the geometry of the domain. Numerical experiments support our conclusions.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1951

Analysis of plane-plastic stress problems with axial symmetry in strain-hardening range

A simple method is developed for solving plane-plastic-stress problems with axial symmetry in the strain-hardening range which is based on the deformation theory of plasticity employing the finite-strain concept. The equations defining the problems are first reduced to two simultaneous nonlinear differential equations involving two dependent variables: (a) the octahedral shear strain, and (b) a parameter indicating the ratio of principal stresses. By multiplying the load and dividing the radius by an arbitrary constant, it is possible to solve these problems without iteration for any value of the modified load. The constant is determined by the boundary condition. This method is applied to a circular membrane under pressure, a rotating disk without and with a central hole, and an infinite plate with a circular hole. Two materials, inconel x and 16-25-6, the octahedral shear stress-strain relations of which do not follow the power law, are used. Distributions of octahedral shear strain, as well as of principal stresses and strains, are obtained. These results are compared with the results of the same problems in the elastic range.

Wu, M H Lee↗
NASA NTRS1987

Nonlinear programming extensions to rational function approximations of unsteady aerodynamics

This paper deals with approximating unsteady generalized aerodynamic forces in the equations of motion of a flexible aircraft. Two methods of formulating these approximations are extended to include both the same flexibility in constraining them and the same methodology in optimizing nonlinear parameters as another currently used 'extended least-squares' method. Optimal selection of 'nonlinear' parameters is made in each of the three methods by use of the same nonlinear (nongradient) optimizer. The objective of the nonlinear optimization is to obtain rational approximations to the unsteady aerodynamics whose state-space realization is of lower order than that required when no optimization of the nonlinear terms is performed. The free 'linear' parameters are determined using least-squares matrix techniques on a Lagrange multiplier formulation of an objective function which incorporates selected linear equality constraints. State-space mathematical models resulting from the different approaches are described, and results are presented which show comparative evaluations from application of each of the extended methods to a numerical example. The results obtained for the example problem show a significant (up to 63 percent) reduction in the number of differential equations used to represent the unsteady aerodynamic forces in linear time-invariant equations of motion as compared to a conventional method in which nonlinear terms are not optimized.

Tiffany, Sherwood H.↗
NASA NTRS1992

Estimating new production in the equatorial Pacific Ocean at 150 deg W

A major goal of the WEC88 cruise of the R/V Wecoma to the equatorial Pacific (made in February-March 1988) was to establish rates of new production along a meridional section at 150 deg W and to compare these measured rates with the relatively high values for the equatorial Pacific that had been reported previously using indirect methods and models. Production values were obtained from the traditional approach using N-15 labeled nitrate uptake, and by using C-14 fixation values multiplied by f (proportion of new production) from various sources: from N-15 data, from a C-14 fixation-versus-f relationship, or from a nitrate-versus-f relationship. The ratios of directly measured nitrate and carbon uptake and the ratios of nitrate to nitrate plus ammonium uptake, i.e., values of f, agree well; values of f calculated from carbon uptake or from nitrate concentration are overestimates for the equatorial upwelling region. Carbon-to-nitrogen uptake ratios measured with C-14 and N-15, respectively, approximate the Redfield molar ratio, 6.6 C:N. The overall mean value of f (0.17) helps confirm the view that the low primary production in the enriched eastern equatorial Pacific is due to failure of the nitrate-uptake system.

Dugdale, Richard C.↗
NASA NTRS1992

Development of a fractional-step method for the unsteady incompressible Navier-Stokes equations in generalized coordinate systems

A fractional step method is developed for solving the time-dependent three-dimensional incompressible Navier-Stokes equations in generalized coordinate systems. The primitive variable formulation uses the pressure, defined at the center of the computational cell, and the volume fluxes across the faces of the cells as the dependent variables, instead of the Cartesian components of the velocity. This choice is equivalent to using the contravariant velocity components in a staggered grid multiplied by the volume of the computational cell. The governing equations are discretized by finite volumes using a staggered mesh system. The solution of the continuity equation is decoupled from the momentum equations by a fractional step method which enforces mass conservation by solving a Poisson equation. This procedure, combined with the consistent approximations of the geometric quantities, is done to satisfy the discretized mass conservation equation to machine accuracy, as well as to gain the favorable convergence properties of the Poisson solver. The momentum equations are solved by an approximate factorization method, and a novel ZEBRA scheme with four-color ordering is devised for the efficient solution of the Poisson equation. Several two- and three-dimensional laminar test cases are computed and compared with other numerical and experimental results to validate the solution method. Good agreement is obtained in all cases.

Rosenfeld, Moshe↗
NASA NTRS1987

Transient response for interaction of two dynamic bodies

During the launch sequence of any space vehicle complicated boundary interactions occur between the vehicle and the launch stand. At the start of the sequence large forces exist between the two; contact is then broken in a short but finite time which depends on the release mechanism. The resulting vehicle response produces loads which are very high and often form the design case. It is known that the treatment of the launch pad as a second dynamic body is significant for an accurate prediction of launch response. A technique was developed for obtaining loads generated by the launch transient with the effect of pad dynamics included. The method solves uncoupled vehicle and pad equations of motion. The use of uncoupled models allows the simulation of vehicle launch in a single computer run. Modal formulation allows a closed-form solution to be written, eliminating any need for a numerical integration algorithm. When the vehicle is on the pad the uncoupled pad and vehicle equations have to be modified to account for the constraints they impose on each other. This necessitates the use of an iterative procedure to converge to a solution, using Lagrange multipliers to apply the required constraints. As the vehicle lifts off the pad the coupling between the vehicle and the pad is eliminated point by point until the vehicle flies free. Results obtained by this method were shown to be in good agreement with observed loads and other analysis methods. The resulting computer program is general, and was used without modification to solve a variety of contact problems.

Prabhakar, A.↗
NASA NTRS1985

Parameter identification in distributed systems

This paper describes a method for the identification of the parameters entering into the equations of motion of distributed systems. Because the motion of distributed systems is described in terms of partial differential equations, these parameters are in general continuous functions of the spatial variables. For vibrating systems, these parameters ordinarily represent the mass, stiffness and damping distributions. In this paper, these distributions are expanded in terms of finite series of known functions of the spatial variables multiplied by undetermined coefficients. It is assumed that the nature of the equations of motion is known and that a limited number of eigenvalues and eigenfunctions is identified in advance. Use is then made of the least squares method, in conjunction with the eigenfunctions' orthogonality, to compute the undetermined coefficients, thus identifying the system distributed parameters. A method for the identification of the eigensolution is also presented. The procedure for the identification of the eigensolution and of the system parameters is demonstrated via a numerical example.

Baruh, H.↗
NASA NTRS1970

Rapid method for interconversion of binary and decimal numbers

Decoding tree consisting of 40-bit semiconductor read-only memories interconverts binary and decimal numbers 50 to 100 times faster than current methods. Decimal-to-binary conversion algorithm is based on a divided-by-2 iterative equation, binary-to-decimal conversion algorithm utilizes multiplied-by-2 iterative equation.

Lim, R. S.↗
NASA NTRS1980

Specifying and calibrating instrumentations for wideband electronic power measurements

The wideband electric power measurement related topics of electronic wattmeter calibration and specification are discussed. Tested calibration techniques are described in detail. Analytical methods used to determine the bandwidth requirements of instrumentation for switching circuit waveforms are presented and illustrated with examples from electric vehicle type applications. Analog multiplier wattmeters, digital wattmeters and calculating digital oscilloscopes are compared. The instrumentation characteristics which are critical to accurate wideband power measurement are described.

Lesco, D. J.↗