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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRS1985

An implicit time-marching method for studying unsteady flow with massive separation

A fully implicit time-marching method is developed such that all spatial derivatives are approximated using central differences, but no use is made of any artificial dissipation. The numerical method solves the discretized equations using Alternating Direction Implicit-Block Gaussian Elimination technique. The method is implemented in the unsteady analysis, which solves the incompressible Navier-Stokes equations in terms of vorticity and stream function in generalized orthogonal coordinates. A clustered conformal C-grid is employed, and every effort is made to resolve the various length scales in the flow problem. The metric discontinuity at the branch-cut is treated appropriately using analytic continuation. Introduction of the BGE reordering permits implicit treatment of the branch cut in the numerical method. The vorticity singularity at the cusped trailing edge is also appropriately treated. This accurate and efficient implicit method is used to study flow at Re = 1000, past a 12-percent thick symmetric Joukowski airfoil at high angle of attack 30 and 53 deg.

Osswald, G. A.↗
NASA NTRS1987

Hard-sphere fluid in infinite dimensions

The exact Mayer series is derived for the pressure, density, and equation of state of a hard-sphere fluid in the limit of infinite space dimensions D. The Mayer series can be analytically continued into the full (cut) activity plane, and there is no sign of a phase transition. The quantum-mechanical system is also treated. For D = infinity, the fluid behaves like an ideal Bose gas and undergoes a Bose-Einstein condensation.

Wyler, D.↗
NASA NTRS1988

Filling the expertise gap

Aeroelasticity, structures, and stability-and-control specialists can work in concerted fashion during the early design phases of future aircraft to achieve active control of naturally unstable configurations. In order to take full advantage of active control, attention must be given by designers to control-law synthesis, and to tools for the efficient synthesis and analysis of complex flexible-aircraft control systems. Analysts must consider interfaces among unsteady aerodynamics, structures, and control theory, as explored by the theory of analytic continuation for unsteady aerodynamics.

Abel, Irving↗
NASA NTRS1989

Fourier transform methods in local gravity modeling

New algorithms were derived for computing terrain corrections, all components of the attraction of the topography at the topographic surface and the gradients of these attractions. These algoriithms utilize fast Fourier transforms, but, in contrast to methods currently in use, all divergences of the integrals are removed during the analysis. Sequential methods employing a smooth intermediate reference surface were developed to avoid the very large transforms necessary when making computations at high resolution over a wide area. A new method for the numerical solution of Molodensky's problem was developed to mitigate the convergence difficulties that occur at short wavelengths with methods based on a Taylor series expansion. A trial field on a level surface is continued analytically to the topographic surface, and compared with that predicted from gravity observations. The difference is used to compute a correction to the trial field and the process iterated. Special techniques are employed to speed convergence and prevent oscillations. Three different spectral methods for fitting a point-mass set to a gravity field given on a regular grid at constant elevation are described. Two of the methods differ in the way that the spectrum of the point-mass set, which extends to infinite wave number, is matched to that of the gravity field which is band-limited. The third method is essentially a space-domain technique in which Fourier methods are used to solve a set of simultaneous equations.

Harrison, J. C.↗
NASA NTRS1989

Thin-plate spline quadrature of geodetic integrals

Thin-plate spline functions (known for their flexibility and fidelity in representing experimental data) are especially well-suited for the numerical integration of geodetic integrals in the area where the integration is most sensitive to the data, i.e., in the immediate vicinity of the evaluation point. Spline quadrature rules are derived for the contribution of a circular innermost zone to Stoke's formula, to the formulae of Vening Meinesz, and to the recursively evaluated operator L(n) in the analytical continuation solution of Molodensky's problem. These rules are exact for interpolating thin-plate splines. In cases where the integration data are distributed irregularly, a system of linear equations needs to be solved for the quadrature coefficients. Formulae are given for the terms appearing in these equations. In case the data are regularly distributed, the coefficients may be determined once-and-for-all. Examples are given of some fixed-point rules. With such rules successive evaluation, within a circular disk, of the terms in Molodensky's series becomes relatively easy. The spline quadrature technique presented complements other techniques such as ring integration for intermediate integration zones.

Vangysen, Herman↗
NASA NTRS1989

Monte Carlo calculation of dynamical properties of the two-dimensional Hubbard model

