Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interior 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.

At least 325 records · Page 18
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Lasergate: A windowless gas target for enhanced laser preheat in magnetized liner inertial fusion

At the Z Facility at Sandia National Laboratories, the magnetized liner inertial fusion (MagLIF) program aims to study the inertial confinement fusion in deuterium-filled gas cells by implementing a three-step process on the fuel: premagnetization, laser preheat, and Z-pinch compression. In the laser preheat stage, the Z-Beamlet laser focuses through a thin polyimide window to enter the gas cell and heat the fusion fuel. However, it is known that the presence of the few μm thick window reduces the amount of laser energy that enters the gas and causes window material to mix into the fuel. These effects are detrimental to achieving fusion; therefore, a windowless target is desired. The Lasergate concept is designed to accomplish this by “cutting” the window and allowing the interior gas pressure to push the window material out of the beam path just before the heating laser arrives. In this work, we present the proof-of-principle experiments to evaluate a laser-cutting approach to Lasergate and explore the subsequent window and gas dynamics. Further, an experimental comparison of gas preheat with and without Lasergate gives clear indications of an energy deposition advantage using the Lasergate concept, as well as other observed and hypothesized benefits. While Lasergate was conceived with MagLIF in mind, the method is applicable to any laser or diagnostic application requiring direct line of sight to the interior of gas cell targets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Source term method for binary neutron stars initial data

The initial condition problem for a binary neutron star system requires a Poisson equation solver for the velocity potential with a Neumann-like boundary condition on the surface of the star. Difficulties that arise in this boundary value problem are: (a) the boundary is not known a priori, but constitutes part of the solution of the problem; (b) various terms become singular at the boundary. In this work, we present a new method to solve the fluid Poisson equation for irrotational/spinning binary neutron stars. The advantage of the new method is that it does not require complex fluid surface fitted coordinates and it can be implemented in a Cartesian grid, which is a standard choice in numerical relativity calculations. This is accomplished by employing the source term method proposed by Towers, where the boundary condition is treated as a jump condition and is incorporated as additional source terms in the Poisson equation, which is then solved iteratively. The issue of singular terms caused by vanishing density on the surface is resolved with an additional separation that shifts the computation boundary to the interior of the star. We present two-dimensional tests to show the convergence of the source term method, and we further apply this solver to a realistic three-dimensional binary neutron star problem. By comparing our solution with the one coming from the initial data solver cocal, we demonstrate agreement to approximately 1%. We report our method can be used in other problems with non-smooth solutions like in magnetized neutron stars.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1988

A new diagnostic method for separating airborne and structureborne noise radiated by plates with applications for propeller driven aircraft

The anticipated application of advanced turboprop propulsion systems is expected to increase the interior noise of future aircraft to unacceptably high levels. The absence of technically and economically feasible noise source-path diagnostic tools has been a prime obstacle in the development of efficient noise control treatments for propeller-driven aircraft. A new diagnostic method that permits the separation and prediction of the fully coherent airborne and structureborne components of the sound radiated by plates or thin shells has been developed. Analytical and experimental studies of the proposed method were performed on an aluminum plate. The results of the study indicate that the proposed method could be used in flight, and has fewer encumbrances than the other diagnostic tools currently available.

Mcgary, Michael C.↗
NASA NTRS1981

Fully implicit shock tracking

A tracking method for solving hyperbolic free boundary problems arising in fluid dynamics is presented. The principal feature of the new approach is that the implicit form is used to simultaneously treat both interior points and boundary conditions. In particular, the location of the free boundary (shock) surface is treated implicitly and coupled with all other unknowns. The method is presented in the context of unsteady one-dimensional flow in a variable area duct with an internal shock wave. The fully implicit method and other strategies for advancing the shock are compared for computing a steady solution via a time asymptotic approach. Issues regarding extension of the method to multiple dimensions are also discussed.

Bell, J. B.↗
DOE OSTI2025

Mitigation of Pitting on Nitrogen-Doped Niobium Surfaces through Two-Step Electropolishing

