Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “JuMP model”

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 37 records · Page 2

Integration of the Total Lightning Jump Algorithm into Current Operational Warning Environment Conceptual Models

The presence and rates of total lightning are both correlated to and physically dependent upon storm updraft strength, mixed phase precipitation volume and the size of the charging zone. The updraft modulates the ingredients necessary for electrification within a thunderstorm, while the updraft also plays a critical role in the development of severe and hazardous weather. Therefore utilizing this relationship, the monitoring of lightning rates and jumps provides an additional piece of information on the evolution of a thunderstorm, more often than not, at higher temporal resolution than current operational radar systems. This correlation is the basis for the total lightning jump algorithm that has been developed in recent years. In order to become a viable option for operational forecasters to incorporate into their severe storm monitoring process, the total lightning jump must be placed into the framework of several severe storm conceptual models (e.g., radar evolution, storm morphology) which forecasters have built through training and experience. Thus, one of the goals of this study is to examine and relate the lightning jump concept to often used radar parameters (e.g., dBZ vertical structure, VIL, MESH, MESO/shear) in the warning environment. Tying lightning trends and lightning jump occurrences to these radar based parameters will provide forecasters with an additional tool that they can use to build an accurate realtime depiction as to what is going on in a given environment. Furthermore, relating the lightning jump concept to these parameters could also increase confidence in a warning decision they have already made, help tip the scales on whether or not to warn on a given storm, or to draw the forecaster s attention to a particular storm that is rapidly developing. Furthermore the lightning information will add vital storm scale information in regions that are not well covered by radar, or when radar failures occur. The physical basis for the lightning jump algorithm in relation to severe storm dynamics and microphysics is a key component that must be further explored. Many radar studies have examined flash rates and their relation to updraft strength, updraft volume, precipitation -sized ice mass, etc.; however, very few have related the concept of the lightning jump and manifestation of severe weather to storm dynamics and microphysics using multi -Doppler and polarimetric radar techniques. Therefore, the second half of this study will combine the lightning jump algorithm and these radar techniques in order to place the lightning jump concept into a physical and dynamical framework. This analysis includes examining such parameters as mixed phase precipitation volume, charging zone, updraft strength and updraft volume. Such a study should provide increased understanding of and confidence in the strengths and limitations of the lightning jump algorithm in the storm warning process.

Schultz, Christopher J.↗

Pressure distributions and shock shapes for a bent-nose biconic at incidence

Initial results of a NASA study of the lift-drag characteristics of a 12.84/7 deg biconic model intended for airbraking during atmospheric entry of probes to Mars, Venus, Saturn, and Titan are reported. Pressure distributions and shock shapes were measured in the Langley 20 in. Mach 6 tunnel with the spherically blunted bent-nose model set at angles from 0-25 deg. Pressure distributions and shock shapes where computed using the STEIN flowfield code, which features a MacCormack scheme to integrate the three-dimensional Euler equations, the Rankine-Hugoniot jump conditions to model shock waves as discontinuities, and requires a supersonic condition at every step. A comparison was made between measured and predicted values. The leeward shock angle was found to be predictable to within 3% for all angles of attack, while parabolized Navier-Stokes equations are regarded as offering more accurate results than the STEIN code for surface pressure distributions.

Miller, C. G., III↗

Application of unsteady aerodynamic methods for transonic aeroelastic analysis

Aerodynamic methods for aeroelastic analysis are applied to various flow problems. These methods include those that solve the three dimensional transonic small disturbance (TSD) potential equation, the two dimensional (2-D) full potential (FP) equation, and the 2-D thin layer Navier-Stokes equations. Flutter analysis performed using TSD aerodynamics show that such methods can be used to analyze some aeroelastic phenomena. For thicker bodies and larger amplitude motions, a nonisentropic FP method is presented. The unsteady FP equation is modified to model the entropy jumps across shock waves. The conservative form of the modified equation is solved in generalized coordinates using an implicit, approximate factorization method. Pressures calculated on the NLR 7301 and NACA 64A010A airfoils using the nonisentropic FP method are presented. It is shown that modeling shock generated entropy extends the range of validity of the FP method. A Navier-Stokes code is correlated with pressures measured on a supercritical airfoil at transonic speeds. When corrections are made for wind tunnel wall effects, the calculations correlate well with the measured data.

