Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “convergence rate”

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,171 records · Page 65

California-Kepler Survey. IX. Revisiting the Minimum-mass Extrasolar Nebula with Precise Stellar Parameters

We investigate a possible correlation between the solid surface density Σ of the minimum-mass extrasolar nebula (MMEN) and the host star mass M {sub ⋆} and metallicity [Fe/H]. Leveraging on the precise host star properties from the California-Kepler Survey (CKS), we found that Σ=50{sub −20}{sup +33} g cm{sup −2} (a/1 au){sup −1.75±0.07} (M {sub ⋆}/M {sub ⊙}){sup 1.04±0.22} 10{sup 0.22±0.05[Fe/H]} for Kepler-like systems (1–4R {sub ⊕}; a < 1 au). The strong M {sub ⋆} dependence is reminiscent of previous dust continuum results that the solid disk mass scales with M {sub ⋆}. The weaker [Fe/H] dependence shows that sub-Neptune planets, unlike giant planets, form readily in lower metallicity environment. The innermost region (a < 0.1 au) of an MMEN maintains a smooth profile despite a steep decline of planet occurrence rate: a result that favors the truncation of disks by corotating magnetospheres with a range of rotation periods, rather than the sublimation of dust. The Σ of Kepler multitransiting systems shows a much stronger correlation with M {sub ⋆} and [Fe/H] than singles. This suggests that the dynamically hot evolution that produced single systems also partially removed the memory of formation in disks. Radial-velocity planets yielded a MMEN very similar to CKS planets; transit-timing-variation planets’ postulated convergent migration history is supported by their poorly constrained MMEN. We found that lower mass stars have a higher efficiency of forming/retaining planets: for Sun-like stars, about 20% of the solid mass within ∼1 au are converted/preserved as sub-Neptunes, compared to 70% for late-K to early-M stars. This may be due to the lower binary fraction, lower giant-planet occurrence, or the longer disk lifetime of lower mass stars.

79 ASTRONOMY AND ASTROPHYSICS↗

IGA-MPM: The Isogeometric Material Point Method

In this work, we propose the use of Isogeometric Analysis (IGA) within the context of the Material Point Method (MPM), and refer to the approach as IGA-MPM. We use the idea of IGA, and its instantiation based on Non-Uniform Rational B-Splines (NURBS), to build higher-order accurate and smooth approximation for MPM. Higher-order smoothness yields a continuous representation of the strain rate, and, as a result, prevents jumps in the stress and other history variables as the material points cross the element boundaries. Furthermore, NURBS can exactly represent all conic sections and the corresponding symmetries in the solution, which may be important in some applications. Several numerical examples of increasing complexity are presented, and show the ability of IGA-MPM to eliminate the well known cell-crossing instability of the conventional MPM. In addition, the examples presented demonstrate improved accuracy, convergence, and symmetry preservation of IGA-MPM compared to the conventional MPM, both for rectilinear and curved geometries.

42 ENGINEERING↗

Catalytic Resonance Theory: Forecasting the Flow of Programmable Catalytic Loops

Chemical transformations on catalyst surfaces occur through series and parallel reaction pathways. These complex networks and their behavior can be most simply evaluated through a three-species surface reaction loop (A* to B* to C* to A*) that is internal to the overall chemical reaction. Application of an oscillating dynamic catalyst to this reactive loop has been shown to exhibit one of three types of behavior: (1) a positive net flux of molecules about the loop in the clockwise direction, (2) a negative net flux of molecules about the loop in the counterclockwise direction, or (3) negligible flux of molecules about the loop at the limit cycle of reaction. Three-species surface loops were simulated with microkinetic modeling to assess the reaction loop behavior resulting from a catalytic surface oscillating between two or more catalyst surface energy states. Selected input parameters for the simulations spanned an 11-dimensional parameter space using 127 688 different parameter combinations. Their converged limit cycle solutions were analyzed for their loop turnover frequencies, the majority of which were found to be approximately zero. Classification and regression machine learning models were trained to predict the sign and magnitude of the loop turnover frequency and successfully performed above accessible baselines. Notably, the classification models exhibited a baseline weighted F1 score of 0.49, whereas trained models achieved weighted F1 scores of 0.94 and 0.96 when trained on the parameters used to define the simulations and derived rate constants, respectively. The trained models successfully predicted catalytic loop behavior, and interpretation of these models revealed all input parameters to be important for the prediction and performance of each model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Dynamics and energetics of the South Pacific Convergence Zone during FGGE SOP-1

