Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical computation”

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,621 records · Page 90

Thermal Dispersion Within a Porous Medium Near a Solid Wall

The regenerator is a key component to Stirling cycle machine efficiency. Typical regenerators are of sintered fine wires or layers of fine-wire screens. Such porous materials are contained within solid-waH casings. Thermal energy exchange between the regenerator and the casing is important to cycle performance for the matrix and casing would not have the same axial temperature profile in an actual machine. Exchange from one to the other may allow shunting of thermal energy, reducing cycle efficiency. In this paper, temperature profiles within the near-wall region of the matrix are measured and thermal energy transport, termed thermal dispersion, is inferred. The data show how the wall affects thermal transport. Transport normal to the mean flow direction is by conduction within the solid and fluid and by advective transport within the matrix. In the near-wall region, both may be interrupted from their normal in-core pattern. Solid conduction paths are broken and scales of advective transport are damped. An equation is presented which describes this change for a wire screen mesh. The near-wall layer typically acts as an insulating layer. This should be considered in design or analysis. Effective thermal conductivity within the core is uniform. In-core transverse thermal effective conductivity values are compared to direct and indirect measurements reported elsewhere and to 3D numerical simulation results, computed previously and reported elsewhere. The 3-D CFD model is composed of six cylinders in cross flow, staggered in arrangement to match the dimensions and porosity of the matrix used in the experiments. The commercial code FLUENT is used to obtain the flow and thermal fields. The thermal dispersion and effective thermal conductivities for the matrix are computed from the results.

Simon, T.↗

SMART: The Future of Spaceflight Avionics

A novel avionics approach is necessary to meet the future needs of low cost space and lunar missions that require low mass and low power electronics. The current state of the art for avionics systems are centralized electronic units that perform the required spacecraft functions. These electronic units are usually custom-designed for each application and the approach compels avionics designers to have in-depth system knowledge before design can commence. The overall design, development, test and evaluation (DDT&E) cycle for this conventional approach requires long delivery times for space flight electronics and is very expensive. The Small Multi-purpose Advanced Reconfigurable Technology (SMART) concept is currently being developed to overcome the limitations of traditional avionics design. The SMART concept is based upon two multi-functional modules that can be reconfigured to drive and sense a variety of mechanical and electrical components. The SMART units are key to a distributed avionics architecture whereby the modules are located close to or right at the desired application point. The drive module, SMART-D, receives commands from the main computer and controls the spacecraft mechanisms and devices with localized feedback. The sensor module, SMART-S, is used to sense the environmental sensors and offload local limit checking from the main computer. There are numerous benefits that are realized by implementing the SMART system. Localized sensor signal conditioning electronics reduces signal loss and overall wiring mass. Localized drive electronics increase control bandwidth and minimize time lags for critical functions. These benefits in-turn reduce the main processor overhead functions. Since SMART units are standard flight qualified units, DDT&E is reduced and system design can commence much earlier in the design cycle. Increased production scale lowers individual piece part cost and using standard modules also reduces non-recurring costs. The benefit list continues, but the overall message is already evident: the SMART concept is an evolution in spacecraft avionics. SMART devices have the potential to change the design paradigm for future satellites, spacecraft and even commercial applications.

Alhorn, Dean C.↗

Validation of Model Forecasts of the Ambient Solar Wind

Independent and automated validation is a vital step in the progression of models from the research community into operational forecasting use. In this paper we describe a program in development at the CCMC to provide just such a comprehensive validation for models of the ambient solar wind in the inner heliosphere. We have built upon previous efforts published in the community, sharpened their definitions, and completed a baseline study. We also provide first results from this program of the comparative performance of the MHD models available at the CCMC against that of the Wang-Sheeley-Arge (WSA) model. An important goal of this effort is to provide a consistent validation to all available models. Clearly exposing the relative strengths and weaknesses of the different models will enable forecasters to craft more reliable ensemble forecasting strategies. Models of the ambient solar wind are developing rapidly as a result of improvements in data supply, numerical techniques, and computing resources. It is anticipated that in the next five to ten years, the MHD based models will supplant semi-empirical potential based models such as the WSA model, as the best available forecast models. We anticipate that this validation effort will track this evolution and so assist policy makers in gauging the value of past and future investment in modeling support.

