Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solve time”

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 973 records · Page 54

Computational Study of Flow Establishment in a Ram Accelerator

The temporal evolution of the combustion process established during projectile transition from the launch tube into the ram accelerator section containing an explosive hydrogen-oxygen-argon gas mixture is studied. The Navier-Stokes equations for chemically reacting flow are solved in a fully coupled manner, using an implicit, time accurate algorithm. The solution procedure is based on a spatially second order total variation diminishing scheme and a temporally second order, variable-step, backward differentiation formula method. The hydrogen-oxygen chemistry is modeled with a 9-species, 19-step mechanism. The accuracy of the solution method is first demonstrated by several benchmark calculations. Numerical simulations of two ram accelerator configurations are then presented. In particular, the temporal developments of shock-induced combustion and thrust forces are followed. Positive thrust is established in both cases; however, in one of the ram accelerator configurations studied, combustion in the boundary layer enhances its separation, ultimately resulting in unstart.

Yungster, S.↗

Nozzle Flow with Vibrational Nonequilibrium

Flow of nitrogen gas through a converging-diverging nozzle is simulated. The flow is modeled using the Navier-Stokes equations that have been modified for vibrational nonequilibrium. The energy equation is replaced by two equations. One equation accounts for energy effects due to the translational and rotational degrees of freedom, and the other accounts for the affects due to the vibrational degree of freedom. The energy equations are coupled by a relaxation time which measures the time required for the vibrational energy component to equilibrate with the translational and rotational energy components. An improved relaxation time is used in this thesis. The equations are solved numerically using the Steger-Warming flux vector splitting method and the Implicit MacCormack method. The results show that uniform flow is produced outside of the boundary layer. Nonequilibrium exists in both the converging and diverging nozzle sections. The boundary layer region is characterized by a marked increase in translational-rotational temperature. The vibrational temperature remains frozen downstream of the nozzle, except in the boundary layer.

Landry, John Gary↗

Computational Study of Flow Establishment in Hypersonic Pulse Facilities

This paper presents a study of the temporal evolution of the combustion flowfield established by the interaction of ram-accelerator-type projectiles with an explosive gas mixture accelerated to hypersonic speeds in an expansion tube. The Navier-Stokes equations for a chemically reacting gas are solved in a fully coupled manner using an implicit, time accurate algorithm. The solution procedure is based on a spatially second order, total variation diminishing (TVD) scheme and a temporally second order, variable-step, backward differentiation formula method. The hydrogen-oxygen chemistry is modeled with a 9-species, 19-step mechanism. The accuracy of the solution method is first demonstrated by several benchmark calculations. Numerical simulations of expansion tube flowfields are then presented for two different configurations. In particular, the development of the shock-induced combustion process is followed. In one case, designed to ensure ignition only in the boundary layer, the lateral extent of the combustion front during the initial transient phase was surprisingly large. The time histories of the calculated thrust and drag forces on the ram accelerator projectile are also presented.

Yungster, S.↗

Numerical Modeling of Pressurization of a Propellant Tank

An unsteady finite volume procedure has been developed to predict the history of pressure, temperature and mass flow rate of the pressurant and propellant during the expulsion of the propellant from a tank. The time dependent mass, momentum and energy conservation equations are solved at the ullage space. The model accounts for the change in the ullage volume due to expulsion of the propellant. It also accounts for the heat transfer from the tank wall and propellant to the ullage gas. The procedure was incorporated in the Generalized Fluid System Simulation Program (GFSSP). The results of several test cases were then compared with a published correlation of pressurant requirements for a given displacement of propellant. The agreement between the predictions and the correlation was found to be satisfactory.

Majumdar, Alok↗

Numerical Modeling of Pressurization of a Propellant Tank

An unsteady finite volume procedure has been developed to predict the history o pressure, temperature and mass flow rate of the pressurant and propellant during the expulsion of the propellant from a tan. The time dependent mass, momentum and energy conservation equations are solved at the ullage space. The model accounts for the change in the ullage volume due to expulsion of the propellant. It also accounts for the heat transfer from the tank wall and propellant to the ullage gas. The procedure was incorporated in the Generalized Fluid System Simulation Program (GFSSP). The results of several test cases were then compared with a published correlation of pressurant requirements for a given displacement of propellant. The agreement between the predictions and the correlation was found to be satisfactory.

