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

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

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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

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↗

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

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

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↗

Trend Detection of Atmospheric Time Series: Incorporating Appropriate Uncertainty Estimates and Handling Extreme Events

This paper is aimed at atmospheric scientists without formal training in statistical theory. Its goal is to, 1) provide a critical review of the rationale for trend analysis of the time series typically encountered in the field of atmospheric chemistry; 2) describe a range of trend-detection methods; and 3) demonstrate effective means of conveying the results to a general audience. Trend detections in atmospheric chemical composition data are often challenged by a variety of sources of uncertainty, which often behave differently to other environmental phenomena such as temperature, precipitation rate, or stream flow, and may require specific methods depending on the science questions to be addressed. Some sources of uncertainty can be explicitly included in the model specification, such as autocorrelation and seasonality, but some inherent uncertainties are difficult to quantify, such as data heterogeneity and measurement uncertainty due to the combined effect of short- and long-term natural variability, instrumental stability, and aggregation of data from sparse sampling frequency. Failure to account for these uncertainties might result in an inappropriate inference of the trends and their estimation errors. On the other hand, the variation in extreme events might be interesting for different scientific questions, for example, the frequency of extremely high surface ozone events and their relevance to human health. In this study we aim to, 1) review trend detection methods for addressing different levels of data complexity in different chemical species; 2) demonstrate that the incorporation of scientifically interpretable covariates can outperform pure numerical curve fitting techniques in terms of uncertainty reduction and improved predictability; 3) illustrate the study of trends based on extreme quantiles that can provide insight beyond standard mean or median based trend estimates; and 4) present an advanced method of quantifying regional trends based on the inter-site correlations of multi-site data. All demonstrations are based on time series of observed trace gases relevant to atmospheric chemistry, but the methods can be applied to other environmental data sets.

Trace gas↗

The CARMENES Search for Exoplanets Around M Dwarfs: Revisiting the GJ 581 Multi-Planetary System With New Doppler Measurements From CARMENES, HARPS, and HIRES

Context.GJ 581 is a nearby M dwarf known to host a packed multiple planet system composed of two super-Earths and a Neptune-mass planet. We present new orbital analyses of the GJ 581 system, utilizing recent radial velocity (RV) data obtained from the CARMENES spectrograph combined with newly reprocessed archival data from the HARPS and HIRES spectrographs. Aims.Our aim was to analyze the post-discovery spectroscopic data of GJ 581, which were obtained with CARMENES. In addition, we used publicly available HIRES and HARPS spectroscopic data to seek evidence of the known and disputed exoplanets in this system. We aimed to investigate the stellar activity of GJ 581 and update the planetary system’s orbital parameters using state-of-the-art numerical models and techniques. Methods.We performed a periodogram analysis of the available precise CARMENES, HIRES, and HARPS RVs and of stellar activity indicators. We conducted detailed orbital analyses by testing various orbital configurations consistent with the RV data. We studied the posterior probability distribution of the parameters fit to the data and we explored the long-term stability and overall orbital dynamics of the GJ 581 system. Results. We refined the orbital parameters of the GJ 581 system using the most precise and complete set of Doppler data available. Consistent with the existing literature, our analysis confirms that the system is unequivocally composed of only three planets detectable in the present data, dismissing the putative planet GJ 581 d as an artifact of stellar activity. Our N-body fit reveals that the system’s inclination is i=47.0+14.6−13.0deg, which implies that the planets could be up to 30% more massive than their previously reported minimum masses. Furthermore, we report that the GJ 581 system exhibits long-term stability, as indicated by the posterior probability distribution, characterized by secular dynamical interactions without the involvement of mean motion resonances.

Satellites↗

Turbulence Modeling near Gas-Liquid Interfaces for Cryogenic Fluid Management

In cryogenic fluid management (CFM), gas-liquid interface turbulence transport (ITT) heavily influences crucial design parameters of spaceflight missions. ITT and phase change processes require length-scales and time-scales that are around 4-8 orders of magnitude smaller than the scales of CFM operations in spaceflights. A sharp interface method based parsimonious yet elegant turbulence modeling approach was developed from experimental observations and high-fidelity numerical studies. The technique was validated using flight-scale CFM tests.

cryogenic fluid management↗

Viscous shock-layer analysis of hypersonic flows over long slender bodies

Numerical solutions from the time-steady viscous shock-layer equations are presented for the hypersonic laminar and turbulent flow of a perfect gas over long slender bodies. These results are obtained from a spatial-marching implicit finite-difference technique. Detailed comparisons have been made with other predictions and experimental data to assess the accuracy of the present numerical technique, especially for slender-body flows. Results from the present method show that coupling the normal momentum and continuity equations and the use of the Vigneron pressure condition in the subsonic nose region give quite accurate and stable results. These results compare (with some exceptions) favorably with those obtained using PNS and other VSL methods.

Gupta, R. N.↗