Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Numerical techniques”

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 1,333 records · Page 74

Motorized Beam Alignment of a Commercial X-ray Diffractometer

X-ray diffraction (XRD) is a powerful analysis method that allows researchers to noninvasively probe the crystalline structure of a material. This includes the ability to determine the crystalline phases present, quantify surface residual stresses, and measure the distribution of crystallographic orientations. The Structures and Materials Division at the NASA Glenn Research Center (GRC) heavily uses the on-site XRD lab to characterize advanced metal alloys, ceramics, and polymers. One of the x-ray diffractometers in the XRD lab (Bruker D8 Discover) uses three different x-ray tubes (Cu, Cr, and Mn) for optimal performance over numerous material types and various experimental techniques. This requires that the tubes be switched out and aligned between experiments. This alignment maximizes the x-ray tube s output through an iterative process involving four set screws. However, the output of the x-ray tube cannot be monitored during the adjustment process due to standard radiation safety engineering controls that prevent exposure to the x-ray beam when the diffractometer doors are open. Therefore, the adjustment process is a very tedious series of blind adjustments, each followed by measurement of the output beam using a PIN diode after the enclosure doors are shut. This process can take up to 4 hr to perform. This technical memorandum documents an in-house project to motorize this alignment process. Unlike a human, motors are not harmed by x-ray radiation of the energy range used in this instrument. Therefore, using motors to adjust the set screws will allow the researcher to monitor the x-ray tube s output while making interactive adjustments from outside the diffractometer. The motorized alignment system consists of four motors, a motor controller, and a hand-held user interface module. Our goal was to reduce the alignment time to less than 30 min. The time available was the 10-week span of the Lewis' Educational and Research Collaborative Internship Project (LERCIP) summer internship program and the budget goal was $1200. In this report, we will describe our motorization design and discuss the results of its implementation.

Van Zandt, Noah R.↗

Further Development of Rotating Rake Mode Measurement Data Analysis

The Rotating Rake mode measurement system was designed to measure acoustic duct modes generated by a fan stage. After analysis of the measured data, the mode amplitudes and phases were quantified. For low-speed fans within axisymmetric ducts, mode power levels computed from rotating rake measured data would agree with the far-field power levels on a tone by tone basis. However, this agreement required that the sound from the noise sources within the duct propagated outward from the duct exit without reflection at the exit and previous studies suggested conditions could exist where significant reflections could occur. To directly measure the modes propagating in both directions within a duct, a second rake was mounted to the rotating system with an offset in both the axial and the azimuthal directions. The rotating rake data analysis technique was extended to include the data measured by the second rake. The analysis resulted in a set of circumferential mode levels at each of the two rake microphone locations. Radial basis functions were then least-squares fit to this data to obtain the radial mode amplitudes for the modes propagating in both directions within the duct. The fit equations were also modified to allow evanescent mode amplitudes to be computed. This extension of the rotating rake data analysis technique was tested using simulated data, numerical code produced data, and preliminary in-duct measured data.

Dahl, Milo D.↗

Development of a Perfectly Matched Layer Technique for a Discontinuous-Galerkin Spectral-Element Method

The perfectly matched layer (PML) technique is developed in the context of a high- order spectral-element Discontinuous-Galerkin (DG) method. The technique is applied to a range of test cases and is shown to be superior compared to other approaches, such as those based on using characteristic boundary conditions and sponge layers, for treating the inflow and outflow boundaries of computational domains. In general, the PML technique improves the quality of the numerical results for simulations of practical flow configurations, but it also exhibits some instabilities for large perturbations. A preliminary analysis that attempts to understand the source of these instabilities is discussed.

Spectral-Element↗

Impedance Eduction in a Duct Using the Linearized Euler Equations

An impedance eduction technique that is based upon a numerical solution to the linearized Euler equations coupled with the minimization of an objective function is presented. The linearized Euler equations are solved numerically using a cubic finite element method and the minimum of the objective function is obtained via a gradient based optimizer. The mean flow is allowed to have a gradient in two coordinate directions and the impedance of the liner may vary arbitrarily over the liner surface. A new regularized wall impedance boundary condition presented by Rienstra which contains the effects of the mean boundary layer thickness is implemented. Impedances are educed for a perforate over honeycomb liner using test data acquired in the Langley Grazing Flow Impedance Tube. The impedance of the liner is educed both with an upstream source and a downstream source. The method is validated by comparing the educed impedance to a benchmark method. The primary conclusion of the study is that the impedances educed for upstream and downstream sources are not well matched and that the primary effect of the boundary thickness is to reduce the resistance of the liner.