Majumdar, Alok↗

Fe(+) Emission from Eta Carinae B and D Blobs: Theory and HST Observations

We present comparisons between numerical simulations of Fe II emission spectra and HST high-resolution and signal-to-noise spectra of the Eta Carinae Weigelt B & D blobs. The spectra show a wealth of FeII and [FeII] lines. We investigate the dependence of the spectrum on electron density and on pumping by blackbody radiation. Using the photoionization code Cloudy we were able successfully reproduce the fluxes of few hundred FeII & [Fe II] lines in the range from 4000A to 6000A. Fe II emission analysis in Eta Car ejecta is very important because it provides information about excitation conditions, the electron temperature and density solved in a selfconsistent way. Moreover, for the first time we were able to use Fe II emission analysis itself to put limits on CNO abundances in emitting regions.

Verner, E.↗

Robust Airfoil Optimization in High Resolution Design Space

The robust airfoil shape optimization is a direct method for drag reduction over a given range of operating conditions and has three advantages: (1) it prevents severe degradation in the off-design performance by using a smart descent direction in each optimization iteration, (2) it uses a large number of B-spline control points as design variables yet the resulting airfoil shape is fairly smooth, and (3) it allows the user to make a trade-off between the level of optimization and the amount of computing time consumed. The robust optimization method is demonstrated by solving a lift-constrained drag minimization problem for a two-dimensional airfoil in viscous flow with a large number of geometric design variables. Our experience with robust optimization indicates that our strategy produces reasonable airfoil shapes that are similar to the original airfoils, but these new shapes provide drag reduction over the specified range of Mach numbers. We have tested this strategy on a number of advanced airfoil models produced by knowledgeable aerodynamic design team members and found that our strategy produces airfoils better or equal to any designs produced by traditional design methods.

Li, Wu↗

Airfoil Design and Optimization by the One-Shot Method

An efficient numerical approach for the design of optimal aerodynamic shapes is presented in this paper. The objective of any optimization problem is to find the optimum of a cost function subject to a certain state equation (governing equation of the flow field) and certain side constraints. As in classical optimal control methods, the present approach introduces a costate variable (Lagrange multiplier) to evaluate the gradient of the cost function. High efficiency in reaching the optimum solution is achieved by using a multigrid technique and updating the shape in a hierarchical manner such that smooth (low-frequency) changes are done separately from high-frequency changes. Thus, the design variables are changed on a grid where their changes produce nonsmooth (high-frequency) perturbations that can be damped efficiently by the multigrid. The cost of solving the optimization problem is approximately two to three times the cost of the equivalent analysis problem.

Kuruvila, G.↗

An Implicit Upwind Algorithm for Computing Turbulent Flows on Unstructured Grids

An implicit, Navier-Stokes solution algorithm is presented for the computation of turbulent flow on unstructured grids. The inviscid fluxes are computed using an upwind algorithm and the solution is advanced in time using a backward-Euler time-stepping scheme. At each time step, the linear system of equations is approximately solved with a point-implicit relaxation scheme. This methodology provides a viable and robust algorithm for computing turbulent flows on unstructured meshes. Results are shown for subsonic flow over a NACA 0012 airfoil and for transonic flow over a RAE 2822 airfoil exhibiting a strong upper-surface shock. In addition, results are shown for 3 element and 4 element airfoil configurations. For the calculations, two one equation turbulence models are utilized. For the NACA 0012 airfoil, a pressure distribution and force data are compared with other computational results as well as with experiment. Comparisons of computed pressure distributions and velocity profiles with experimental data are shown for the RAE airfoil and for the 3 element configuration. For the 4 element case, comparisons of surface pressure distributions with experiment are made. In general, the agreement between the computations and the experiment is good.

Anerson, W. Kyle↗

Non-coincident Inter-instrument Comparisons of Ozone Measurements Using Quasi-conservative Coordinates