Surface processing of niobium superconducting RF cavities used in particle accelerators is crucial for determining their performance. The nitrogen (N)-doping method, developed at Fermilab, enhances the cavities’ quality factor by doping the cavity’s interior surface with N atoms at 800 °C. After N-doping, electropolishing (EP) is employed to remove the niobium nitride phases and attain the desired concentration of interstitial N atoms. This study explores the effects of EP conditions on N-doped niobium samples. The samples were electropolished at 5 to 40 °C and voltages between 3 and 22 V. Different EP voltage ranges produced two types of pits, classified based on their sizes. Larger pits, linked to higher voltages, are also affected by the reaction rate. Pitting could be reduced at a lower temperature and adequate EP voltage. However, to effectively minimize gas evolution forming large-sized pitting, peeling off the top nitride layer before EP might be beneficial. This was achieved through precise low-voltage etching, forming a two-step EP process. The process involves removing the top 100–200 nm layer at a low voltage, followed by standard EP to remove the required material. This two-step EP process effectively mitigated the pitting risk associated with gas evolution on N-doped surfaces.

43 PARTICLE ACCELERATORS↗
NASA NTRS1973

The gravitational fields of the major planets.

Review of the constraints imposed on models of the interiors of the major planets by the nonspherical components of their gravitational fields. Several methods for determining these nonspherical components are described and evaluated. The nature of the constraints and their effect on the models are discussed.

Peale, S. J.↗
NASA NTRS1978

An improved higher order panel method for linearized supersonic flow

An improved higher order panel method for linearized supersonic flow is described. Each panel, defined by four points on the surface, is divided into eight subpanels in such a way that all subpanel and panel edges are contiguous. By prescribing a quadratic distribution of the doublet on each subpanel, the doublet strength is made strictly continuous on the paneled surface. A linear source distribution is also used. Numerical results are smoother and in better agreement with experiment than the previous method with less strict continuity. A brief discussion of superinclined panels used to eliminate interior interference in nacelles is included.

Ehlers, F. E.↗
NASA NTRS1984

Transformation of two and three-dimensional regions by elliptic systems

The research during this period continued to expand the class of numerical algorithms that can be accurately and efficiently implemented on overlapping grids. Whereas previous calculations have been used to solve elliptic equations and to find the steady-state solution of parabolic equations, the present work is aimed towards developing time-accurate solution techniques for parabolic and hyperbolic equations. The primary difficulty here is in the correct treatment of the interior boundary nodes that must be updated at each iteration. The implementation of explicit methods is straightforward. However, the common practice of lagging these values when using an implicit methods leads to inconsistencies in the difference equation. One way to avoid this problem is to alternately calculate with an implicit and an explicit method on each subgrid. With this procedure, the explicit method generates boundary values at the next time level which are then used by the implicit step. It can be shown that when a backward implicit method is combined with a forward explicit method, the composite method is second order accurate and unconditionally stable for linear problems. A second area in which progress can be reported is in the distribution of grid points on curves and surfaces.

Mastin, E. W.↗
NASA NTRS1985

Chemical evolution of the Earth: Equilibrium or disequilibrium process?

To explain the apparent chemical incompatibility of the Earth's core and mantle or the disequilibrium process, various core forming mechanisms have been proposed, i.e., rapid disequilibrium sinking of molten iron, an oxidized core or protocore materials, and meteorite contamination of the upper mantle after separation from the core. Adopting concepts used in steady state thermodynamics, a method is devised for evaluating how elements should distribute stable in the Earth's interior for the present gradients of temperature, pressure, and gravitational acceleration. Thermochemical modeling gives useful insights into the nature of chemical evolution of the Earth without overly speculative assumptions. Further work must be done to reconcile siderophile elements, rare gases, and possible light elements in the outer core.

Sato, M.↗
NASA NTRS1985

On optimal strategies in event-constrained differential games

Combat games are formulated as zero-sum differential games with unilateral event constraints. An interior penalty function approach is employed to approximate optimal strategies for the players. The method is very attractive computationally and possesses suitable approximation and convergence properties.