Watson, Willie R.↗

Solid Fuel Ignition and Extinction (SoFIE) Project on ISS

The Solid Fuel Ignition and Extinction (SoFIE) project studies ignition and flammability of solid spacecraft materials (fuels) in practical geometries and realistic atmospheric conditions. It is an experiment insert designed for use within the existing Combustion Integrated Rack (CIR). The CIR chamber provides a level of containment and permits testing at variable oxygen concentrations and pressures representative of current and planned NASA Space Exploration Atmospheres. The applications of SoFIE include: (1) Determining safer selection of cabin materials and validating NASA materials flammability selection using 1-g test protocols for low-gravity fires, (2) Improving understanding of early fire growth behavior, (3) Validating material flammability numerical models, (4) Determining optimal suppression techniques for burning materials by diluents, flow reduction, and venting, (5) Obtaining burning behavior of actual engineering materials planned for spacecraft, (6) Developing corresponding models of microgravity flame spread, flammability, and extinction, and use the results to improve normal gravity combustion models for terrestrial applications. The hardware permits a wide range of solid-material combustion and fire suppression studies. It supports multiple investigations using common infrastructure including sample holders, flow control, test sections, external radiant heaters, igniters, and diagnostics. SoFIE has been developed to meet the requirements of five unique investigations. It is currently being built and slated to begin operations on the ISS in July 2021. Given the general capabilities of the hardware insert, it is intended to be used as a facility for future researchers who can propose to NASA for related solid combustion studies.

Ferkul, Paul↗

Pterodactyl: 6-DOF Integration of Guidance and Control Algorithms in Genesis

The NASA-funded Pterodactyl project seeks to advance the current state-of-the-art for entry vehicles by developing novel guidance and control technologies for Deployable Entry Vehicles. This paper builds upon the Pterodactyl architecture that employed eight individually articulating flaps with two options for guidance, bank angle modulation with the Fully Numerical Predictor Corrector Entry Guidance (FNPEG) technique and angle of attack and sideslip modulation with FNPEG uncoupled range control. These, with a linear quadratic regulator controller, have previously been presented in separate 3-DOF trajectory simulations, one for translational motion and another for rotational dynamics. This work will show results for integrated 6-DOF simulations, leveraging recent advancements in trajectory simulation tools, namely the Julia-based Genesis package. Results show similar performance for both schemes when compared to the previously presented 3-DOF results once the controllers were adequately tuned. Angle of attack and sideslip modulation necessitated a controller designed at high dynamic pressure conditions, and bank angle modulation necessitated a controller designed at lower dynamic pressure conditions. Notably, angle of attack and sideslip modulation could not achieve the parachute deploy point target Mach number of 2, when only tuning the controller gains, thus future work will be needed to tune guidance specific parameters to achieve acceptable range targeting and guidance command tracking.

DEV↗

Pterodactyl: 6-DOF Integration of Guidance and Control Algorithms in Genesis

The NASA-funded Pterodactyl project seeks to advance the current state-of-the-art for entry vehicles by developing novel guidance and control technologies for Deployable Entry Vehicles. This paper builds upon the Pterodactyl architecture that employed eight individually articulating flaps with two options for guidance, bank angle modulation with the Fully Numerical Predictor Corrector Entry Guidance (FNPEG) technique and angle of attack and sideslip modulation with FNPEG uncoupled range control. These, with a linear quadratic regulator controller, have previously been presented in separate 3-DOF trajectory simulations, one for translational motion and another for rotational dynamics. This work will show results for integrated 6-DOF simulations, leveraging recent advancements in trajectory simulation tools, namely the Julia-based Genesis flight simulation tool. Results show similar performance for both schemes when compared to the previously presented 3-DOF results once the controllers were adequately tuned. Angle of attack and sideslip angle modulation necessitated a controller designed at high dynamic pressure conditions, and bank angle modulation necessitated a controller designed at lower dynamic pressure conditions. Notably, the angle of attack and sideslip modulation approach could not achieve the parachute deploy point target Mach number of 2, when only tuning the controller gains, so future work will be needed to improve guidance and control tracking.