The major research objectives are to diagnose the physical processes responsible for the maintenance of the South Pacific Convergence Zone (SPCZ) and to examine the role of the SPCZ in the large-scale circulation patterns of the Southern Hemisphere. To accomplish these objectives researchers used several data sources which include: a modified set of Level III-b upper air analyses, originally produced by ECMWF (Vincent, 1982); subjectively analyzed surface analyses for the South Pacific based on island station reports (Vincent, 1985); outgoing longwave radiation values supplied to us by NOAA/NESDIS; and equivalent black body temperatures and precipitation rates derived by Robertson. In the past year researchers found that wave number four plays an inportant role in the Southern Hemisphere tropics during the 15-day period when the sPCZ was a dominant feature, particularly with regard to the baroclinic conversion of potential to kinetic energy (Huang and Vincent, 1985). The convectively-active SPCZ area was found to make a significant contribution to this conversion process; thus, it appears that baroclinic effects and latent heating are important in maintaining the SPCZ. Recently efforts concentrated on two research tasks, an examination of cyclone activity within the SPCZ (Kann, 1985; Vincent, 1985; Vincent and Kann, 1985) and a study of the heat and moisture budgets in the South Pacific (Miller, et al., 1985). It was found that cyclonic disturbances occurred with regularity in the Zone from 10 to 17 January.

Vincent, D. G.↗

New higher-order Godunov code for modelling performance of two-stage light gas guns

A new quasi-one-dimensional Godunov code for modeling two-stage light gas guns is described. The code is third-order accurate in space and second-order accurate in time. A very accurate Riemann solver is used. Friction and heat transfer to the tube wall for gases and dense media are modeled and a simple nonequilibrium turbulence model is used for gas flows. The code also models gunpowder burn in the first-stage breech. Realistic equations of state (EOS) are used for all media. The code was validated against exact solutions of Riemann's shock-tube problem, impact of dense media slabs at velocities up to 20 km/sec, flow through a supersonic convergent-divergent nozzle and burning of gunpowder in a closed bomb. Excellent validation results were obtained. The code was then used to predict the performance of two light gas guns (1.5 in. and 0.28 in.) in service at the Ames Research Center. The code predictions were compared with measured pressure histories in the powder chamber and pump tube and with measured piston and projectile velocities. Very good agreement between computational fluid dynamics (CFD) predictions and measurements was obtained. Actual powder-burn rates in the gun were found to be considerably higher (60-90 percent) than predicted by the manufacturer and the behavior of the piston upon yielding appears to differ greatly from that suggested by low-strain rate tests.

Bogdanoff, D. W.↗

A Comparison of Risk Sensitive Path Planning Methods for Aircraft Emergency Landing

Determining the best site to land a damaged aircraft presents some interesting challenges for standard path planning techniques. There are multiple possible locations to consider, the space is 3-dimensional with dynamics, the criteria for a good path is determined by overall risk rather than distance or time, and optimization really matters, since an improved path corresponds to greater expected survival rate. We have investigated a number of different path planning methods for solving this problem, including cell decomposition, visibility graphs, probabilistic road maps (PRMs), and local search techniques. In their pure form, none of these techniques have proven to be entirely satisfactory - some are too slow or unpredictable, some produce highly non-optimal paths or do not find certain types of paths, and some do not cope well with the dynamic constraints when controllability is limited. In the end, we are converging towards a hybrid technique that involves seeding a roadmap with a layered visibility graph, using PRM to extend that roadmap, and using local search to further optimize the resulting paths. We describe the techniques we have investigated, report on our experiments with these techniques, and discuss when and why various techniques were unsatisfactory.

Meuleau, Nicolas↗

Characterization, Operation and Analysis of Test Motors Containing Aluminized Hybrid Fuels

