Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “analytic continuation”

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 55 records · Page 3

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↗

Trajectory Calculator for Finite-Radius Cutter on a Lathe

A computer program calculates the two-dimensional trajectory (radial vs. axial position) of a finite-radius-of-curvature cutting tool on a lathe so as to cut a workpiece to a piecewise-continuous, analytically defined surface of revolution. (In the original intended application, the tool is a diamond cutter, and the workpiece is made of a crystalline material and is to be formed into an optical resonator disk.) The program also calculates an optimum cutting speed as F/L, where F is a material-dependent empirical factor and L is the effective instantaneous length of the cutting edge.

Savchenkov, Anatoliy↗

Investigation of Aerodynamic and Aerodynamic and Radiometric Land Surface Temperatures

The overall goal of the project was to reconcile the difference between T(sub s,r) and T(sub aero), while maintaining consistency within models and with theory and data. The project involved collaboration between researchers at Bucknell University, Boston University, University of mode Island, and the USDNARS Hydrology Laboratory. This report focuses on the work done at Bucknell, which used an analytical continuous-source flux model developed by Crago (1998), based on work by Brutsaert and Sugita (1996) to generate fluxes at all levels of the canopy. Named ALARM [Analytical Land- Atmosphere-Radiometer Model] by Suleiman and Crago (2002), the model assumes the foliage has an exponential vertical temperature profile. The same profile is felt by the within-canopy turbulence and 'seen" by a radiometer viewing the surface from any zenith view angle. ALARM converts radiometric surface temperatures taken from any view angle into a clearly-defined version of Taero called the equivalent isothermal surface temperature T(sub s,j), and then calculates the sensible heat flux H using Monin-Obukhov similarity theory. This allows remotely sensed Ts,r measurements to be used to produce high quality sensible and latent heat flux estimates, or to validate or update the surface temperature produced by SVATs in climate or mesoscale models.

Crago, Richard D.↗

Downward continuation of the free-air gravity anomalies to the ellipsoid using the gradient solution and terrain correction: An attempt of global numerical computations

The formulas for the determination of the coefficients of the spherical harmonic expansion of the disturbing potential of the earth are defined for data given on a sphere. In order to determine the spherical harmonic coefficients, the gravity anomalies have to be analytically downward continued from the earth's surface to a sphere-at least to the ellipsoid. The goal is to continue the gravity anomalies from the earth's surface downward to the ellipsoid using recent elevation models. The basic method for the downward continuation is the gradient solution (the g sub 1 term). The terrain correction was also computed because of the role it can play as a correction term when calculating harmonic coefficients from surface gravity data. The fast Fourier transformation was applied to the computations.

Wang, Y. M.↗

On the Gibbs phenomenon 4: Recovering exponential accuracy in a sub-interval from a Gegenbauer partial sum of a piecewise analytic function

We continue our investigation of overcoming Gibbs phenomenon, i.e., to obtain exponential accuracy at all points (including at the discontinuities themselves), from the knowledge of a spectral partial sum of a discontinuous but piecewise analytic function. We show that if we are given the first N Gegenbauer expansion coefficients, based on the Gegenbauer polynomials C(sub k)(sup mu)(x) with the weight function (1 - x(exp 2))(exp mu - 1/2) for any constant mu is greater than or equal to 0, of an L(sub 1) function f(x), we can construct an exponentially convergent approximation to the point values of f(x) in any subinterval in which the function is analytic. The proof covers the cases of Chebyshev or Legendre partial sums, which are most common in applications.

Gottlieb, David↗

The reduction correction in North America

An inverse Poisson integral technique was used to determine a gravity field on the geoid which, when continued by analytic free space methods to the topographic surface, agrees with the observed field. The computation is performed in three stages, each stage refining the previous solution using data at progressively increasing resolution (1 x 1 deg, 5 x 5', 5/8 x 5/8') from a decreasing area of integration. Reduction corrections are computed at 5/8 x 5/8' granularity by differencing the geoidal and surface values, smoothed by low-pass filtering and sub-sampled at 5' intervals. The 1 x 1 deg averages of the reduction corrections thus obtained for 172 1 x 1 deg squares in western North America are discussed. The 1 x 1 deg mean reduction corrections are predominantly positive, varying from -3 to +15 mgal, with values in excess of 5 mgal for 26 squares. Their mean and rms values are +2.4 and 3.6 mgal respectively and they correlate well with the mean terrain corrections. The mean and rms contributions from the three stages of computation are: 1 x 1 deg stage +0.15 and 0.7 mgal; 5 x 5' stage + 1.0 and 1.6 mgal; and 5/8 x 5/8' stage +1.3 and 1.8 mgal. These results reflect a tendency for the contributions to become larger and more systematically positive as the wavelengths involved become shorter. The results are discussed in terms of two mechanisms; the first is a tendency for the absolute values of both positive and negative anomalies to become larger when continued downwards and, the second, a non-linear rectification, due to the correlation between gravity anomaly and topographic height, which results in the values continued to a level surface being systematically more positive than those on the topography.