EDL↗

Numerical calculations of velocity and pressure distribution around oscillating airfoils

An analytical procedure based on the Navier-Stokes equations was developed for analyzing and representing properties of unsteady viscous flow around oscillating obstacles. A variational formulation of the vorticity transport equation was discretized in finite element form and integrated numerically. At each time step of the numerical integration, the velocity field around the obstacle was determined for the instantaneous vorticity distribution from the finite element solution of Poisson's equation. The time-dependent boundary conditions around the oscillating obstacle were introduced as external constraints, using the Lagrangian Multiplier Technique, at each time step of the numerical integration. The procedure was then applied for determining pressures around obstacles oscillating in unsteady flow. The obtained results for a cylinder and an airfoil were illustrated in the form of streamlines and vorticity and pressure distributions.

Bratanow, T.↗

Encoding Y,I,Q component estimates of an NTSC composite signal

Several new techniques for reducing the entropy of a Y,I,Q component digital transmission are presented. A discussion of each technique is presented along with objective and subjective analysis. For a set of four test pictures, the average Y,I,Q entropy/pixel was reduced by nearly 80 percent. An initial component system entropy of 11.08 bits/pixel was reduced to 2.33 bits/pixel. Although exact numeric results are dependent on the separation technique used, the methods presented may be applied to systems using other component separation methods.

Border, David A.↗

Complete flowfields over low and wide angle AOTV conceptual configurations

Flowfields over two conceptual aeroassisted orbital transfer vehicles, a moderate lift-to-drag ratio biconic, and an axisymmetric zero lift aerobrake are investigated. The solution procedures employ a compressible Navier-Stokes time-asymptotic technique for blunt body flow and a Parabolized Navier-Stokes spatial marching technique for the primarily supersonic shock layer flow over low-angle bodies. Emphasis is placed on the laminar convective heating predictions and comparisons to experimental data where available. Code robustness and solution times are discussed and guidelines are given for the specification of various numerical parameters. The Parabolized Navier-Stokes technique has a high probability of solution success for the range of applications considered herein, and very good comparisons with experimental heat-transfer data have been obtained. The blunt body Navier-Stokes code is more sensitive to parameter specification for the aerobrake calculations. Still, coarse grid solutions can be obtained quickly to get a good first cut at flowfield definition. Fair to good comparisons have been obtained between experimental heating and pressure data and coarse and fine grid solutions.

Gnoffo, P. A.↗

Numerical solutions for laminar and turbulent viscous flow over single and multi-element airfoils using body-fitted coordinate systems

The technique of body-fitted coordinate systems is applied in numerical solutions of the complete time-dependent compressible and incompressible Navier-Stokes equations for laminar flow and to the time-dependent mean turbulent equations closed by modified Kolmogorov hypotheses for turbulent flow. Coordinate lines are automatically concentrated near to the bodies at higher Reynolds number so that accurate resolution of the large gradients near the solid boundaries is achieved. Two-dimensional bodies of arbitrary shapes are treated, the body contour(s) being simply input to the program. The complication of the body shape is thus removed from the problem.

Thompson, J. F.↗

Development of a boundary-layer-type solver based on simultaneous iteration technique for axisymmetric separated flows

A boundary-layer-type solver is developed for the numerical solution of axisymmetric separated flows. A new fully implicit coupling scheme for the viscous and inviscid regions is demonstrated. This fully implicit coupling technique is similar to the work of Carter, Veldman, and an extension of an earlier work of Halim and Hafez. A comparison is made for the convergence rate using this new fully implicit coupling technique and the semiimplicit coupling of Halim and Hafez. Numerical results using the fully implicit coupling are obtained for laminar incompressible separated flows, including a boattail and a series of trough geometries. Also, the near-wake flow problem is considered using the present formulation. A clear conclusion of this investigation is that the present scheme using the fully implicit coupling method converges at a faster rate than the semiimplicit coupling and the partially parabolized Navier-Stokes (PPNS) procedures.