Macneice, P. J.↗

Legacy of the Space Shuttle from an Aerodynamic and Aerothermodynamic Perspective

The development of the Space Shuttle Orbiter thermal protection system heating environment is described from a design stand point that began in the early 1970s. The desire for a light weight, reusable heat shield required the development of new technology, relative to previous manned spacecraft, and a systems approach to the design of the vehicle, entry guidance, and thermal protection system. Several unanticipated issues had to be resolved in both the entry and ascent phases of flight, which are discussed at a high level. During the life of the Program, significant improvements in computing power and numerical methods have been applied to Space Shuttle aerodynamic and aerothermodynamic issues, with the Shuttle Program often being the motivation, and or sponsor of the analysis development.

Martin, Fred W.↗

CO2 Insulation for Thermal Control of the Mars Science Laboratory

The National Aeronautics and Space Administration (NASA) is sending a large (>850 kg) rover as part of the Mars Science Laboratory (MSL) mission to Mars in 2011. The rover's primary power source is a Multi-Mission Radioisotope Thermoelectric Generator (MMRTG) that generates roughly 2000 W of heat, which is converted to approximately 110 W of electrical power for use by the rover electronics, science instruments, and mechanism-actuators. The large rover size and extreme thermal environments (cold and hot) for which the rover is designed for led to a sophisticated thermal control system to keep it within allowable temperature limits. The pre-existing Martian atmosphere of low thermal conductivity CO2 gas (8 Torr) is used to thermally protect the rover and its components from the extremely cold Martian environment (temperatures as low as -130 deg C). Conventional vacuum based insulation like Multi Layer Insulation (MLI) is not effective in a gaseous atmosphere, so engineered gaps between the warm rover internal components and the cold rover external structure were employed to implement this thermal isolation. Large gaps would lead to more thermal isolation, but would also require more of the precious volume available within the rover. Therefore, a balance of the degree of thermal isolation achieved vs. the volume of rover utilized is required to reach an acceptable design. The temperature differences between the controlled components and the rover structure vary from location to location so each gap has to be evaluated on a case-by-case basis to arrive at an optimal thickness. For every configuration and temperature difference, there is a critical thickness below which the heat transfer mechanism is dominated by simple gaseous thermal conduction. For larger gaps, the mechanism is dominated by natural convection. In general, convection leads to a poorer level of thermal isolation as compared to conduction. All these considerations play important roles in the optimization process. A three-step process was utilized to design this insulation. The first step is to come up with a simple, textbook based, closed-form equation assessment of gap thickness vs. resultant thermal isolation achieved. The second step is a more sophisticated numerical assessment using Computational Fluid Dynamics (CFD) software to investigate the effect of complicated geometries and temperature contours along them to arrive at the effective thermal isolation in a CO2 atmosphere. The third step is to test samples of representative geometries in a CO2 filled chamber to measure the thermal isolation achieved. The results of these assessments along with the consistency checks across these methods leads to the formulation of design-guidelines for gap implementation within the rover geometry. Finally, based on the geometric and functional constraints within the real rover system, a detailed design that accommodates all these factors is arrived at. This paper will describe in detail this entire process, the results of these assessments and the final design that was implemented.

CO2↗

Ultrasonic Measurement and Monitoring of Loads in Bolts used in Structural Joints

The paper is an overview of work by the author in measuring and monitoring loads in bolts using an ultrasonic extensometer. A number of cases of bolted joints are covered. These include, a clamped joint with clearance fit between the bolt and hole, a clamped joint with the bolt in an interference fit with the hole, a flanged joint which allows the flange and bolt to bend; and a shear joint in a clevis and tang configuration. These applications were initially developed for measuring and monitoring preload in the NASA Space Shuttle Orbiter critical joints but are also applicable for monitoring loads in other critical bolted joints of structures such as transportation bridges and other aerospace structures. The paper explains how to set‐up a model to estimate the load factor and accuracy for the ultrasonic preload application in a clamped joint with clearance fit. The ultrasonic preload application for clamped joint with bolt in an interference fit can also be used to measure diametrical interference between the bolt shank and hole; and interference pressure on the bolt shank. Model and experimental data are given to demonstrate use of ultrasonic measurements in a shear joint. A bolt in a flanged joint experiences both tensile and bending loads. This application involves measurement of bending and tensile preload in a bolt. The ultrasonic beam bends due to the bending load on the bolt. A numerical technique to compute the trace of ultrasonic ray is presented.