A new method is introduced for analytically continuing imaginary-time data from quantum Monte Carlo calculations to the real-frequency axis. The method is based on a least-squares-fitting procedure with constraints of positivity and smoothness on the real-frequency quantities. Results are shown for the single-particle spectral-weight function and density of states for the half-filled, two-dimensional Hubbard model.

White, S. R.↗
NASA NTRS1993

Establishing a relationship between maximum torque production of isolated joints to simulate EVA ratchet push-pull maneuver: A case study

As manned exploration of space continues, analytical evaluation of human strength characteristics is critical. These extraterrestrial environments will spawn issues of human performance which will impact the designs of tools, work spaces, and space vehicles. Computer modeling is an effective method of correlating human biomechanical and anthropometric data with models of space structures and human work spaces. The aim of this study is to provide biomechanical data from isolated joints to be utilized in a computer modeling system for calculating torque resulting from any upper extremity motions: in this study, the ratchet wrench push-pull operation (a typical extravehicular activity task). Established here are mathematical relationships used to calculate maximum torque production of isolated upper extremity joints. These relationships are a function of joint angle and joint velocity.

Pandya, Abhilash↗
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↗
NASA NTRS2009

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↗
DOE OSTI2023

Demonstrating Computational Equivalence Between Continuous and Discrete Adjoint Methods by Calculating Time-Dependent Adjoint Solutions with Neutron Diffusion Models

The continuous adjoint method and the discrete adjoint method are two alternative approaches used to calculate adjoint solutions for adjoint systems. The continuous adjoint method derives adjoint equations analytically from continuous forward equations and then solves the adjoint equations either analytically or numerically in a discretized form whereas the discrete adjoint method calculates the adjoint solutions directly from the discretized forward equations. With regard to the methodology development and calculation procedure, distinct differences are well recognized between the two methods. For certain reasons, both methods are exclusively preferred and commonly used by different computational communities, but limited studies clarify the connections between the two adjoint methods from either of the communities. Herein, this paper demonstrates the computational equivalence between the continuous and discrete adjoint methods by investigating time-dependent adjoint solutions to the two-group neutron diffusion model in nuclear reactor analysis problems using both methods. Adjoint solutions can be used to estimate system parameters for reactor safety analysis. Appropriate final state conditions for the adjoint systems are specified in both of the methods, and the conditions are clarified with proper physical explanations. With the help of an event-based case study on neutron diffusion models, the accuracy of the time-dependent adjoint fluxes obtained from both methods is verified, and the pros and cons of both adjoint methods are examined. More importantly, the computational equivalence of both methods is demonstrated when they are applied to multigroup neutron diffusion systems. The advantage of calculating time-dependent adjoint fluxes by directly solving time-dependent adjoint systems rather than taking steady-state approximations as in common practice is also demonstrated.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2003

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

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

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↗
NASA NTRS1989

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

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↗
DOE OSTI2024

Automated Framework for Groundwater Monitoring Using DWT with LSTM and Transformers

Environmental monitoring is critical for safeguarding public health and ecological well-being. Traditional data structuring and workflow monitoring methods consume significant time and effort, hindering timely insights and effective decision-making. Our study addresses this challenge by presenting an AI framework that automates data cleaning, structuring, and modeling processes, specifically targeting applications in groundwater monitoring. By leveraging automation for data processing and model training, our framework establishes a novel and efficient paradigm for environmental monitoring, with its potential application to the vast network of over a hundred Department of Energy Environmental Management (DoE-EM) cleanup sites across the country. It analyzes data streams from a network of groundwater Internet-of-Things (IoT) sensors deployed at the Savannah River Site (SRS) for prediction modeling. This allows human experts to focus on analysis and decision-making, ultimately leading to better environmental outcomes.The framework employs multivariate time-series forecasting methods to study and model the behavior of varying chemical analytes. The continuous learning process is enabled by utilizing deep learning techniques. It allows the framework to become more nuanced in its analysis over time, adapting to the specific characteristics of the environmental site and the evolving nature of contaminant behavior. Deep learning models known for sequence modeling, LSTM, and Transformers are employed for time series forecasting. Data processing and structuring are essential components significantly impacting the final model's performance. This hypothesis was proven by presenting a comparative analysis of model performance with processed and unprocessed data. The feature engineering approach utilized was the Discrete Wavelet Transform, which works well with time series data.

Discrete Wavelet Transform (DWT)↗
NASA NTRS1978

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

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