Halim, A. A. M.↗

Algorithms for computing the geopotential using a simple density layer

Several algorithms have been developed for computing the potential and attraction of a simple density layer. These are numerical cubature, Taylor series, and a mixed analytic and numerical integration using a singularity-matching technique. A computer program has been written to combine these techniques for computing the disturbing acceleration on an artificial earth satellite. A total of 1640 equal-area, constant surface density blocks on an oblate spheroid are used. The singularity-matching algorithm is used in the subsatellite region, Taylor series in the surrounding zone, and numerical cubature on the rest of the earth.

Morrison, F.↗

Numerical simulation of a disk-shaped electron accelerating electrostatic probe

This paper presents a method utilizing the numerical plasma simulation particle-in-cell technique for studying electron accelerating probe characteristics. A probe configuration other than the commonly used plane, cylindrical, or spherical case will be studied; the probe configuration used is that of a disk. Cylindrical geometry is used in the simulation model. One of the advantages of a numerical simulation is the ease with which the details of the system can be examined. Therefore, in addition to the current-voltage characteristics of the static probe, other results are presented such as the time evolution of the plasma properties and studies of the sheath region surrounding the probe.

Nonnast, J. H.↗

Numerical study of three-dimensional separation and flow control at a wing/body junction

The problem of three-dimensional separation and flow control at a wing/body junction has been investigated numerically using a three-dimensional Navier-Stokes code. The numerical code employs an algebraic grid generation technique for generating the grid for unmodified junction and an elliptic grid generation technique for filleted fin junction. The results for laminar flow past a blunt fin/flat plate junction demonstrate that after grid refinement, the computations agree with experiment and reveal a strong dependency of the number of vortices at the junction on Mach number and Reynolds number. The numerical results for pressure distribution, particle paths and limiting streamlines for turbulent flow past a swept fin show a decrease in the peak pressure and in the extent of the separated flow region compared to the laminar case. The results for a filleted juncture indicate that the streamline patterns lose much of their vortical character with proper filleting. Fillets with a radius of three and one-half times the fin leading edge diameter or two times the incoming boundary layer thickness, significantly weaken the usual necklace interaction vortex for the Mach number and Reynolds number considered in the present study.

Ash, Robert L.↗

Numerical simulation of the control of the three-dimensional transition process in boundary layers

Surface heating techniques to control the three-dimensional laminar-turbulent transition process are numerically investigated for a water boundary layer. The Navier-Stokes and energy equations are solved using a fully implicit finite difference/spectral method. The spatially evolving boundary layer is simulated. Results of both passive and active methods of control are shown for small amplitude two-dimensional and three-dimensional disturbance waves. Control is also applied to the early stages of the secondary instability process using passive or active control techniques.

Kral, L. D.↗

Three-dimensional Euler solutions for long-duct nacelles

A three-dimensional Euler-equation computational technique has been developed to solve for the transonic flow past flow-through nacelles. The technique employs an approximately-factored alternating-direction implicit numerical algorithm and a radiation treatment of the outflow boundary. Studies are presented which show that the radiation treatment gives better numerical convergence than the condition of specifying the pressure at the outflow boundary. Calculations made with the technique are presented for a long-duct turbofan engine nacelle at a Mach number of 0.80 and angles of attack of 0 deg and 4 deg. Good agreement is shown between the computational results and wind-tunnel data. Problem areas are identified and recommendations are made for further numerical studies.

Compton, W. B., III↗

A two-dimensional finite difference solution for the transient thermal behavior of tubular solar collector

A numerical procedure was established using the finite-difference technique in the determination of the time-varying temperature distribution of a tubular solar collector under changing solar radiancy and ambient temperature. Three types of spatial discretization processes were considered and compared for their accuracy of computations and for selection of the shortest computer time and cost. The stability criteria of this technique was analyzed in detail to give the critical time increment to ensure stable computations. The results of the numerical analysis were in good agreement with the analytical solution previously reported. The numerical method proved to be a powerful tool in the investigation of the collector sensitivity to two different flow patterns and several flow control mechanisms.

Lansing, F. L.↗