Martzen, P. D.↗

PYROLYTIC GRAPHITE ROCKET THRUST CHAMBER INVESTIGATION PROJECT. MONTHLY LETTER

This is the Monthly Letter Progress Report for the eighth month of a twelve month technical study directed toward the following objectives: A. Determining analytically and experimentally the feasibility and attractiveness of using free-standing pyrolytic graphite for a radiation cooled liquid rocket thrust chamber and exit nozzle for space application. B. Establishing through analysis and experimental demonstration, a set of design principles to facilitate the design of a reliable pyrolytic graphite thrust chamber assembly for a specific range of operating conditions. C. Demonstrate by means of test firings of complete thrust chamber assemblies, the validity of the design data and principles developed during the program. Preparations have been completed for test firing the two 100 lb. thrust Boron Pyralloy chambers at Marquardt during the first week in March 1962. These two chambers will first be static pressure checked with water to a pressure of 135 psig. The chambers will then be test fired with N204/N2H4-MMH propellants . Test samples of eleven different pyrolytic carbide and nitride materials have been ordered for high temperature oxidation tests and should be delivered by the end of March for evaluation in the plasma torch facility. The specific materials ordered in free- standing condition in thicknesses from 0.030 inch to 0.065 inch were as follows: (a) Pyrolytic graphite (on hand) (b) Pyrolytic graphite-boron alloy (on hand) (c) Pyrolytic graphite-tungsten alloy (d)Boron nitride (e) Titanium nitride (f) Silicon carbide (g) Titanium carbide (h) Zirconium carbide (i) Hafnium carbide (j) Niobium carbide (k) Tantalum carbide (l) High density (1 1-90) graphite (m) High density graphite (heat treated) These materials were ordered from Raytheon, HTM, and American Meta 1 Products . A summary of the available data on the oxidation rates of pyrolytic materials under slow air and rocket motor conditions is shown in Figure 1. High Temperature Materials, Inc. has been contacted relative to providing Pyrographite thrust chambers with extra-thick walls . They believe they can deliver chambers with substantially thicker walls than were possible a short time ago. Such a configuration, with acceptable residual stresses, would provide a longer firing life since it is ultimately limited by the oxidation rate of the chamber wall material. Stress and structural studies have continued. The analytical stress studies point out some of the problems associated with the experimental determination of the elastic constants of this highly anisotropic material. Further analysis may shed light on improved techniques for making these experimental determinations. The large potential weight saving available with the use of pyrolytic graphite as a thrust chamber wall material in a radiation cooled rocket motor has been calculated and is shown graphically in Figure 2. Figure A presents the program progress schedule and the actual expenditures as of 7 February 1962.

ROCKET ENGINE↗

Continued development and correlation of analytically based weight estimation codes for wings and fuselages

The implementation of the changes to the program for Wing Aeroelastic Design and the development of a program to estimate aircraft fuselage weights are described. The equations to implement the modified planform description, the stiffened panel skin representation, the trim loads calculation, and the flutter constraint approximation are presented. A comparison of the wing model with the actual F-5A weight material distributions and loads is given. The equations and program techniques used for the estimation of aircraft fuselage weights are described. These equations were incorporated as a computer code. The weight predictions of this program are compared with data from the C-141.

Mullen, J., Jr.↗

Design and true Reynolds number 2-D testing of an advanced technology airfoil

A NASA-industry program has been conducted to determine the accuracy of available 2-D airfoil analysis procedures over a wide range of Reynolds numbers. The program also served to develop and demonstrate effective wind tunnel model designs for use in a cryogenic environment. A Lockheed design, CRYO 12X, supercritical, shockfree airfoil was configured using a continuous curvature analytical definition of the ordinates. Test results show a very close ordinate tolerance was necessary to realize the intended pressure distribution. Correlation of test with Korn-Garabedian 2-D analysis pressure data were generally good. GRUMFOIL analysis with a sidewall correction gave a better correlation.