Ozone measurements from ozonesondes, AROTAL, DIAL, and POAM III instruments during the SOLVE-2/VINTERSOL period are composited in a time-varying, flow-following quasi-conservative (PV-6) coordinate space; the resulting composites from each instrument are mapped onto the other instruments locations and times. The mapped data are then used to intercompare data from the different instruments. Overall, the four data sets are found to be in good agreement. AROTAL shows somewhat lower values below 16 km, and DIAL has a positive bias at the upper limits of its altitude range. These intercomparisons are consistent with those obtained from more conventional near-coincident profiles, where available. Although the PV-theta mapping technique entails larger uncertainties of individual profile differences compared to direct near-coincident comparisons, the ability to include much larger numbers of comparisons can make this technique advantageous.

Lait, L. R.↗

Frequency Domain Calculations Of Acoustic Propagation

Two complex geometry problems are solved using the linearized Euler equations. The impedance mismatch method1 is used to impose the solid surfaces without the need to use a body-fitted grid. The problem is solved in the frequency domain to avoid long run times. Although the harmonic assumption eliminates all time dependence, a pseudo-time term is added to allow conventional iterative methods to be employed. A Jameson type, Runge-Kutta scheme is used to advance the solution in pseudo time. The spatial operator is based on a seven-point, sixth-order finite difference. Constant coefficient, sixth-derivative artificial dissipation is used throughout the domain. A buffer zone technique employing a complex frequency to damp all waves near the boundaries is used to minimize reflections. The results show that the method is capable of capturing the salient features of the scattering, but an excessive number of grid points are required to resolve the phenomena in the vicinity of the solid bodies because the wavelength of the acoustics is relatively short compared with the size of the bodies. Smoothly transitioning into the immersed boundary condition alleviates the difficulties, but a fine mesh is still required.

Lockard, David P.↗

NASA's Spaceflight Visual Impairment and Intracranial Hypertension Research Plan: An accelerated Research Collaboration

To date six long duration astronauts have experienced in flight visual changes and post flight signs of optic disc edema, globe flattening, choroidal folds, hyperoptic shifts and or raised intracranial pressure. In some cases the changes were transient while in others they are persistent with varying degrees of visual impairment. Given that all astronauts exposed to microgravity experience a cephalad fluid shift, and that both symptomatic and asymptomatic patients have exhibited optic nerve sheath edema on MRI, there is a high probability that all astronauts develop in-flight idiopathic intracranial hypertension to some degree. Those who are susceptible, have an increased likelihood of developing treatment resistant papilledema resulting in visual impairment and possible long-term vision loss. Such an acquired disability would have a profound mission impact and would be detrimental to the long term health of the astronaut. The visual impairment and increased intracranial pressure phenomenon appears to have multiple contributing factors. Consequently, the working "physiological fault bush" with elevated intracranial pressure at its center, is divided into ocular effects, and CNS and other effects. Some of these variables have been documented and or measured through operational data gathering, while others are unknown, undocumented and or hypothetical. Both the complexity of the problem and the urgency to find a solution require that a unique, non-traditional research model be employed such as the Accelerated Research Collaboration(TM) (ARC) model that has been pioneered by the Myelin Repair Foundation. In the ARC model a single entity facilitates and manages all aspects of the basic, translational, and clinical research, providing expert oversight for both scientific and managerial efforts. The result is a comprehensive research plan executed by a multidisciplinary team and the elimination of stove-piped research. The ARC model emphasizes efficient and effective communication between management and investigators; and real-time sharing of scientific discoveries in an effort to solve complex problems.

Otto, Christian↗

Conjugate Compressible Fluid Flow and Heat Transfer in Ducts