NASA Marshall Space Flight Center's Materials and Processes Department, with support from the Propulsion Systems Department, has renewed the development and maintenance of a hybrid test bed for exposing ablative thermal protection materials to an environment similar to that seen in solid rocket motors (SRM). The Solid Fuel Torch (SFT), operated during the Space Shuttle program, utilized gaseous oxygen for oxidizer and an aluminized hydroxyl-terminated polybutadiene (HTPB) fuel grain to expose a converging section of phenolic material to a 400 psi, 2-phase flow combustion environment. The configuration allows for up to a 2 foot long, 5 inch diameter fuel grain cartridge. Wanting to now test rubber insulation materials with a turn-back feature to mimic the geometry of an aft dome being impinged by alumina particles, the throat area has now been increased by several times to afford flow similarity. Combined with the desire to maintain a higher operating pressure, the oxidizer flow rate is being increased by a factor of 10. Out of these changes has arisen the need to characterize the fuel/oxidizer combination in a higher mass flux condition than has been previously tested at MSFC, and at which the literature has little to no reporting as well. Testing for fuel regression rate comprised a two-level, full factorial design available over Aluminum loading level, mass flow rate, pressure, and diameter. The data taken significantly surpasses the previous available data on regression rate of aluminized HTPB fuel burning with gaseous oxygen. It encompasses higher mass fluxes, and appears to generate more consistent data. The good test article and facility design and testing work of the Penn State HPCL combined with careful analysis of the data and good planning has made this possible. This should be able to assist with developing rate laws that are useful both for research planning and for developing flight system sizing relationships that can help optimize hybrid rocket concepts for trade studies. The successful approach of this DOE and test setup is applicable to other propellant combinations as well.

Kibbey, Timothy P.↗

Characterization, Operation and Analysis of Test Motors Containing Aluminized Hybrid Fuels

NASA Marshall Space Flight Center's Materials and Processes Department, with support from the Propulsion Systems Department, has renewed the development and maintenance of a hybrid test bed for exposing ablative thermal protection materials to an environment similar to that seen in solid rocket motors (SRM). The Solid Fuel Torch (SFT), operated during the Space Shuttle program, utilized gaseous oxygen for oxidizer and an aluminized hydroxyl-terminated polybutadiene (HTPB) fuel grain to expose a converging section of phenolic material to a 400 psi, 2-phase flow combustion environment. The configuration allows for up to a 2 foot long, 5 inch diameter fuel grain cartridge. Wanting to now test rubber insulation materials with a turn-back feature to mimic the geometry of an aft dome being impinged by alumina particles, the throat area has now been increased by several times to afford flow similarity. Combined with the desire to maintain a higher operating pressure, the oxidizer flow rate is being increased by a factor of 10. Out of these changes has arisen the need to characterize the fuel/oxidizer combination in a higher mass flux condition than has been previously tested at MSFC, and at which the literature has little to no reporting as well. Testing for fuel regression rate comprised a two-level, full factorial design available over Aluminum loading level, mass flow rate, pressure, and diameter. The data taken significantly surpasses the previous available data on regression rate of aluminized HTPB fuel burning with gaseous oxygen. It encompasses higher mass fluxes, and appears to generate more consistent data. The good test article and facility design and testing work of the Penn State HPCL combined with careful analysis of the data and good planning has made this possible. This should be able to assist with developing rate laws that are useful both for research planning and for developing flight system sizing relationships that can help optimize hybrid rocket concepts for trade studies. The successful approach of this DOE and test setup is applicable to other propellant combinations as well.

Kibbey, Timothy P.↗

Characterization and Detailed Analysis of Regression Behavior for HTPB Solid Fuels Containing High Aluminum Loadings