Koshti, Ajay↗

An Examination of Parameters Affecting Large Eddy Simulations of Flow Past a Square Cylinder

Separated flow over a bluff body is analyzed via large eddy simulations. The turbulent flow around a square cylinder features a variety of complex flow phenomena such as highly unsteady vortical structures, reverse flow in the near wall region, and wake turbulence. The formation of spanwise vortices is often times artificially suppressed in computations by either insufficient depth or a coarse spanwise resolution. As the resolution is refined and the domain extended, the artificial turbulent energy exchange between spanwise and streamwise turbulence is eliminated within the wake region. A parametric study is performed highlighting the effects of spanwise vortices where the spanwise computational domain's resolution and depth are varied. For Re=22,000, the mean and turbulent statistics computed from the numerical large eddy simulations (NLES) are in good agreement with experimental data. Von-Karman shedding is observed in the wake of the cylinder. Mesh independence is illustrated by comparing a mesh resolution of 2 million to 16 million. Sensitivities to time stepping were minimized and sampling frequency sensitivities were nonpresent. While increasing the spanwise depth and resolution can be costly, this practice was found to be necessary to eliminating the artificial turbulent energy exchange.

Mankbadi, M. R.↗

Effects of Distortion on Mass Flow Plug Calibration

A numerical, and experimental investigation to study the effects of flow distortion on a Mass Flow Plug (MFP) used to control and measure mass-flow during an inlet test has been conducted. The MFP was first calibrated using the WIND-US flow solver for uniform (undistorted) inflow conditions. These results are shown to compare favorably with an experimental calibration under similar conditions. The effects of distortion were investigated by imposing distorted flow conditions taken from an actual inlet test to the inflow plane of the numerical simulation. The computational fluid dynamic (CFD) based distortion study only showed the general trend in mass flow rate. The study used only total pressure as the upstream boundary condition, which was not enough to define the flow. A better simulation requires knowledge of the turbulence structure and a specific distortion pattern over a range of plug positions. It is recommended that future distortion studies utilize a rake with at least the same amount of pitot tubes as the AIP rake.

Inlets↗

Implementation of Kane's Method for a Spacecraft Composed of Multiple Rigid Bodies

Equations of motion are derived for a general spacecraft composed of rigid bodies connected via rotary (spherical or gimballed) joints in a tree topology. Several supporting concepts are developed in depth. Basis dyads aid in the transition from basis-free vector equations to component-wise equations. Joint partials allow abstraction of 1-DOF, 2-DOF, 3-DOF gimballed and spherical rotational joints to a common notation. The basic building block consisting of an "inner" body and an "outer" body connected by a joint enables efficient organization of arbitrary tree structures. Kane's equation is recast in a form which facilitates systematic assembly of large systems of equations, and exposes a relationship of Kane's equation to Newton and Euler's equations which is obscured by the usual presentation. The resulting system of dynamic equations is of minimum dimension, and is suitable for numerical solution by computer. Implementation is ·discussed, and illustrative simulation results are presented.

Stoneking, Eric T.↗

Structured Overlapping Grid Simulations of Contra-rotating Open Rotor Noise

Computational simulations using structured overlapping grids with the Launch Ascent and Vehicle Aerodynamics (LAVA) solver framework are presented for predicting tonal noise generated by a contra-rotating open rotor (CROR) propulsion system. A coupled Computational Fluid Dynamics (CFD) and Computational AeroAcoustics (CAA) numerical approach is applied. Three-dimensional time-accurate hybrid Reynolds Averaged Navier-Stokes/Large Eddy Simulation (RANS/LES) CFD simulations are performed in the inertial frame, including dynamic moving grids, using a higher-order accurate finite difference discretization on structured overlapping grids. A higher-order accurate free-stream preserving metric discretization with discrete enforcement of the Geometric Conservation Law (GCL) on moving curvilinear grids is used to create an accurate, efficient, and stable numerical scheme. The aeroacoustic analysis is based on a permeable surface Ffowcs Williams-Hawkings (FW-H) approach, evaluated in the frequency domain. A time-step sensitivity study was performed using only the forward row of blades to determine an adequate time-step. The numerical approach is validated against existing wind tunnel measurements.