A computational approach to modeling transient, compressible fluid flow with heat transfer in long, narrow ducts is presented. The primary application of the model is for analyzing fluid flow and heat transfer in solid propellant rocket motor nozzle joints during motor start-up, but the approach is relevant to a wide range of analyses involving rapid pressurization and filling of ducts. Fluid flow is modeled through solution of the spatially one-dimensional, transient Euler equations. Source terms are included in the governing equations to account for the effects of wall friction and heat transfer. The equation solver is fully-implicit, thus providing greater flexibility than an explicit solver. This approach allows for resolution of pressure wave effects on the flow as well as for fast calculation of the steady-state solution when a quasi-steady approach is sufficient. Solution of the one-dimensional Euler equations with source terms significantly reduces computational run times compared to general purpose computational fluid dynamics packages solving the Navier-Stokes equations with resolved boundary layers. In addition, conjugate heat transfer is more readily implemented using the approach described in this paper than with most general purpose computational fluid dynamics packages. The compressible flow code has been integrated with a transient heat transfer solver to analyze heat transfer between the fluid and surrounding structure. Conjugate fluid flow and heat transfer solutions are presented. The author is unaware of any previous work available in the open literature which uses the same approach described in this paper.

Cross, M. F.↗

Using Markov Models of Fault Growth Physics and Environmental Stresses to Optimize Control Actions

A generalized Markov chain representation of fault dynamics is presented for the case that available modeling of fault growth physics and future environmental stresses can be represented by two independent stochastic process models. A contrived but representatively challenging example will be presented and analyzed, in which uncertainty in the modeling of fault growth physics is represented by a uniformly distributed dice throwing process, and a discrete random walk is used to represent uncertain modeling of future exogenous loading demands to be placed on the system. A finite horizon dynamic programming algorithm is used to solve for an optimal control policy over a finite time window for the case that stochastic models representing physics of failure and future environmental stresses are known, and the states of both stochastic processes are observable by implemented control routines. The fundamental limitations of optimization performed in the presence of uncertain modeling information are examined by comparing the outcomes obtained from simulations of an optimizing control policy with the outcomes that would be achievable if all modeling uncertainties were removed from the system.

Bole, Brian↗

A Summary of the Space-Time Conservation Element and Solution Element (CESE) Method

The space-time Conservation Element and Solution Element (CESE) method for solving conservation laws is examined for its development motivation and design requirements. The characteristics of the resulting scheme are discussed. The discretization of the Euler equations is presented to show readers how to construct a scheme based on the CESE method. The differences and similarities between the CESE method and other traditional methods are discussed. The strengths and weaknesses of the method are also addressed.

Fluid Dynamics↗

Unexpected Nonlinear Effects in Superconducting Transition-Edge Sensors

When a normal metal transitions into the superconducting state the DC resistance drops from a finite value to zero over some finite transition width in temperature, current, and magnetic field. Superconducting transition-edge sensors (TESs) operate within this transition region and uses resistive changes to measure deposited thermal energy. This resistive transition is not perfectly smooth and a wide range of TES designs and materials show sub-structure in the resistive transition (as seen in smooth nonmonotonic behavior, jump discontinuities, and hysteresis in the devices current-voltage relation and derivatives of the resistance with respect to temperature, bias current, and magnetic field). TES technology has advanced to the point where for many applications this structure is the limiting factor in performance and optimization consists of finding operating points away from these structures. For example, operating at or near this structure can lead to nonlinearity in the detectors response and gain scale, limit the spectral range of the detector by limiting the usable resistive range, and degrade energy resolution. The origin of much of this substructure is unknown. This presentation investigates a number of possible sources in turn. First we model the TES as a superconducting weak-link and solve for the characteristic differential equations current and voltage time dependence. We find:(1) measured DC biased current-voltage relationship is the time-average of a much higher frequency limit cycle solution.(2) We calculate the fundamental frequency and estimate the power radiated from the TES treating the bias leads as an antennae.(3) The solution for a set of circuit parameters becomes multivalued leading to current transitions between levels.(4)The circuit parameters can change the measure resistance and mask the true critical current. As a consequence the TES resistance surface is not just a function of temperature, current, and magnetic field but is also a function of the circuit elements (such as shunt resistor, SQUID inductance, and capacitor values). In other words, same device measured in different electrical circuits will have a different resistive surface in temperature, current, and magnetic field. Next we consider that at the transition temperature of a superconductor both the magnetic penetration depth and coherence length are divergent. As a consequence these important characteristic length scales are changing with operating point. We present measurements on devices showing commensurate behavior between these characteristic lengths and the length scale of added normal metal structures. Reordering of proximity vortices leads to discontinuities and irreversibility of the current-voltage curves. Last we consider a weak-link TES including both thermal activated resistance effects and the effect of the magnetic penetration depth being a function of temperature and magnetic field. We derive its impact on the resistive transition surface and the important device parameters a and b.