Heymann, M.↗
NASA NTRS1991

Full complex modulation with two one-parameter SLMs

Although the action of a spatial light modulator (SLM) is usually restricted to certain locations on the operating curve of the complex plane, NASA is planning to use architectures that allow two continuously variable SLMs to function jointly so as to access the full interior of a closed curve in the complex plane. This paper describes three fundamental methods for attaining full complex modulation. The mathematics for two of these methods is presented, and signal decomposition in their terms is outlined.

Juday, Richard D.↗
NASA NTRS1994

Nucleation of Crystals From Solution in Microgravity (USML-1 Glovebox (GBX) Investigation)

A new method for initiating nucleation from solutions in microgravity which avoids nucleation on container walls and other surfaces is described. This method consists of injecting a small quantity of highly concentrated, heated solution into the interior of a lightly supersaturated, cooler host gowth solution. It was tested successfully on USML-I, producing a large number of LAP crystals whose longest dimension averaged 1 mm.

Kroes, Roger L.↗
NASA NTRS2004

Onset of Convection in Ice I with Composite Newtonian and Non-Newtonian Rheology: Application to the Icy Galilean Satellites

Ice I exhibits a complex rheology at temperature and pressure conditions appropriate for the interiors of the ice I shells of Europa, Ganymede, and Callisto. We use numerical methods and existing parameterizations of the critical Rayleigh number to determine the conditions required to trigger convection in an ice I shell with the stress-, temperature- and grain size- dependent rheology measured in laboratory experiments by Goldsby and Kohlstedt [2001]. The critical Rayleigh number depends on the ice grain size and the amplitude and wavelength of temperature perturbation issued to an initially conductive ice I shell. If the shells have an assumed uniform grain size less than 0.4 mm, deformation during initial plume growth is accommodated by Newtonian volume diffusion. If the ice grain size is between 0.4 mm and 3 cm, deformation during plume growth is accommodated by weakly non-Newtonian grain boundary sliding, where the critical ice shell thickness for convection depends on the amplitude of temperature perturbation to the _0.5 power. If the ice grain size exceeds 2 cm, convection can not occur in the ice I shells of the Galilean satellites regardless of the amplitude or wavelength of temperature perturbation. If the grain size in a convecting ice I shell evolves to effective values greater than 2 cm, convection will cease. If the ice shell has a grain size large enough to permit flow by dislocation creep, the ice is too stiff to permit convection, even in the thickest possible ice I shell. Consideration of the composite rheology implies that estimates of the grain size in the satellites and knowledge of their initial thermal states are required when judging the convective instability of their ice I shells.

Barr, A. C.↗
NASA NTRS1983

On a finite-difference method for solving transient viscous flow problems

A method has been developed to solve the unsteady, compressible Navier-Stokes equation with the property of consistency and the ability of minimizing the equation stiffness. It relies on innovative extensions of the state-of-the-art finite-difference techniques and is composed of: (1) the upwind scheme for split-flux and the central scheme for conventional flux terms in the inviscid and viscous regions, respectively; (2) the characteristic treatment of both inviscid and viscous boundaries; (3) an ADI procedure compatible with interior and boundary points; and (4) a scalar matrix coefficient including viscous terms. The performance of this method is assessed with four sample problems; namely, a standing shock in the Laval duct, a shock reflected from the wall, the shock-induced boundary-layer separation, and a transient internal nozzle flow. The results from the present method, an existing hybrid block method, and a well-known two-step explicit method are compared and discussed. It is concluded that this method has an optimal trade-off between the solution accuracy and computational economy, and other desirable properties for analyzing transient viscous flow problems.

Li, C. P.↗
NASA NTRS2000

An Assessment of the Ozone Loss During the 1999-2000 SOLVE Campaign