NASA Marshall Space Flight Center's Materials and Processes Department, with support from the Propulsion Systems Department, has renewed the development and maintenance of a hybrid test bed for exposing ablative thermal protection materials to an environment similar to that seen in solid rocket motors (SRM). The Solid Fuel Torch (SFT), operated during the Space Shuttle program, utilized gaseous oxygen for oxidizer and an aluminized hydroxyl-terminated polybutadiene (HTPB) fuel grain to expose a converging section of phenolic material to a 400 psi, 2-phase flow combustion environment. The configuration allows for up to a 2 foot long, 5 inch diameter fuel grain cartridge. Wanting to now test rubber insulation materials with a turn-back feature to mimic the geometry of an aft dome being impinged by alumina particles, the throat area has now been increased by several times to afford flow similarity. Combined with the desire to maintain a higher operating pressure, the oxidizer flow rate is being increased by a factor of 10. Out of these changes has arisen the need to characterize the fuel/oxidizer combination in a higher mass flux condition than has been previously tested at MSFC, and at which the literature has little to no reporting as well. For (especially) metalized fuels, hybrid references have pointed out possible dependence of fuel regression rate on a number of variables: mass flux, G - oxidizer only (G0), or - total mass flux (Gtot), Length, L, Pressure, P, and Diameter, D.

Kibbey, Timothy P.↗

Sensitivity of Simulations of Double-detonation Type Ia Supernovae to Integration Methodology

Abstract We study the coupling of hydrodynamics and reactions in simulations of the double-detonation model for Type Ia supernovae. When assessing the convergence of simulations, the focus is usually on spatial resolution; however, the method of coupling the physics together as well as the tolerances used in integrating a reaction network also play an important role. In this paper, we explore how the choices made in both coupling and integrating the reaction portion of a simulation (operator/Strang splitting versus the simplified spectral deferred corrections method we introduced previously) influences the accuracy, efficiency, and nucleosynthesis of simulations of double detonations. We find no need to limit reaction rates or reduce the simulation time step to the reaction timescale. The entire simulation methodology used here is GPU-accelerated and made freely available as part of the Castro simulation code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

The effect of thicker oceanic crust in the Archaean on the growth of continental crust through time

Present crustal evolution models fail to account for the generation of the large volume of continental crust in the required time intervals. All Archaean plate tectonic models, whether invoking faster spreading rates, similar to today's spreading rates, or longer ridge lengths, essentially propose that continental crust has grown by island arc accretion due to the subduction of oceanic crust. The petrological differences that characterize the Archaean from later terrains result from the subduction of hotter oceanic crust into a hotter mantle. If the oceanic crust was appreciably thicker in the Archaean, as geothermal models would indicate, this thicker crust is surely going to have an effect on tectonic processes. A more valid approach is to compare the possible styles of convergence of thick oceanic crust with modern convergence zones. The best modern analog occurs where thick continental crust is colliding with thick continental crust. Oceanic crustal collision on the scale of the present-day Himalayan continental collision zone may have been a frequent occurrence in the Archaean, resulting in extensive partial melting of the hydrous underthrust oceanic crust to produce voluminous tonalite melts, leaving a depleted stabilized basic residuum. Present-day island arc accretion may not have been the dominant mechanism for the growth of the early Archaean crust.

Wilks, M. E.↗

Flow measurements in a model of the mildly curved femoral artery of man

The effects of curvature on the flow rate near the wall in the vicinity of the mildly curved femoral artery of man, and on the pressure distributions along the curved segment, were investigated using glass and tygon flow models constructed to conform to the shape of the femoral angiogram of a human subject. The test fluid was 33 percent aqueous sucrose. Steady flow observations, made using a dye flow visualization system, revealed a flow pattern like that observed in coiled pipes. A double helical type flow was found to develop, with converging streamlines in the wall vicinity from the upper and lower plane of curvature merging asymptotically along the inner curvature in a stable manner. Pressure measurements for steady flow revealed progressively larger pressure drops with distance along the entrance region of the curved segment, relative to that for a straight lumen.

Back, L. H.↗

Numerical Solution of the Steady-State Network Flow Equations for a Non-Ideal Gas