Superconducting Detectors↗

Rapid Assessment of a Lunar Surface Laboratory Module

Johnson Space Center is piloting an innovation team known as the Forge to perform rapid turnaround, concurrent engineering studies. The Forge is intended to be able to pull together a team of strategically selected domain experts to innovatively solve a specified problem within a short period of time. The Forge’s inaugural project being the development of a Lunar Surface Laboratory Module concept. The lab module is intended as an option for the Artemis program’s Lunar Architecture Team. The lab module may allow for increased science productivity and return on investment than a surface architecture with only the habitat. This paper will discuss the configuration of the Lunar Laboratory module, including the pressure vessel, its placement relative to the Surface Habitat and Pressurized Rover, and internal architecture. In order to determine the suitability of such a design, a human factors evaluation is conducted. Historically, a human factors evaluation of a habitable volume could require a month or longer, including time needed to prepare a test plan, conduct the evaluation, compile the data, and produce a report. The entire Forge effort will only last for five weeks, forcing the creation of a rapid assessment process. This rapid assessment includes three evaluation studies, all conducted and processed in the space of a single day, to determine the acceptability of the resulting design. The first evaluation is a card sort evaluation of the laboratory instruments. A brief discussion of the card sort evaluation is provided. This card sort is focused on the relative positioning of the laboratory instruments, determining whether a given instrument should be adjacent to, or separated from, another. It will use human factors personnel and both life and physical scientists as subject matter experts. The second evaluation is a science productivity tabletop evaluation. This evaluation focuses on the selection of instruments and the number of instruments, attempting to assess the amount of work the crew will be able to use the facility to do. It uses life and physical science personnel as subject matter experts. The third evaluation is a CAD or Virtual Reality human factors walkthrough evaluation. This evaluation focuses on the layout of the lab module and its connectivity to the basecamp’s pressurized elements. If possible, astronauts will be used as subject matter experts. Each of these evaluations will focus primarily on objective data with limited subjective comments, reducing the time required for processing responses. The results of each evaluation and their implications will be discussed. Based on the evaluation data, an overall acceptability of the Lunar Surface Laboratory Module produced by the Forge study can be determined. Lessons learned and recommendations for future use of the human factors rapid assessment process in early design trades will also be provided.

Lunar Laboratory↗

A Rapid Assessment of a Lunar Surface Laboratory Module

Johnson Space Center has piloted an innovation team known as the Forge to perform rapid turnaround, concurrent engineering studies. The Forge is intended to be able to pull together a team of strategically selected domain experts to innovatively solve a specified problem within a short period of time. The Forge’s inaugural project was the development of a Lunar Surface Laboratory Module concept. The lab module is intended as an option for the NASA Lunar Architecture Team and may allow for increased science productivity and return on investment than with only the Surface Habitat. The Forge team developed two concepts, a stationary laboratory and a mobile laboratory. Human factors evaluations have been conducted to determine the acceptability of the resulting design concepts. The human factors evaluations conducted include a card sort evaluation of the laboratories by science discipline and the System Usability Scale (SUS). The card sort focused subject matter expert expectations for the relative positioning of the laboratory functional areas, determining whether a given lab area should be adjacent to, or separated from, another. The SUS focused on layout efficiency. Both used human factors personnel and both life and physical scientists as subject matter experts. Each of these evaluations will focus primarily on objective data with limited subjective comments, reducing the time required for processing responses. The results of the card sort and SUS evaluations and their implications will be discussed. Based on the evaluation data, an overall acceptability of the Lunar Surface Laboratory Module produced by the Forge study can be determined. Lessons learned and recommendations for future use of human factors rapid assessment processes in early design trades will also be provided.

Forge↗