Whitlow, Woodrow, Jr.↗

From S-N to the Paris Law with a New Mixed-Mode Cohesive Fatigue Model

The relationship between fatigue life and fatigue crack propagation rate is explored with a new cohesive damage model. The parameters of the model are obtained from idealizations of S-N diagrams used in engineering design. The model is based on the hypothesis that both stable tearing damage and damage due to cyclic loading are representations of a density of microcracks and, therefore, a single damage variable can describe the state of damage. This assumption implies that the quasi-static cohesive law that describes tearing is also the envelope of the fatigue damage. Fatigue damage within the cohesive envelope is assumed to accumulate at a rate that depends on the displacement jumps. The fatigue model was implemented as a UMAT subroutine for Abaqus cohesive elements by adding fatigue damage accumulation to a cohesive model based on the Turon mixed-mode cohesive laws. The analyses were conducted using a simplified cyclic loading procedure in which the maximum applied load is kept constant and the computational expense of cycling the load is avoided. The predicted propagation rates in mode I and mixed mode were compared to experimental results for IM7/8552 graphite/epoxy tape. Several aspects of the results were investigated, including the effect of R-curves, the stress ratio R, and the difference between displacement control and force control.

Davila, Carlos G.↗

A Cohesive Fatigue Model Based on the S-N Diagram

The relationship between fatigue life and fatigue crack propagation rate is explored with a new cohesive damage model. The parameters of the model are obtained from idealizations of S-N diagrams used in engineering design. The model is based on the hypothesis that both stable tearing damage and damage due to cyclic loading are representations of a density of microcracks and, therefore, a single damage variable can describe the state of damage. This assumption implies that the quasi-static cohesive law that describes tearing is also the envelope of the fatigue damage. Fatigue damage within the cohesive envelope is assumed to accumulate at a rate that depends on the displacement jumps. The fatigue model was implemented as a UMAT subroutine for Abaqus cohesive elements by adding fatigue damage accumulation to a cohesive model based on the Turon quasi-static cohesive law. The analyses were conducted using a simplified cyclic loading procedure in which the maximum applied load is kept constant and the computational expense of cycling the load is avoided. The predicted propagation rates for a double cantilever beam (DCB) specimen were compared to experimental results for IM7/8552 graphite/epoxy tape.

Dávila, Carlos G.↗

Application of a nonisentropic full potential method to AGARD standard airfoils

An entropy-correction method for the unsteady full potential equation is presented. The unsteady potential equation is modified to model the entropy jumps across shock waves. The conservative form of the modified equation is solved in generalized coordinates using an implicit, approximate factorization method. A flux-biasing differencing method, which generates the proper amounts of artificial viscosity in supersonic regions, is used to discretize the flow equations in space. Calculated results are presented for the NLR 7301, NACA 0012, and NACA 64A010A airfoils. Comparisons of the present method and solutions of the Euler equations are presented for the NLR 7301 airfoil, and comparisons of the present method and experimental data are presented for all three airfoils. The comparisons show that the present method more accurately models solutions of the Euler equations and experiment than does the isentropic potential formulation. In addition, it is shown that modeling shock-generated entropy extends the range of validity of the full potential method.

Whitlow, Woodrow, Jr.↗

Application of a nonisentropic full potential method to AGARD standard airfoils

An entropy-correction method for the unsteady full potential equation is presented. The unsteady potential equation is modified to model the entropy jumps across shock waves. The conservative form of the modified equation is solved in generalized coordinates using an implicit, approximate factorization method. A flux-biasing differencing method, which generates the proper amounts of artificial viscosity in supersonic regions, is used to discretize the flow equations in space. Calculated results are presented for the NLR 7301, NACA 0012, and NACA 64A010A airfoils. Comparisons of the present method and solutions of the Euler equations are presented for the NLR 7301 airfoil, and comparisons of the present method and experimental data are presented for all three airfoils. The comparisons show that the present method more accurately models solutions of the Euler equations and experiment than does the isentropic potential formulation. In addition, it is shown that modeling shock-generated entropy extends the range of validity of the full potential method.

Whitlow, Woodrow, Jr.↗