Grid↗

An Analytical Method to Calculate Effective Elastic Properties of General Multifunctional Honeycomb Cores in Sandwich Composites

Sandwich composite structures are ideal configurations in which to incorporate additional functionality beyond load-carrying capabilities. The inner core-walls can be layered to incorporate other functions such as power storage for a battery. In this work we investigate an assemblage of analytical tools to compute effective properties that allow complex layered core architectures to be homogenized into a single continuum layer. This provides a great increase in computational efficiency to numerically simulate the structural response of multifunctional sandwich structures under applied loads. We present a coupled analytical method including an extensive numerical verification of the accuracy of this method.

Sather, Erik↗

Estimation of Effective Elastic Properties of General Multifunctional Honeycomb Structures Using a Unit Cell Method

Sandwich composite structures are ideal configurations in which to incorporate additional functionality beyond load carrying capabilities. The inner core can be layered to facilitate other functions such as power storage for a battery. In this work we investigate an assemblage of analytical tools to compute effective properties that allow complex layered core architectures to be homogenized into a single continuum layer. This provides a great increase in computational efficiency to numerically simulate the structural response of multifunctional sandwich structures under applied loads.

Krishnamurthy, T.↗

Planetary-Induced Nutation of the Earth - Direct Terms

Small torques by the planets can alter the direction of the Earth's rotation axis in space. These direct planetary-induced nutations have been computed using a numerical technique. 55 nutation terms with periods between 0.2 and 243 yr are presented. The maximum nutation is 0.5 milliarcseconds. Also given are planetary-induced rates and accelerations of precession and obliquity.

Planetary-Induced↗

Numerical investigations of background pressure effects and channel erosion in the SPT-140 Hall thruster for the Psyche mission

Four SPT-140 Hall thrusters are the main components of the electric propulsion system for the Psyche mission. We employ the 2-D (r-z) code Hall2De to conduct numerical simulations of the SPT-140 with the goal of estimating erosion of the acceleration channel rings for the Psyche mission profile. The simulations are informed by direct measurements of the plasma conditions in the acceleration channel that were obtained using laser-induced fluorescence. We validate our simulation results with additional plasma measurements, wear test erosion rates, and performance data as a function of background pressure. We then make use of a series of numerical simulations to compute the progressive erosion of the channel walls, using as input parameters a profile for the Psyche mission that includes operation at 4.5, 2.5 and 1.0 kW with duration of all segments increased by a 50% margin. We find that the channel walls are not completely eroded at the end of the mission (plus 50% margin). We also observe a progressive decrease in thrust as the channel walls erode that is consistent with the wear test and flight data, with a rate of ~4.5 mN/kh for the first 2,500 h and constant thrust after that. The reasons for these thrust trends are explained.

Lenguito, Giovanni↗

Extension and Validation of the Pendulum Model for Longitudinal Solar Prominence Oscillations

Context. Longitudinal oscillations in prominences are common phenomena on the Sun. These oscillations can be used to infer the geometry and intensity of the filament magnetic field. Previous theoretical studies of longitudinal oscillations made two simplifying assumptions: uniform gravity and semicircular dips on the supporting flux tubes. However, the gravity is not uniform and realistic dips are not semicircular. Aims. Our aim is to understand the effects of including the nonuniform solar gravity on longitudinal oscillations and explore the validity of the pendulum model with different flux-tube geometries. Methods. We first derived the equation describing the motion of the plasma along the flux tube including the effects of nonuniform gravity, yielding corrections to the original pendulum model. We also computed the full numerical solutions for the normal modes and compared them with the new pendulum approximation. Results. We find that the nonuniform gravity introduces a significant modification in the pendulum model. We also found a cut-off period; i.e., the longitudinal oscillations cannot have a period longer than 167 min. In addition, considering different tube geometries, the period depends almost exclusively on the radius of curvature at the bottom of the dip. Conclusions. We conclude that nonuniform gravity significantly modifies the pendulum model. These corrections are important for prominence seismology, because the inferred values of the radius of curvature and minimum magnetic-field strength differ substantially from those of the old model. However, we find that the corrected pendulum model is quite robust and is still valid for noncircular dips.

