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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,171 records · Page 65

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Flow and acoustic field due to an inclined plate with a downstream splitter

In the present work, the high Reynolds number flow past an inclined plate with a splitter plate placed in its wake is considered numerically. A numerical conformal mapping technique is employed to transform the two-plate system into the same number of cylinders: the flow field is assumed to be two-dimensional. The vortex shedding from the inclined plate is modelled using the discrete vortex method. It is shown that the splitter plate has a profound effect on the development of the flow over a range of values of a suitably defined offset parameter and for a range of positions of the leading edge of the splitter plate. The acoustic field is also calculated and the spectrum reflects the flow results.

DISCRETE VORTEX METHOD↗

NASA NTRS1996

Artificial Boundary Conditions for Computation of Oscillating External Flows

In this paper, we propose a new technique for the numerical treatment of external flow problems with oscillatory behavior of the solution in time. Specifically, we consider the case of unbounded compressible viscous plane flow past a finite body (airfoil). Oscillations of the flow in time may be caused by the time-periodic injection of fluid into the boundary layer, which in accordance with experimental data, may essentially increase the performance of the airfoil. To conduct the actual computations, we have to somehow restrict the original unbounded domain, that is, to introduce an artificial (external) boundary and to further consider only a finite computational domain. Consequently, we will need to formulate some artificial boundary conditions (ABC's) at the introduced external boundary. The ABC's we are aiming to obtain must meet a fundamental requirement. One should be able to uniquely complement the solution calculated inside the finite computational domain to its infinite exterior so that the original problem is solved within the desired accuracy. Our construction of such ABC's for oscillating flows is based on an essential assumption: the Navier-Stokes equations can be linearized in the far field against the free-stream back- ground. To actually compute the ABC's, we represent the far-field solution as a Fourier series in time and then apply the Difference Potentials Method (DPM) of V. S. Ryaben'kii. This paper contains a general theoretical description of the algorithm for setting the DPM-based ABC's for time-periodic external flows. Based on our experience in implementing analogous ABC's for steady-state problems (a simpler case), we expect that these boundary conditions will become an effective tool for constructing robust numerical methods to calculate oscillatory flows.

Tsynkov, S. V.↗

NASA NTRS1996

Helicopter Rotor Blade Computation in Unsteady Flows Using Moving Overset Grids

An overset grid thin-layer Navier-Stokes code has been extended to include dynamic motion of helicopter rotor blades through relative grid motion. The unsteady flowfield and airloads on an AH-IG rotor in forward flight were computed to verify the methodology and to demonstrate the method's potential usefulness towards comprehensive helicopter codes. In addition, the method uses the blade's first harmonics measured in the flight test to prescribe the blade motion. The solution was impulsively started and became periodic in less than three rotor revolutions. Detailed unsteady numerical flow visualization techniques were applied to the entire unsteady data set of five rotor revolutions and exhibited flowfield features such as blade vortex interaction and wake roll-up. The unsteady blade loads and surface pressures compare well against those from flight measurements. Details of the method, a discussion of the resulting predicted flowfield, and requirements for future work are presented. Overall, given the proper blade dynamics, this method can compute the unsteady flowfield of a general helicopter rotor in forward flight.

Ahmad, Jasim↗

NASA NTRS2000

Control of Convection by Dynamic Magnetic Fields for VB, FZ and THM Crystal Growth Application

The Traveling Magnetic Field (TMF) technique, recently proposed to control meridional flow in electrically conducting melts is reviewed. In particular, the natural convection damping capability of this technique has been numerically demonstrated with the implication of significantly improving crystal quality. Advantages of the traveling magnetic field, in comparison to the more mature rotating magnetic field method, are discussed. Finally, a detailed one-dimensional planar TMF model has been developed and is presented.

Mazuruk, Konstantin↗

NASA NTRS2000

Resin Film Infusion (RFI) Process Modeling for Large Transport Aircraft Wing Structures

This investigation completed the verification of a three-dimensional resin transfer molding/resin film infusion (RTM/RFI) process simulation model. The model incorporates resin flow through an anisotropic carbon fiber preform, cure kinetics of the resin, and heat transfer within the preform/tool assembly. The computer model can predict the flow front location, resin pressure distribution, and thermal profiles in the modeled part. The formulation for the flow model is given using the finite element/control volume (FE/CV) technique based on Darcy's Law of creeping flow through a porous media. The FE/CV technique is a numerically efficient method for finding the flow front location and the fluid pressure. The heat transfer model is based on the three-dimensional, transient heat conduction equation, including heat generation. Boundary conditions include specified temperature and convection. The code was designed with a modular approach so the flow and/or the thermal module may be turned on or off as desired. Both models are solved sequentially in a quasi-steady state fashion. A mesh refinement study was completed on a one-element thick model to determine the recommended size of elements that would result in a converged model for a typical RFI analysis. Guidelines are established for checking the convergence of a model, and the recommended element sizes are listed. Several experiments were conducted and computer simulations of the experiments were run to verify the simulation model. Isothermal, non-reacting flow in a T-stiffened section was simulated to verify the flow module. Predicted infiltration times were within 12-20% of measured times. The predicted pressures were approximately 50% of the measured pressures. A study was performed to attempt to explain the difference in pressures. Non-isothermal experiments with a reactive resin were modeled to verify the thermal module and the resin model. Two panels were manufactured using the RFI process. One was a stepped panel and the other was a panel with two 'T' stiffeners. The difference between the predicted infiltration times and the experimental times was 4% to 23%.

Loos, Alfred C.↗

NASA NTRS2000

Convection Induced by Traveling Magnetic Fields in Semiconductor Melts

Axisymmetric traveling magnetic fields (TMF) can be beneficial for crystal growth applications. such as the vertical Bridgman, float zone or traveling heater methods. TMF induces a basic flow in the form of a single roll. This type of flow can enhance mass and heat transfer to the growing crystal. More importantly, the TMF Lorentz body force induced in the system can counterbalance the buoyancy forces, so the resulting convection can be much smaller and even the direction of it can be changed. In this presentation, we display basic features of this novel technique. In particular, numerical calculations of the Lorentz force for arbitrary frequencies will be presented along with induced steady-state fluid flow profiles. Also, numerical modeling of the TMF counter-balancing natural convection in vertical Bridgman systems will be demonstrated.

Konstantin, Mazuruk↗

NASA NTRS2001

A Leo Satellite Navigation Algorithm Based on GPS and Magnetometer Data

The Global Positioning System (GPS) has become a standard method for low cost onboard satellite orbit determination. The use of a GPS receiver as an attitude and rate sensor has also been developed in the recent past. Additionally, focus has been given to attitude and orbit estimation using the magnetometer, a low cost, reliable sensor. Combining measurements from both GPS and a magnetometer can provide a robust navigation system that takes advantage of the estimation qualities of both measurements. Ultimately, a low cost, accurate navigation system can result, potentially eliminating the need for more costly sensors, including gyroscopes. This work presents the development of a technique to eliminate numerical differentiation of the GPS phase measurements and also compares the use of one versus two GPS satellites.

Deutschmann, Julie↗

NASA NTRS2002

Magnetic Field Applications in Semiconductor Crystal Growth and Metallurgy

The Traveling Magnetic Field (TMF) technique, recently proposed to control meridional flow in electrically conducting melts, is reviewed. In particular, the natural convection damping capability of this technique has been numerically demonstrated with the implication of significantly improving crystal quality. Advantages of the traveling magnetic field, in comparison to the more mature rotating magnetic field method, are discussed. Finally, results of experiments with mixing metallic alloys in long ampoules using TMF is presented

Mazuruk, Konstantin↗

NASA NTRS2003

Formation Flying Design and Applications in Weak Stability Boundary Regions

Weak Stability regions serve as superior locations for interferometric scientific investigations. These regions are often selected to minimize environmental disturbances and maximize observing efficiency. Design of formations in these regions are becoming ever more challenging as more complex missions are envisioned. The development of algorithms to enable the capability for formation design must be further enabled to incorporate better understanding of WSB solution space. This development will improve the efficiency and expand the capabilities of current approaches. The Goddard Space Flight Center (GSFC) is currently supporting multiple formation missions in WSB regions. This end-to-end support consists of mission operations, trajectory design, and control. It also includes both algorithm and software development. The Constellation-X, Maxim, and Stellar Imager missions are examples of the use of improved numerical methods for attaining constrained formation geometries and controlling their dynamical evolution. This paper presents a survey of formation missions in the WSB regions and a brief description of the formation design using numerical and dynamical techniques.

Folta, David↗

NASA NTRS2004

Computational Simulation of Acoustic Modes in Rocket Combustors

A combination of computational fluid dynamic analysis and analytical solutions is being used to characterize the dominant modes in liquid rocket engines in conjunction with laboratory experiments. The analytical solutions are based on simplified geometries and flow conditions and are used for careful validation of the numerical formulation. The validated computational model is then extended to realistic geometries and flow conditions to test the effects of various parameters on chamber modes, to guide and interpret companion laboratory experiments in simplified combustors, and to scale the measurements to engine operating conditions. In turn, the experiments are used to validate and improve the model. The present paper gives an overview of the numerical and analytical techniques along with comparisons illustrating the accuracy of the computations as a function of grid resolution. A representative parametric study of the effect of combustor mean flow Mach number and combustor aspect ratio on the chamber modes is then presented for both transverse and longitudinal modes. The results show that higher mean flow Mach numbers drive the modes to lower frequencies. Estimates of transverse wave mechanics in a high aspect ratio combustor are then contrasted with longitudinal modes in a long and narrow combustor to provide understanding of potential experimental simulations.

Harper, Brent↗

NASA NTRS2006

Towards an Airframe Noise Prediction Methodology: Survey of Current Approaches

In this paper, we present a critical survey of the current airframe noise (AFN) prediction methodologies. Four methodologies are recognized. These are the fully analytic method, CFD combined with the acoustic analogy, the semi-empirical method and fully numerical method. It is argued that for the immediate need of the aircraft industry, the semi-empirical method based on recent high quality acoustic database is the best available method. The method based on CFD and the Ffowcs William- Hawkings (FW-H) equation with penetrable data surface (FW-Hpds ) has advanced considerably and much experience has been gained in its use. However, more research is needed in the near future particularly in the area of turbulence simulation. The fully numerical method will take longer to reach maturity. Based on the current trends, it is predicted that this method will eventually develop into the method of choice. Both the turbulence simulation and propagation methods need to develop more for this method to become useful. Nonetheless, the authors propose that the method based on a combination of numerical and analytical techniques, e.g., CFD combined with FW-H equation, should also be worked on. In this effort, the current symbolic algebra software will allow more analytical approaches to be incorporated into AFN prediction methods.

Farassat, Fereidoun↗

NASA NTRS2007

GLAST Large Area Telescope - Design and Optimization of Anti-Coincidence Detector

The gamma-ray Large Area Telescope (GLAST) with LAT as its main instrument is now in the final stage of integration into the spacecraft and is scheduled for launch in November 2007. The Anti-Coincidence Detector (ACD) is one of the key telescope systems, providing much of the charged particle background rejection. Taking into account the required large ACD detecting area (more than 8 sq m) and high efficiency of relativistic single charged particle detection (greater than 0.9997 over the entire ACD area), the design was very challenging. It included the choice of conceptual design, extensive Monte Carlo simulations and trade-offs of the experimental techniques, supported by numerous bench tests, custom electronics, careful mechanical design and extensive mechanical/environmental tests, precision fabrication and finally the decisive step- final integration and system-level tests. Currently the LAT (with ACD integrated) is integrated into the spacecraft and undergoing the numerous spacecraft-level tests

Moiseev, Alexander↗

NASA NTRS2007

Ultra-High Temperature Materials Characterization for Propulsion Applications

Propulsion system efficiency increases as operating temperatures are increased. Some very high-temperature materials are being developed, including refractory metal alloys, carbides, borides, and silicides. System design requires data for materials properties at operating temperatures. Materials property data are not available for many materials of interest at the desired operating temperatures (up to approx. 3000 K). The objective of this work is to provide important physical property data at ultra-high temperatures. The MSFC Electrostatic levitation (ESL) facility can provide measurements of thermophysical properties which include: creep strength, density and thermal expansion for materials being developed for propulsion applications. The ESL facility uses electrostatic fields to position samples between electrodes during processing and characterization studies. Because the samples float between the electrodes during studies, they are free from any contact with a container or test apparatus. This provides a high purity environment for the study of high-temperature, reactive materials. ESL can be used to process a wide variety of materials including metals, alloys, ceramics, glasses and semiconductors. The MSFC ESL has provided non-contact measurements of properties of materials up to 3400 C. Density and thermal expansion are measured by analyzing digital images of the sample at different temperatures. Our novel, non-contact method for measuring creep uses rapid rotation to deform the sample. Digital images of the deformed samples are analyzed to obtain the creep properties, which match those obtained using ASTM Standard E-139 for Nb at 1985 C. Data from selected ESL-based characterization studies will be presented. The ESL technique could support numerous propulsion technologies by advancing the knowledge base and the technology readiness level for ultra-high temperature materials. Applications include non-eroding nozzle materials and lightweight, high-temperature alloys for turbines and structures.

Rogers, Jan↗

NASA NTRS2010

Comprehensive T-Matrix Reference Database: A 2007-2009 Update

The T-matrix method is among the most versatile, efficient, and widely used theoretical techniques for the numerically exact computation of electromagnetic scattering by homogeneous and composite particles, clusters of particles, discrete random media, and particles in the vicinity of an interface separating two half-spaces with different refractive indices. This paper presents an update to the comprehensive database of T-matrix publications compiled by us previously and includes the publications that appeared since 2007. It also lists several earlier publications not included in the original database.

Mishchenko, Michael I.↗

NASA NTRS2012

A Review of Low Frequency Electromagnetic Wave Phenomena Related to Tropospheric-Ionospheric Coupling Mechanisms

Investigation of coupling mechanisms between the troposphere and the ionosphere requires a multidisciplinary approach involving several branches of atmospheric sciences, from meteorology, atmospheric chemistry, and fulminology to aeronomy, plasma physics, and space weather. In this work, we review low frequency electromagnetic wave propagation in the Earth-ionosphere cavity from a troposphere-ionosphere coupling perspective. We discuss electromagnetic wave generation, propagation, and resonance phenomena, considering atmospheric, ionospheric and magnetospheric sources, from lightning and transient luminous events at low altitude to Alfven waves and particle precipitation related to solar and magnetospheric processes. We review in situ ionospheric processes as well as surface and space weather phenomena that drive troposphere-ionosphere dynamics. Effects of aerosols, water vapor distribution, thermodynamic parameters, and cloud charge separation and electrification processes on atmospheric electricity and electromagnetic waves are reviewed. We also briefly revisit ionospheric irregularities such as spread-F and explosive spread-F, sporadic-E, traveling ionospheric disturbances, Trimpi effect, and hiss and plasma turbulence. Regarding the role of the lower boundary of the cavity, we review transient surface phenomena, including seismic activity, earthquakes, volcanic processes and dust electrification. The role of surface and atmospheric gravity waves in ionospheric dynamics is also briefly addressed. We summarize analytical and numerical tools and techniques to model low frequency electromagnetic wave propagation and solving inverse problems and summarize in a final section a few challenging subjects that are important for a better understanding of tropospheric-ionospheric coupling mechanisms.

Simoes, Fernando↗

NASA NTRS2016

Infrared Observations of the Orion Capsule During EFT-1 Hypersonic Reentry

High-resolution infrared observations of the Orion capsule during its atmospheric reentry on December 5, 2015 were made from a US Navy NP-3D. This aircraft, equipped with a long-range optical sensor system, tracked the capsule from Mach 10 to 7 from a distance of approximately 60 nmi. Global surface temperatures of the capsule's thermal heatshield were derived from near infrared intensity measurements. The global surface temperature measurements complemented onboard instrumentation and were invaluable to the interpretation of the in-depth thermocouple measurements which rely on inverse heat transfer methods and material response codes to infer the desired surface temperature from the sub-surface measurements. The full paper will address the motivations behind the NASA Engineering Safety Center sponsored observation and highlight premission planning processes with an emphasis on aircraft placement, optimal instrument configuration and sensor calibrations. Critical aspects of mission operations coordinated from the NASA Johnson Spaceflight Center and integration with the JSC Flight Test Management Office will be discussed. A summary of the imagery that was obtained and processed to global surface temperature will be presented. At the capsule's point of closest approach relative to the imaging system, the spatial resolution was estimated to be approximately 15-inches per pixel and was sufficient to identify localized temperature increases associated with compression pad support hardware on the heatshield. The full paper will discuss the synergy of the quantitative imagery derived temperature maps with in-situ thermocouple measurements. Comparison of limited onboard surface thermocouple data to the image derived surface temperature will be presented. The two complimentary measurements serve as an example of the effective leveraging of resources to advance the understanding of high Mach number environments associated with an ablated heatshield and provide unique data for the validation of design tools and numerical flight simulation techniques. Collaborative opportunities and technology investments in support of planned observations of NASA's next Orion flight test in 2018 will be explored in the full manuscript.

Horvath, Thomas J.↗

NASA NTRS2017

Geometrically Flexible and Efficient Flow Analysis of High Speed Vehicles Via Domain Decomposition, Part 1: Unstructured-Grid Solver for High Speed Flows

The ability to solve the equations governing the hypersonic turbulent flow of a real gas on unstructured grids using a spatially-elliptic, 2nd-order accurate, cell-centered, finite-volume method has been recently implemented in the VULCAN-CFD code. This paper describes the key numerical methods and techniques that were found to be required to robustly obtain accurate solutions to hypersonic flows on non-hex-dominant unstructured grids. The methods and techniques described include: an augmented stencil, weighted linear least squares, cell-average gradient method, a robust multidimensional cell-average gradient-limiter process that is consistent with the augmented stencil of the cell-average gradient method and a cell-face gradient method that contains a cell skewness sensitive damping term derived using hyperbolic diffusion based concepts. A data-parallel matrix-based symmetric Gauss-Seidel point-implicit scheme, used to solve the governing equations, is described and shown to be more robust and efficient than a matrix-free alternative. In addition, a y+ adaptive turbulent wall boundary condition methodology is presented. This boundary condition methodology is deigned to automatically switch between a solve-to-the-wall and a wall-matching-function boundary condition based on the local y+ of the 1st cell center off the wall. The aforementioned methods and techniques are then applied to a series of hypersonic and supersonic turbulent flat plate unit tests to examine the efficiency, robustness and convergence behavior of the implicit scheme and to determine the ability of the solve-to-the-wall and y+ adaptive turbulent wall boundary conditions to reproduce the turbulent law-of-the-wall. Finally, the thermally perfect, chemically frozen, Mach 7.8 turbulent flow of air through a scramjet flow-path is computed and compared with experimental data to demonstrate the robustness, accuracy and convergence behavior of the unstructured-grid solver for a realistic 3-D geometry on a non-hex-dominant grid.

White, Jeffery A.↗

NASA NTRS2018

Multi-Scale Modeling of Low-Density Carbon-Phenolic Ablators

Protecting a spacecraft during atmospheric entry is one of highest risk factors that needs to be mitigated during design of a space exploration mission. At entry speeds from space, air turns into high-temperature plasma, and spacecraft Thermal Protection Systems (TPS) are needed to protect the vehicle payload. Modern successful material architectures of spacecraft shields use a porous carbon fiber substrate impregnated with phenolic as an ablator material. In the lecture, efforts to build a Predictive Material Modeling framework for porous ablators from micro-scale to macro-scale will be presented. Several numerical methods and techniques will be summarized that use voxelized images to compute geometrical properties of the porous substrate. These computed properties include porosity, specific surface area and tortuosity that are otherwise indirectly measured through experimental techniques. Direct simulation Monte Carlo (DSMC), a particle-based method for approximating the Boltzmann equation, is used to compute the permeability coefficient of the porous substrate based on its digitized representation. The method computes the flow within the microstructure, where the size of the pores may approach the mean-free-path of the flow. Finally, a high-fidelity model implemented in PATO (Porous-material Analysis Toolbox) is discussed, and some examples of ablative material response are presented, including for the first time 3D simulations of the full tiled heat shield for the Mars Science Laboratory (MSL) capsule.

Ablators↗