Solar wind stream interfaces

Results are presented for a superposed epoch analysis of discontinuous solar wind interfaces. The average time-space profiles of stream interfaces are discussed with reference to fluid properties (flow speed, pressure ridge, density, electron and proton temperatures) and kinetic properties (electron core and halo, flow speed fluctuations, electron heat flux, alpha particles). Other aspects of stream interfaces are described, such as the persistence of individual interfaces, shock associations, the sector boundaries of the interplanetary magnetic field, and sudden impulses in the geomagnetic field. Interface position is considered in terms of the observed temperature jump. A conceptual model of high-speed stream evolution is proposed.

Gosling, J. T.↗

Post-glacial relaxation of a viscously stratified compressible mantle

The postglacial relaxation phenomenon of a viscously compressible nonself-gravitating spherical-shell model is investigated. Analytical solutions to an exponentially depth-dependent viscosity are developed with various types of compressible density models such as those with an exponential dependence of the radius and an algebraic-root dependence of the radius. The solutions of a viscously compressible multilayered model with constant thermodynamic properties and viscosity are developed by the propagator matrix method. The results show that inferences of deep mantle viscosity from postglacial rebound would be hampered by mantle compressibility for long-wavelength harmonics because of the smaller excitation of compressible eigenfunctions in the lower mantle. The relaxation times and velocity fields are more sensitive to higher viscosity contrast for the exponentially varying viscosity than for models with discrete jumps in the viscosity structure.

Wu, J.↗

Stochastic Stability of Sampled Data Systems with a Jump Linear Controller

In this paper an equivalence between the stochastic stability of a sampled-data system and its associated discrete-time representation is established. The sampled-data system consists of a deterministic, linear, time-invariant, continuous-time plant and a stochastic, linear, time-invariant, discrete-time, jump linear controller. The jump linear controller models computer systems and communication networks that are subject to stochastic upsets or disruptions. This sampled-data model has been used in the analysis and design of fault-tolerant systems and computer-control systems with random communication delays without taking into account the inter-sample response. This paper shows that the known equivalence between the stability of a deterministic sampled-data system and the associated discrete-time representation holds even in a stochastic framework.

Gonzalez, Oscar R.↗

3D Material Response of the MSL Heatshield Using NuSil-Coated PICA

The Mars Science Laboratory (MSL) Entry, Descent and Landing Instrumentation (MEDLI) collected in-flight data largely used by the ablation community to verify and validate physics-based models for the response of the Phenolic Impregnated Carbon Ablator (PICA) material [1-4]. MEDLI data were recently used to guide the development of NASA’s high-fidelity material response models for PICA implemented in the Porous material Analysis Toolbox based on OpenFOAM (PATO) software [5-6]. A follow-up instrumentation suite, MEDLI2, was used for the Mars 2020 mission [7] after the large scientific impact of MEDLI. Recent analyses performed as part of MEDLI2 development drew the attention to significant effects of a protective coating to the aerothermal response of PICA. NuSil, a silicone-based overcoat sprayed onto the MSL heatshield, including the MEDLI plugs, as contamination control, is currently neglected in PICA ablation models. To mitigate the spread of phenolic dust from PICA, NuSil was applied to the entire MSL and Mars 2020 heatshields, including the MEDLI and MEDLI2 plugs. Ground testing of PICA-NuSil (PICA-N) models exhibited surface temperature jumps due to oxide scale formation and subsequent NuSil burn- off. It is therefore critical to include a model for the aerothermal response of the coating in ongoing code development and validation efforts. Figure 1 illustrates the PICA-N model implemented in PATO with a test case that shows a temperature jump and a change of recession rate at the wall after 20 seconds of simulation when the NuSil coating was fully removed.

Aerospace↗

A Computational Model for Rotor-Fuselage Interactional Aerodynamics

A novel unsteady rotor-fuselage interactional aerodynamics model has been developed. This model loosely couples a Generalized Dynamic Wake Theory (GDWT) to a thin-layer Navier-Stokes solution procedure. This coupling is achieved using an unsteady pressure jump boundary condition in the Navier-Stokes model. The new unsteady pressure jump boundary condition models each rotor blade as a moving pressure jump which travels around the rotor azimuth and is applied between two adjacent planes in a cylindrical, non-rotating grid. Comparisons are made between measured and predicted time-averaged and time-accurate rotor inflow ratios. Additional comparisons are made between measured and predicted unsteady surface pressures on the top centerline and sides of the fuselage.