Sun: corona↗

Salt Distribution from Freezing Intrusions in Ice Shells on Ocean Worlds: Application to Europa

Several icy moons and dwarf planets appear to have hosted subsurface liquid water. Liquid water intruding upwards into the icy outer shells of these worlds freezes, forming ice and (from ocean solutes) non-ice solids. Here, we model concentrated aqueous solutions below 273 K to simulate the compositional evolution of freezing spherical intrusions. Starting solutions are based on five previously reported compositional end members for Europa's ocean. For moderate-pH end members dominated by chloride, sulfate, and/or carbonate, the solids formed include Ca-, Mg-, and Na-sulfates and -carbonates, as well as Na- and K-chlorides. For silica-rich, high-pH end members, abundant amorphous silica forms with, potentially, similarly abundant NaOH and KOH. We further develop a new numerical model to compute the spatial distribution of the formed solids and residual brine as freezing progresses. If non-ice solids settle to the bottom, their deposits tend to have stacked hourglass shapes, widening each time the crystallization temperature of a new solid is reached. We discuss the applicability of this model to vertical fractures and global freezing of a subsurface ocean. These results inform (i) how compositional heterogeneities may affect the thermophysical properties of ice shells, which in turn influence convective and cryovolcanic transport, (ii) the compatibility of brine pockets with physicochemical conditions suitable for microbial life, and (iii) possible measurements of compositional heterogeneities within ice shells by spacecraft such as NASA's Europa Clipper and ESA's JUICE missions. The methodology developed here is applicable to other ice-covered ocean worlds.

Planetary structure↗

Development of an Extended Reality (XR) Tool for Earth Science Visualization

In this presentation, we will discuss our work in adapting the NASA open source XR software, the Mixed Reality Exploration Toolkit (MRET), to an earth science domain. MRET is a NASA open source XR software for rapidly building extended reality (XR) environments for NASA domain problems, e.g., pulling in CAD models of thermal vac chamber and Roman Space Telescope to do fit checks. Primarily used for hardware integration & test, we have been adapting and extending MRET for science problems. Traditionally, scientists view and analyze the result of calculated or measured observables with static 1-D, 2-D or 3-D plots. It can be difficult to identify, track and understand the evolution of key features due to poor viewing angles and the nature of flat computer screens. Additionally, numerical models, such as the NASA GEOS climate model, are almost exclusively formulated and analyzed on Eulerian grids with points fixed in space and time. However, atmospheric phenomena such as convective clouds, hurricanes and wildfire smoke plumes move with the 3-D flow field, and it is often difficult and unnatural to understand these phenomena in an Eulerian reference frame as opposed to the Lagrangian reference frame in which nature operates. As part of an Earth Science Technology Office (ESTO) proposal, we have been adapting MRET to be a scientific exploration and analysis XR tool with integrated Lagrangian Dynamics (LD) for the Goddard Earth Observing System (GEOS) numerical weather prediction model. We believe this will help scientists identify, track, and understand the evolution of Earth Science phenomena. This presentation will discuss current results in our work in developing this XR tool for Earth Science.

Thomas G Grubb↗

Phased Microphone Array on Aircraft Fuselage

This paper studies some challenging features for phased microphone arrays on aircraft fuselage, including the effects of sound scattering such as refraction by nonuniform boundary layer flows, reflection by curved surfaces, and diffraction by sharp edges and smooth geometry. The objective is to show if and how source locations and source amplitudes can be accurately captured by the array data analysis, in the presence of these features. The study starts with canonical problems to separately examine the individual features, making use of their respective analytical solutions, which provide exact results for all quantities needed for the array data processing. To simulate realistic flight environments, a generic aircraft geometry is then utilized, involving multiple scattering mechanisms and multiple noise sources to model aircraft engine noise. For this application, numerical solutions are computed to provide noise propagation and scattering from the sources to the microphones, as well as the quantities needed for array data analysis. The study reveals the effects of various features on fuselage-mounted arrays and provides guidelines to account for these effects in array data analysis.

Yueping Guo↗