Reaser, J. S.↗

Dispersal of gases generated near a lunar outpost

The dispersal of gases evolved by prospective lunar base operations is presently modeled analytically using continuous (mining and habitat-venting) and impulsive (primarily rocket exhaust) injections of gases. In the case of impulsive injection, the neutral atmosphere and associated ionosphere both decay on time-scales of about 20 min; in that of continuous injection, the atmosphere near the base grows and achieves a steady state after about 20 min. Both direct and diffusive transport mechanisms are considered, and it is concluded that for the injection rates and assumptions presently employed the artificial lunar atmospheres produced are not significantly detrimental to astronomical observations and high-vacuum materials processing operations.

Burns, Jack O.↗

Reliability Considerations for Composite Overwrapped Pressure Vessels on Spacecraft

Composite Overwrapped Pressure Vessels (COPVs) are used to store gases under high pressure onboard spacecraft. These are used for a variety of purposes such as propelling liquid fuel etc, Kevlar, glass, Carbon and other more recent fibers have all been in use to overwrap the vessels. COPVs usually have a thin metallic liner with the primary purpose of containing the gases and prevent any leakage. The liner is overwrapped with filament wound composite such as Kevlar, Carbon or Glass fiber. Although the liner is required to perform in the leak before break mode making the failure a relatively benign mode, the overwrap can fail catastrophically under sustained load due to stress rupture. It is this failure mode that is of major concern as the stored energy of such vessels is often great enough ta cause loss of crew and vehicle. The present paper addresses some of the reliability concerns associated specifically with Kevlar Composite Overwrapped Pressure Vessels. The primary focus of the paper is on how reliability of COPV's are established for the purpose of deciding in general their flight worthiness and continued use. Analytical models based on existing design data will be presented showing how to achieve the required reliability metric to the end of a specific period of performance. Uncertainties in the design parameters and how they affect reliability and confidence intervals will be addressed as well. Some trade studies showing how reliability changes with time during a program period will be presented.

Murthy, Pappu L. N.↗

Benchmarking Bayesian Optimization Frameworks and Acquisition Strategies for Materials Discovery and Autonomous Laboratories

Bayesian optimization (BO) can accelerate materials discovery by guiding expensive experiments toward the most promising processing conditions. We systematically compare five BO surrogate and framework combinations (Gaussian processes in Ax, Gaussian processes and Monte-Carlo neural networks in BayBE, random forests in Lolopy, and tree-structured Parzen (TPE) estimators in Hyperopt) on three benchmarks that mimic common materials design tasks (a discrete solid-electrolyte composition space, a hybrid discrete/continuous laminate-composite design problem solved with micromechanics modeling, and the continuous Ishigami analytic function which is a standard optimization benchmark). Each BO surrogate is paired with posterior mean, probability of improvement, and expected improvement acquisition functions and run for 100 trials from randomized initial samples with uniform random search providing a control. Across five random seeds per setting, BayBE’s Gaussian-process surrogate with expected improvement consistently reached ≥95 % of the known optimum in the fewest evaluations, while Lolopy’s random forest matched or exceeded GP performance on purely categorical or mixed spaces at a higher computational cost. Posterior mean alone often stagnated at local optima, underscoring the need for exploration, whereas probability and expected improvement balanced exploration and exploitation leading to better optimization in fewer trials. Execution times ranged from milliseconds for TPE to minutes for neural-network and random-forest surrogates. These results establish baseline expectations for BO in automated materials laboratories and highlight expected improvement with Gaussian processes as a reliable first choice, with random forests offering a strong alternative when categorical variables dominate. The benchmark suite and code are released to facilitate future surrogate, acquisition, and constraint-handling research in data-driven materials optimization.

Bayesian optimization↗

Solidification interface shape for continuous casting in an offset mold - Two analytical methods

A solution method for finding the unknown solidification interface in manufacturing slab ingots as a continuous casting is presented, which involves a product solution in the potential plane and the use of conjugate harmonic functions. It is argued that the method may be more direct for some geometries than the Cauchy boundary value method. Moreover, the usefulness of the Cauchy boundary value method is demonstrated through the example of a nonsymmetric horizontal mold where the walls are offset to support the lower ingot boundary.

Siegel, R.↗

On the Behavior of Three-dimensional Wave Packets in Viscously Spreading Mixing Layers