D. Douglas Boyd, Jr.↗

PDF turbulence modeling and DNS

The problem of time discontinuity (or jump condition) in the coalescence/dispersion (C/D) mixing model is addressed in probability density function (pdf). A C/D mixing model continuous in time is introduced. With the continuous mixing model, the process of chemical reaction can be fully coupled with mixing. In the case of homogeneous turbulence decay, the new model predicts a pdf very close to a Gaussian distribution, with finite higher moments also close to that of a Gaussian distribution. Results from the continuous mixing model are compared with both experimental data and numerical results from conventional C/D models. The effect of Coriolis forces on compressible homogeneous turbulence is studied using direct numerical simulation (DNS). The numerical method used in this study is an eight order compact difference scheme. Contrary to the conclusions reached by previous DNS studies on incompressible isotropic turbulence, the present results show that the Coriolis force increases the dissipation rate of turbulent kinetic energy, and that anisotropy develops as the Coriolis force increases. The Taylor-Proudman theory does apply since the derivatives in the direction of the rotation axis vanishes rapidly. A closer analysis reveals that the dissipation rate of the incompressible component of the turbulent kinetic energy indeed decreases with a higher rotation rate, consistent with incompressible flow simulations (Bardina), while the dissipation rate of the compressible part increases; the net gain is positive. Inertial waves are observed in the simulation results.

Hsu, A. T.↗

Modeling and Simulation of Tank Pressure Control using Zero-Boil Off Active Thermal Control for LOXSAT Technology Demonstration Mission