Herein we formulate a steady-state network flow problem for non-ideal gas that relates injection rates and nodal pressures in the network to flows in pipes. For this problem, we present and prove a theorem on uniqueness of generalized solution for a broad class of non-ideal pressure-density relations that satisfy a monotonicity property. Further, we develop a Newton-Raphson algorithm for numerical solution of the steady-state problem, which is made possible by a systematic non-dimensionalization of the equations. The developed algorithm has been extensively tested on benchmark instances and shown to converge robustly to a generalized solution. Previous results [1]-[4], indicate that the steady-state network flow equations for an ideal gas are difficult to solve by the Newton-Raphson method because of its extreme sensitivity to the initial guess. In contrast, we find that non-dimensionalization of the steady-state problem is key to robust convergence of the Newton-Raphson method. We identify criteria based on the uniqueness of solutions under which the existence of a non-physical generalized solution found by a non-linear solver implies non-existence of a physical solution, i.e., infeasibility of the problem. Finally, we compare pressure and flow solutions based on ideal and non-ideal equations of state to demonstrate the need to apply the latter in practice. The solver developed in this article is open-source and is made available for both the academic and research communities as well as the industry.

97 MATHEMATICS AND COMPUTING↗

Significance of blast wave studies to propulsion.

Brief survey of experimental methods currently used for the study of blast wave phenomena with emphasis on high rate exothermic processes. The experimental techniques have used such devices as divergent test sections in shock or detonation tubes, employment of proper test gases, as in marginal detonations, and a variety of explosion systems from finite source explosion apparatus to devices where virtually point explosions are obtained by local breakdown initiated by means of focused laser irradiation. Other methods used are detonation tubes where pressure waves are generated by accelerating flames or by exothermic reactions developed behind reflected shocks, as well as a variety of converging shock and implosion vessels.

Oppenheim, A. K.↗

Application of program LAURA to thermochemical nonequilibrium flow through a nozzle

Program LAURA (Langley Aerothermodynamic Upwind Relaxation Algorithm) is an upwind-biased, point-implicit relaxation algorithm for obtaining the numerical solution to the governing equations for 3D viscous hypersonic flows in chemical and thermal nonequilibrium. The algorithm is derived using a finite-volume formulation in which the inviscid components of flux across cell walls are described with a modified Roe's averaging and with second-order corrections based on Yee's Symmetric Total Variation Diminishing scheme. The code has been applied to Problem 8.2 of this workshop for the case of thermochemical nonequilibrium flow through a nozzle. Chemical reaction rates are defined with the model of Park (1987). Thermal nonequilibrium is modeled using a two-temperature approximation in which the vibrational energies of all molecules are assumed to be in equilibrium at a single temperature which is generally different from the translational-rotational temperature. Two grids were used to define the flow for the original problem, with a stagnation temperature of 6500 K. A third case with a stagnation temperature of 10,000 K is also presented. The solution domain includes the converging nozzle, subsonic flow domain in which the gas is substantially in thermochemical equilibrium and the diverging nozzle, hypersonic flow domain in which the gas is substantially in thermochemical nonequilibrium.

Gnoffo, Peter A.↗

Actively controlled shaft seals for aerospace applications

This study experimentally investigates an actively controlled mechanical seal for aerospace applications. The seal of interest is a gas seal, which is considerably more compact than previous actively controlled mechanical seals that were developed for industrial use. In a mechanical seal, the radial convergence of the seal interface has a primary effect on the film thickness. Active control of the film thickness is established by controlling the radial convergence of the seal interface with piezoelectric actuator. An actively controlled mechanical seal was initially designed and evaluated using a mathematical model. Based on these results, a seal was fabricated and tested under laboratory conditions. The seal was tested with both helium and air, at rotational speeds up to 3770 rad/sec, and at sealed pressures as high as 1.48 x 10(exp 6) Pa. The seal was operated with both manual control and with a closed-loop control system that used either the leakage rate or face temperature as the feedback. The output of the controller was the voltage applied to the piezoelectric actuator. The seal operated successfully for both short term tests (less than one hour) and for longer term tests (four hours) with a closed-loop control system. The leakage rates were typically 5-15 slm (standard liters per minute), and the face temperatures were generally maintained below 100 C. When leakage rate was used as the feedback signal, the setpoint leakage rate was typically maintained within 1 slm. However, larger deviations occurred during sudden changes in sealed pressure. When face temperature was used as the feedback signal, the setpoint face temperature was generally maintained within 3 C, with larger deviations occurring when the sealed pressure changed suddenly.

Salant, Richard F.↗

Photoevaporation of Clumps in Photodissociation Regions