We consider analytically the evolution of a three-dimensional wave packet generated by an impulsive source in a mixing layer. The base flow is assumed to be spreading due to viscous diffusion. The analysis is restricted to small disturbances (linearized theory). A suitable high-frequency ansatz is used to describe the packet; the key elements of this description are a complex phase and a wave action density. It is found that the product of this density and an infinitesimal material volume convecting at the local group velocity is not conserved: there is a continuous interaction between the base flow and the wave action. This interaction is determined by suitable mode-weighted averages of the second and fourth derivatives of the base-flow velocity profile. Although there is some tendency for the dominant wave number in the packet to shift from the most unstable value toward the neutral value, this shift is quite moderate. In practice, wave packets do not become locally neutral in a diverging base flow (as do instability modes), therefore, they are expected to grow more suddenly than pure instability modes and do not develop critical layers. The group velocity is complex; the full significance of this is realized by analytically continuing the equations for the phase and wave action into a complex domain. The implications of this analytic continuation are discussed vis-a-vis the secondary instabilities of the packet: very small scale perturbations on the phase can grow very rapidly initially, but saturate later because most of the energy in these perturbations is convected away by the group velocity. This remark, as well as the one regarding critical layers, has consequences for the nonlinear theories.

Balsa, Thomas F.↗

System Definition Phase and Acquisition Phase Project Plan for Small Astronomy Satellite SAS-D

This Project Plan for the Small Astronomy Satellite (SAS-D) is Goddard Space Flight Center's plan for accomplishing the project. It is issued in compliance with NASA Handbook NHB 7121.2 and NMI 7120.1, dated May 4, 1970, Subject: Approval and Implementation of Office of Space Science and Applications (OSSA) Research and Development (R&D) Projects. The Project Plan, as approved by GSFC executive management and NASA Headquarters, is the project manager's charter and defines his specific plan of action. System definition phase activities of the SAS-D Project will be a continuation of the analytical work that established the overall feasibility of the project. The system concept derived from that analytical work is defined in the Phase A (Mission Definition Phase) technical report. System definition phase activities will include detailed engineering analysis and preliminary breadboarding necessary to establish a single project approach from the feasible approaches. This document will be revised at the end of the system definition phase to provide greater detail on acquisition phase planning and to specify a specific system design.

Donald A Krueger↗

Analytical design and performance studies of nuclear furnace tests of small nuclear light bulb models

Analytical studies were continued to identify the design and performance characteristics of a small-scale model of a nuclear light bulb unit cell suitable for testing in a nuclear furnace reactor. Emphasis was placed on calculating performance characteristics based on detailed radiant heat transfer analyses, on designing the test assembly for ease of insertion, connection, and withdrawal at the reactor test cell, and on determining instrumentation and test effluent handling requirements. In addition, a review of candidate test reactors for future nuclear light bulb in-reactor tests was conducted.

Latham, T. S.↗

Large space structures controls research and development at Marshall Space Flight Center: Status and future plans

Work performed in Large Space Structures Controls research and development program at Marshall Space Flight Center is described. Studies to develop a multilevel control approach which supports a modular or building block approach to the buildup of space platforms are discussed. A concept has been developed and tested in three-axis computer simulation utilizing a five-body model of a basic space platform module. Analytical efforts have continued to focus on extension of the basic theory and subsequent application. Consideration is also given to specifications to evaluate several algorithms for controlling the shape of Large Space Structures.

Buchanan, H. J.↗

Studies of convection in a solidifying binary mixture at reduced gravity

A great deal of interest was generated recently in the possibility of producing new materials in the reduced gravity environment provided during the forthcoming missions of Spacelab. The range of possibilities extend from producing large crystals of uniform properties to manufacturing materials with unique properties. Most of these processes involve the solidification of materials from the liquid state. Convective motions within the liquid during solidification can influence the local material composite and the shape of the solid-liquid interface which may result in solids with non-uniform properties and crystal defects. The microgravity environment of Spacelab is being viewed as one in which the buoyancy forces are eliminated so that convection driven by thermal gradients does occur, resulting in an improved solidification process. However, convection may occur for other reasons and whether convection is negligible or not during solidification constitutes processing in low-gravity environment. Little information exists presently on convection during solidification under such circumstances. A continuation of an analytical investigation into the nature of convective motion in a binary liquid layer due to surface tension forces during its solidification is reported. The onset of convection will be determined through a stability analysis which is described.

Antar, B. N.↗