To-date, research and modeling of cryogenic fluid management technologies (CFM) for spaceflight has been limited to ground tests, short-duration zero-g simulations (e.g. drop towers), and small-scale experiments on-orbit. There has not been a large-scale flight demonstration of a flight-like system. As future NASA missions to take humans further from Earth will require large, in-space cryogenic propulsion vehicles, it is imperative to begin collecting flight data for these systems to accurately model and design future vehicles. To meet this goal, NASA awarded tipping point technology demonstration awards to Eta Space, Lockheed Martin, Space Exploration Technologies (SpaceX), and United Launch Alliance (ULA) to demonstrate on-orbit storage and transfer of cryogenic propellant. For its award, Eta Space is developing LOXSAT-1. It is a small satellite that will be launched on a Rocket Lab Electron rocket. The spacecraft consists of a Rocket Lab Photon spacecraft bus with a primary payload of a spherical liquid oxygen (LOX) storage tank with thermodynamic control systems. The primary objective of the mission is to demonstrate zero-boil-off storage of liquid oxygen. To accomplish this objective, the payload is equipped with an active thermal control system fluid loop that consists of propellant management device (PMD), positive-displacement pump, heat exchanger connected to a cryocooler, and a spray bar mixing injector. When the fluid loop is operating, fluid is drawn from the tank by the PMD and pumped through the heat exchanger, lowering the fluid temperature below the fluid temperature in the tank. This subcooled liquid is then injected back into the tank through the spray bar. The subcooled injected liquid has two effects. If it is sprayed into the ullage space, the injected liquid will form into jets or droplets and exchange heat with the ullage gas. This will cool and condense the gas, reducing the pressure in the tank. Additionally, the liquid that is not sprayed through the ullage, as well as any remaining liquid spray from the ullage, will rejoin the liquid mass of the tank, lowering the bulk temperature of the liquid. These combines effects provide for zero-boil off pressure control by lowering the tank pressure and the liquid saturation pressure simultaneously, ensuring the liquid stays subcooled. To model these complex mechanics and predict the performance of the active thermal control system, NASA is providing Eta Space with 3 parallel models of the tank. The approach of providing 3 different models allows for cross-checking and comparisons between the three to better understand how different modeling assumptions and selection semi-empirical factors affects the modeling result. Additionally, developing 3 models provides three different schemes for numerical simulation, providing confidence that results depict real physical phenomenon and not numerical quirks of the program. Within the tank thermodynamics, there are two primary areas of heat transfer we concern ourselves with: the heat transfer between the ullage space and the droplet spray, and between the ullage space and bulk liquid. For the droplet heat transfer, there are multiple sets of assumptions that can be made and correlations that can be used. Currently, two working models - the TankSIM model and Easy5 model – provide for an overview of the different approaches available. The TankSIM model and Easy5 model use two different models for droplet heating and evaporation that illustrate how the models use different types of mechanisms to arrive at the same answer. For the TankSIM model, droplets are treated as spheres of constant radius. Heat is transferred from the ullage to the droplet and warms the droplet until it reaches saturation, then the droplet begins evaporating and reducing its radius and mass. To calculate the heat transfer coefficient between the droplet and ullage, the Ranz-Marshall correlation is used. To determine the number of droplets in the ullage, a resident mass approach is used. This approach averages the number of droplets such that residuals at the start-up and shut-off of the spray bar cancel out. This same approach is used in the Easy5 model. The GFSSP model implements a linked list to track individual droplet “nodes” within the model. For the Easy5 model, the droplet is assumed to have an interface at a temperature equal to the saturation temperature corresponding to the pressure of the gas phase. The heat transfer from the gas to the interface and the interface to the droplet bulk is then calculated, and the net mass transfer between the droplet and interface is determined by performing an energy balance across the interface. For the gas side of the interface, the Ranz-Marshall correlation is used. For the liquid side of the interface, a variety of correlations were tried, including Kronig and Brink (1950) and effective conductivity models. As a result of these assumptions, the Easy5 model currently predicts faster depressurization, as at saturated vapor conditions, the heat transfer coeffect on the liquid side for the Easy5 model is greater than the heat transfer coefficient predicted by Ranz-Marshall used in the TankSIM code. This greater heat flux translates into faster condensation of the saturated ullage gas. At the bulk liquid to ullage interface, the models are in much closer agreement. Both models model the ullage as a sphere centered within the bulk liquid in the tank, and both use the energy-jumping boundary condition to model heat and mass transfer across the interface. There are slight differences in how the interface temperature is calculated, however. The Easy5 model assumes the temperature of the interface is equal to the saturation temperature associated with the pressure of the gas phase. The TankSIM model calculates this temperature with Alabovskii’s equation, which provides a correction factor for interface temperatures. Analysis tasks are focused on determining rates of depressurization within the tank during active cooling operation. To maintain net positive suction head at the pump inlet, the tank pressure cannot fall faster than the saturation pressure associated with the temperature of the bulk liquid. Additionally, there is interest in analyzing the performance of the loop at different pump speeds and cryocooler input powers. Adjusting the flowrate affects both the performance of the heat exchanger between the cryocooler and pumped liquid, and the heat transfer between the droplet spray and the ullage. Ideally, a pump speed and cryocooler power can be selected that will allow the tank to operate in zero-boil-off mode with a very narrow range of storage pressure.

Cameron J. Hines↗

Modeling and Simulation of Tank Pressure Control using Zero-Boil Off Active Thermal Control for LOXSAT Technology Demonstration Mission