We present the results of an investigation of the effects of Far Ultraviolet (FUV) radiation (6.0eV < hv < 13.6eV) from hot early type OB stars on clumps in star-forming molecular clouds. Clumps in FUV-illuminated regions (or photodissociation regions or PDRs) undergo external heating and photodissociation as they are exposed to the FUV field, resulting in a loss of cold, molecular lump mass as it is converted to warm atomic gas. The heating, if rapid, creates strong photoevaporative mass flows off the clump surfaces, and drives shocks into the clumps, compressing them to high densities. The clumps lose mass on relatively short timescales. The evolution of an individual clump is found to be sensitive to three dimensionless parameters: Nc0, the ratio of the initial column density of the clump to the column N(0) approx. 10(exp 21) cm(exp -2) of a warm FUV-heated surface region; upsilon, the ratio of the sound speed in the heated surface to that in the cold clump material: and t(FUV)t(c), the ratio of the "turn-on time" t(FUV) of the heating flux on a clump to its initial sound crossing-time t(c). The evolution also depends on whether a confining interclump medium exists, or whether the interclump region has negligible pressure, as is the case for turbulence-generated clumps. In this paper, we use spherical 1-D numerical hydrodynamic models as well as approximate analytical models to study the dependence of clump photoevaporation on the physical parameters of the clump, and to derive the dynamical evolution, mass loss rates and photoevaporative timescales of a clump for a variety of astrophysical situations. Turbulent clumps evolve so that their column densities are equal to a critical value determined by the local FUV field, and typically have short photo evaporation timescales, approx. 10(exp 4-5) years for a 1 M(solar mass) clump in a typical star-forming region (Nc0 = 10, upsilon = 10). Clumps with insufficient magnetic pressure support, and in strong FUV fields may be driven to collapse by the compressional effect of converging shock waves. We also estimate the rocket effect on photoevaporating clumps and find that it is significant only for the smallest clumps, with sizes much less than the extent of the PDR itself. Clumps that are confined by all interclump medium may either get completely photoevaporated, or may preserve a shielded core with a warm, dissociated, protective shell that, absorbs the incident FUV flux. We compare our results with observations of some well studied PDRs: the Orion Bar, M17SW NGC 2023 and the Rosette Nebula. The data are consistent with both interpretations of clump origin. turbulence and pressure confinement, with a slight indication for favouring the turbulent model for clumps over pressure-confined

Gorti, Uma↗

Actively Controlled Shaft Seals for Aerospace Applications

This study experimentally investigates an actively controlled mechanical seal for aerospace applications. The seal of interest is a gas seal, which is considerably more compact than previous actively controlled mechanical seals that were developed for industrial use. In a mechanical seal, the radial convergence of the seal interface has a primary effect on the film thickness. Active control of the film thickness is established by controlling the radial convergence of the seal interface with a piezoelectric actuator. An actively controlled mechanical seal was initially designed and evaluated using a mathematical model. Based on these results, a seal was fabricated and tested under laboratory conditions. The seal was tested with both helium and air, at rotational speeds up to 3770 rad/sec, and at sealed pressures as high as 1.48 x 10(exp 6) Pa. The seal was operated with both manual control and with a closed-loop control system that used either the leakage rate or face temperature as the feedback. The output of the controller was the voltage applied to the piezoelectric actuator. The seal operated successfully for both short term tests (less than one hour) and for longer term tests (four hours) with a closed-loop control system. The leakage rates were typically 5-15 slm (standard liters per minute), and the face temperatures were generally maintained below 100C. When leakage rate was used as the feedback signal, the setpoint leakage rate was typically maintained within 1 slm. However, larger deviations occurred during sudden changes in sealed pressure. When face temperature was used as the feedback signal, the setpoint face temperature was generally maintained within 3 C, with larger deviations occurring when the sealed pressure changes suddenly. the experimental results were compared to the predictions from the mathematical model. The model was successful in predicting the trends in leakage rate that occurred as the balance ratio and sealed pressure changed, although the leakage rates were not quantitatively predicted with a high degree of accuracy. This model could be useful in providing valuable design information for future actively controlled mechanical seals.

Salant, Richard F.↗