Ozone observations from ozonesondes, the DIAL and AROTEL lidars aboard the DC-8, in situ ozone measurements from the ER-2 and satellite ozone measurements from POAM were used to assess ozone loss during the SOLVE 1999-2000 campaign. We compare three different methods of computing the ozone loss. The first method simply compares the time sequence of ozonesondes taken at the same station inside the vortex from December through the end of March. In the second method, ozonesondes from a variety of stations are compared using a variant on the Match technique. This method uses short (approx. 5-10 day) forward diabatic trajectories to connect various sonde launches. In the third method, the measurements are simply injected into a diabatic trajectory model and carried forward in time from December 1 to March 16. Over 60,000 individual measurements were used in the last calculation. Again, ozone loss is estimated by comparing vortex interior measurements made early in the campaign with those made later in the campaign. The diabatic nature of the second and third methods calculation presumably corrects for the normal increase in ozone within the vortex due to downward advection. The three methods agree that the largest ozone loss occurs between 400 and 460 K potential temperatures (approx. 16-20 km) with slightly over 1.5 ppmv lost over the winter period. Between 460 K and 500 K (approx. 22 km) net ozone loss is less than 0.8 ppmv. From 500K to 600K (26 km) net loss is less than 0.5 ppmv.

Schoeberl, M. R.↗
DOE OSTI2023

Ultrasonic measurement of temperature distributions in extreme environments: Electrical power plants testing in utility-scale steam generators

Thermal heterogeneities within energy conversion and storage, material processing, nuclear processes, aerospace, and military applications are often inaccessible to characterization by insertion sensors. When sensor deployment is possible, conventional pointwise temperature probes quickly degrade when inserted into harsh environments typical of such processes. We developed spatially-resolved ultrasonic thermometry to noninvasively measure the spatial distributions of thermal properties in such applications, even when sizable thermal gradients are present. Our method divides the path of ultrasonic propagation into segments bound by echogenic features, which create echoes in pulse-echo mode, encoding the information about interior temperature distributions. We use the acquired ultrasonic responses to estimate the internal temperature distributions by solving an inverse problem or concatenating segmental estimates. This work describes the implementation and industrial testing of the developed method at a coal-fired electrical power generation plant. We inserted an echogenically segmented Inconel 625 waveguide into the combustion zone of the utility-scale boiler and continuously acquired ultrasonic data while keeping sensitive components away from the damaging combustion environment. The accuracy of the time-dependent temperature distributions reconstructed from the ultrasonic measurements was comparable to that of thermocouples. The resiliency of ultrasonic thermometry to harsh combustion conditions was far superior to conventional insertion sensors. The measurements obtained during plant operation captured daily steam generation cycles in response to changing customer demand and intermittent contributions of renewable power sources to the power grid. These measurements have revealed new insights into the relationship between the dynamic power generation load and the conditions inside the steam generator. As a result, the successful industrial testing of spatially-resolved ultrasonic thermometry in solids indicates that the developed technology has matured to become an attractive alternative to conventional sensing in solving challenging problems of long-term thermal characterizations in extreme environments.

42 ENGINEERING↗
NASA NTRS1986

Multigrid-sinc methods

A Galerkin method using Whittaker cardinal or 'sinc' functions as basis functions is described for the solution of boundary-value problems. When the solution is analytic in the interior of the domain, the error of approximation using 2N + 1 points is O(e exp /-gamma sq rt N/) even if derivatives of the solution are singular at the boundaries. A multigrid method with overall complexity O(N log N) is used to solve the discrete equations. This paper contains a description of the multigrid-sinc algorithm along with some preliminary numerical results for two-point boundary-value problems.

Schaffer, Steve↗
NASA NTRS1973

Survey of computational methods for three-dimensional supersonic inviscid flows with shocks

Hyperbolic systems of partial differential equations governing supersonic inviscid flows are discussed and analyzed. Finite-difference analogues for integrating these systems in the interior of fluid domains are described from two points of view: a differential form approach and an integral form approach. The algorithms presented are analyzed for stability and accuracy. The concept of time splitting is discussed and supplied to these methods to achieve increased numerical efficiency. Techniques for treating conditions at the boundaries of the fluid domain and shock-wave discontinuities at surfaces within the domain are described.

Maccormack, R. W.↗