To-date, research and modeling of cryogenic fluid management technologies (CFM) for spaceflight has been limited to ground tests, short-duration zero-g simulations (e.g. drop towers), and small-scale experiments on-orbit. There has not been a large-scale flight demonstration of a flight-like system. As future NASA missions to take humans further from Earth will require large, in-space cryogenic propulsion vehicles, it is imperative to begin collecting flight data for these systems to accurately model and design future vehicles. To meet this goal, NASA awarded tipping point technology demonstration awards to Eta Space, Lockheed Martin, Space Exploration Technologies (SpaceX), and United Launch Alliance (ULA) to demonstrate on-orbit storage and transfer of cryogenic propellant. For its award, Eta Space is developing LOXSAT-1. It is a small satellite that will be launched on a Rocket Lab Electron rocket. The spacecraft consists of a Rocket Lab Photon spacecraft bus with a primary payload of a spherical liquid oxygen (LOX) storage tank with thermodynamic control systems. The primary objective of the mission is to demonstrate zero-boil-off storage of liquid oxygen. To accomplish this objective, the payload is equipped with an active thermal control system fluid loop that consists of propellant management device (PMD), positive-displacement pump, heat exchanger connected to a cryocooler, and a spray bar mixing injector. When the fluid loop is operating, fluid is drawn from the tank by the PMD and pumped through the heat exchanger, lowering the fluid temperature below the fluid temperature in the tank. This subcooled liquid is then injected back into the tank through the spray bar. The subcooled injected liquid has two effects. If it is sprayed into the ullage space, the injected liquid will form into jets or droplets and exchange heat with the ullage gas. This will cool and condense the gas, reducing the pressure in the tank. Additionally, the liquid that is not sprayed through the ullage, as well as any remaining liquid spray from the ullage, will rejoin the liquid mass of the tank, lowering the bulk temperature of the liquid. These combines effects provide for zero-boil off pressure control by lowering the tank pressure and the liquid saturation pressure simultaneously, ensuring the liquid stays subcooled. To model these complex mechanics and predict the performance of the active thermal control system, NASA is providing Eta Space with 3 parallel models of the tank. The approach of providing 3 different models allows for cross-checking and comparisons between the three to better understand how different modeling assumptions and selection semi-empirical factors affects the modeling result. Additionally, developing 3 models provides three different schemes for numerical simulation, providing confidence that results depict real physical phenomenon and not numerical quirks of the program. Within the tank thermodynamics, there are two primary areas of heat transfer we concern ourselves with: the heat transfer between the ullage space and the droplet spray, and between the ullage space and bulk liquid. For the droplet heat transfer, there are multiple sets of assumptions that can be made and correlations that can be used. Currently, two working models - the TankSIM model and Easy5 model – provide for an overview of the different approaches available. The TankSIM model and Easy5 model use two different models for droplet heating and evaporation that illustrate how the models use different types of mechanisms to arrive at the same answer. For the TankSIM model, droplets are treated as spheres of constant radius. Heat is transferred from the ullage to the droplet and warms the droplet until it reaches saturation, then the droplet begins evaporating and reducing its radius and mass. To calculate the heat transfer coefficient between the droplet and ullage, the Ranz-Marshall correlation is used. To determine the number of droplets in the ullage, a resident mass approach is used. This approach averages the number of droplets such that residuals at the start-up and shut-off of the spray bar cancel out. This same approach is used in the Easy5 model. The GFSSP model implements a linked list to track individual droplet “nodes” within the model. For the Easy5 model, the droplet is assumed to have an interface at a temperature equal to the saturation temperature corresponding to the pressure of the gas phase. The heat transfer from the gas to the interface and the interface to the droplet bulk is then calculated, and the net mass transfer between the droplet and interface is determined by performing an energy balance across the interface. For the gas side of the interface, the Ranz-Marshall correlation is used. For the liquid side of the interface, a variety of correlations were tried, including Kronig and Brink (1950) and effective conductivity models. As a result of these assumptions, the Easy5 model currently predicts faster depressurization, as at saturated vapor conditions, the heat transfer coeffect on the liquid side for the Easy5 model is greater than the heat transfer coefficient predicted by Ranz-Marshall used in the TankSIM code. This greater heat flux translates into faster condensation of the saturated ullage gas. At the bulk liquid to ullage interface, the models are in much closer agreement. Both models model the ullage as a sphere centered within the bulk liquid in the tank, and both use the energy-jumping boundary condition to model heat and mass transfer across the interface. There are slight differences in how the interface temperature is calculated, however. The Easy5 model assumes the temperature of the interface is equal to the saturation temperature associated with the pressure of the gas phase. The TankSIM model calculates this temperature with Alabovskii’s equation, which provides a correction factor for interface temperatures. Analysis tasks are focused on determining rates of depressurization within the tank during active cooling operation. To maintain net positive suction head at the pump inlet, the tank pressure cannot fall faster than the saturation pressure associated with the temperature of the bulk liquid. Additionally, there is interest in analyzing the performance of the loop at different pump speeds and cryocooler input powers. Adjusting the flowrate affects both the performance of the heat exchanger between the cryocooler and pumped liquid, and the heat transfer between the droplet spray and the ullage. Ideally, a pump speed and cryocooler power can be selected that will allow the tank to operate in zero-boil-off mode with a very narrow range of storage pressure.

zero boil-off↗