Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Closure 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 973 records · Page 54

Constraint Embedding Technique for Multibody System Dynamics

Multibody dynamics play a critical role in simulation testbeds for space missions. There has been a considerable interest in the development of efficient computational algorithms for solving the dynamics of multibody systems. Mass matrix factorization and inversion techniques and the O(N) class of forward dynamics algorithms developed using a spatial operator algebra stand out as important breakthrough on this front. Techniques such as these provide the efficient algorithms and methods for the application and implementation of such multibody dynamics models. However, these methods are limited only to tree-topology multibody systems. Closed-chain topology systems require different techniques that are not as efficient or as broad as those for tree-topology systems. The closed-chain forward dynamics approach consists of treating the closed-chain topology as a tree-topology system subject to additional closure constraints. The resulting forward dynamics solution consists of: (a) ignoring the closure constraints and using the O(N) algorithm to solve for the free unconstrained accelerations for the system; (b) using the tree-topology solution to compute a correction force to enforce the closure constraints; and (c) correcting the unconstrained accelerations with correction accelerations resulting from the correction forces. This constraint-embedding technique shows how to use direct embedding to eliminate local closure-loops in the system and effectively convert the system back to a tree-topology system. At this point, standard tree-topology techniques can be brought to bear on the problem. The approach uses a spatial operator algebra approach to formulating the equations of motion. The operators are block-partitioned around the local body subgroups to convert them into aggregate bodies. Mass matrix operator factorization and inversion techniques are applied to the reformulated tree-topology system. Thus in essence, the new technique allows conversion of a system with closure-constraints into an equivalent tree-topology system, and thus allows one to take advantage of the host of techniques available to the latter class of systems. This technology is highly suitable for the class of multibody systems where the closure-constraints are local, i.e., where they are confined to small groupings of bodies within the system. Important examples of such local closure-constraints are constraints associated with four-bar linkages, geared motors, differential suspensions, etc. One can eliminate these closure-constraints and convert the system into a tree-topology system by embedding the constraints directly into the system dynamics and effectively replacing the body groupings with virtual aggregate bodies. Once eliminated, one can apply the well-known results and algorithms for tree-topology systems to solve the dynamics of such closed-chain system.

Woo, Simon S.↗

Separation of crack extension modes in composite delamination problems

This work concerns fracture mechanics modeling of composite delamination problems. In order to predict delamination resistance, an applied stress intensity factor, K, or energy release rate, G, must be compared to a mode-dependent critical value of K or G from experiment. In the interfacial fracture analysis of most applications and some tests, the mode of crack extension is not uniquely defined. It is instead a function of distance from the crack tip due to the oscillating singularity existing at the tip. In this work, a consistent method is presented of extracting crack extension modes in such cases. In particular, use of the virtual crack closure technique (VCCT) to extract modes of crack extension is studied for cases of a crack along the interface between two in-plane orthotropic materials. Modes of crack extension extracted from oscillatory analyses using VCCT are a function of the virtual crack extension length, delta. Most existing efforts to obtain delta-independent modes of crack extension involve changing the analysis in order to eliminate its oscillatory nature. One such method involves changing one or more properties of the layers to make the oscillatory exponent parameter, epsilon, equal zero. Standardized application of this method would require consistent criteria for identifying which properties can be altered without changing the physical aspects of the problem. Another method involves inserting a thin homogeneous layer (typically referred to as a resin interlayer) along the interface and placing the crack within it. The drawbacks of this method are that it requires increased modeling effort and introduces the thickness of the interlayer as an additional length parameter. The approach presented here does not attempt to alter the interfacial fracture analysis to eliminate its oscillatory behavior. Instead, the argument is made that the oscillatory behavior is non-physical and that if its effects were separated from VCCT quantities, then consistent, delta-independent modes of crack extension could be defined. Knowledge of the near-tip fields in a planar orthotropic material interfacial fracture analysis is used to determine the explicit delta dependence of VCCT parameters. Once this delta dependence is determined, energy release rates are defined with this delta dependence factored out. This modified VCCT method is applied to results from two finite element test cases. It is shown that, as predicted, delta-independent modes of crack extension result. The modified VCCT approach shows potential as a consistent method of extracting crack extension modes. It uses the same information from a finite element analysis (i.e., nodal forces and displacements) as the traditional VCCT method does. The A-independent modes extracted using the modified VCCT approach can also be used as guides to test the convergence of finite element solutions.

Beuth, Jack L.↗

Effect of propellant deformation on ignition and combustion processes in solid propellant cracks

A comprehensive theoretical model was formulated to study the development of convective burning in a solid propellant crack which continually deforms due to burning and pressure loading. In the theoretical model, the effect of interrelated structural deformation and combustion processes was taken into account by considering (1) transient, one dimensional mass, momentum, and energy conservation equations in the gas phase; (2) a transient, one dimensional heat conduction equation in the solid phase; and (3) quasi-static deformation of the two dimensional, linear viscoelastic propellant crack caused by pressure loading. Partial closures may generate substantial local pressure peaks along the crack, implying a strong coupling between chamber pressurization, crack combustion, and propellant deformation, especially when the cracks are narrow and the chamber pressurization rates high. The maximum pressure in the crack cavity is generally higher than that in the chamber. The initial flame-spreading process is not affected by propellant deformation.

Kumar, M.↗

Glass shell manufacturing in space

Highly-uniform, hollow glass spheres (shells), which are used for inertial confinement fusion targets, were formed from metal-organic gel powder feedstock in a vertical furnace. As a result of the rapid pyrolysis caused by the furnace, the gel is transformed to a shell in five distinct stages: (a) surface closure of the porous gel; (b) generation of a closed-cell foam structure in the gel; (c) spheridization of the gel and further expansion of the foam; (d) coalescence of the closed-cell foam to a single-void shell; and (e) fining of the glass shell. The heat transfer from the furnace to the falling gel particle was modeled to determine the effective heating rate of the gel. The model predicts the temperature history for a particle as a function of mass, dimensions, specific heat, and absorptance as well as furnace temperature profile and thermal conductivity of the furnace gas. A model was developed that predicts the gravity-induced degradation of shell concentricity in falling molten shells as a function of shell characteristics and time.

Downs, R. L.↗

Experimental investigation of shock-cell noise reduction for dual-stream nozzles in simulated flight comprehensive data report. Volume 1: Test nozzles and acoustic data

Parameters which contribute to supersonic jet shock noise were investigated for the purpose of determining means to reduce such noise generation to acceptable levels. Six dual-stream test nozzles with varying flow passage and plug closure designs were evaluated under simulated flight conditions in an anechoic chamber. All nozzles had combined convergent-divergent or convergent flow passages. Acoustic behavior as a function of nozzle flow passage geometry was measured. The acoustic data consist primarily of 1/3 octave band sound pressure levels and overall sound pressure levels. Detailed schematics and geometric characteristics of the six scale model nozzle configurations and acoustic test point definitions are presented. Tabulation of aerodynamic test conditions and a computer listing of the measured acoustic data are displayed.

Yamamoto, K.↗

Transonic Navier-Stokes flow computations over wing-fuselage geometries

The steady-state transonic viscous flows over two wing-fuselage configurations are solved numerically. The 3D compressible thin-layer Navier-Stokes equations are solved by means of multistage timestepping similar to the Runge-Kutta method. The timestepping technique is applied with a multigrid scheme to make the transonic Navier-Stokes computations using the Baldwin-Lomax turbulence for closure and converging the solution on a grid of 1.52 million grid points with a C-O topology. Both generic and fighter wing-fuselage geometries are solved, and the results are compared to experimental studies. Solution convergence is found to be fast, and the fine grid generates accurate results for both configurations by means of the Baldwin-Lomax turbulence model. Results from the Johnson-King model calculations produce results that are more consistent with experimental results for the fighter configuration.

Agrawal, S.↗

Large eddy simulations and direct numerical simulations of high speed turbulent reacting flows

The basic objective of this research is to extend the capabilities of Large Eddy Simulations (LES) and Direct Numerical Simulations (DNS) for the computational analyses of high speed reacting flows. In the efforts related to LES, we were primarily involved with assessing the performance of the various modern methods based on the Probability Density Function (PDF) methods for providing closures for treating the subgrid fluctuation correlations of scalar quantities in reacting turbulent flows. In the work on DNS, we concentrated on understanding some of the relevant physics of compressible reacting flows by means of statistical analysis of the data generated by DNS of such flows. In the research conducted in the second year of this program, our efforts focused on the modeling of homogeneous compressible turbulent flows by PDF methods, and on DNS of non-equilibrium reacting high speed mixing layers. Some preliminary work is also in progress on PDF modeling of shear flows, and also on LES of such flows.

Givi, P.↗

Gear crack propagation investigations

Analytical and experimental studies were performed to investigate the effect of gear rim thickness on crack propagation life. The FRANC (FRacture ANalysis Code) computer program was used to simulate crack propagation. The FRANC program used principles of linear elastic fracture mechanics, finite element modeling, and a unique re-meshing scheme to determine crack tip stress distributions, estimate stress intensity factors, and model crack propagation. Various fatigue crack growth models were used to estimate crack propagation life based on the calculated stress intensity factors. Experimental tests were performed in a gear fatigue rig to validate predicted crack propagation results. Test gears were installed with special crack propagation gages in the tooth fillet region to measure bending fatigue crack growth. Good correlation between predicted and measured crack growth was achieved when the fatigue crack closure concept was introduced into the analysis. As the gear rim thickness decreased, the compressive cyclic stress in the gear tooth fillet region increased. This retarded crack growth and increased the number of crack propagation cycles to failure.

Lewicki, David G.↗

Influence of Finite Element Software on Energy Release Rates Computed Using the Virtual Crack Closure Technique

Strain energy release rates were computed along straight delamination fronts of Double Cantilever Beam, End-Notched Flexure and Single Leg Bending specimens using the Virtual Crack Closure Technique (VCCT). Th e results were based on finite element analyses using ABAQUS# and ANSYS# and were calculated from the finite element results using the same post-processing routine to assure a consistent procedure. Mixed-mode strain energy release rates obtained from post-processing finite elem ent results were in good agreement for all element types used and all specimens modeled. Compared to previous studies, the models made of s olid twenty-node hexahedral elements and solid eight-node incompatible mode elements yielded excellent results. For both codes, models made of standard brick elements and elements with reduced integration did not correctly capture the distribution of the energy release rate acr oss the width of the specimens for the models chosen. The results suggested that element types with similar formulation yield matching results independent of the finite element software used. For comparison, m ixed-mode strain energy release rates were also calculated within ABAQUS#/Standard using the VCCT for ABAQUS# add on. For all specimens mod eled, mixed-mode strain energy release rates obtained from ABAQUS# finite element results using post-processing were almost identical to re sults calculated using the VCCT for ABAQUS# add on.

Krueger, Ronald↗

NEWS Climatology Project: The State of the Water Cycle at Continental to Global Scales

NASA's Energy and Water Cycle Study (NEWS) program fosters collaborative research towards improved quantification and prediction of water and energy cycle consequences of climate change. In order to measure change, it is first necessary to describe current conditions. The goal of the NEWS Water and Energy Cycle Climatology project is to develop "state of the global water cycle" and "state of the global energy cycle" assessments based on data from modern ground and space based observing systems and data integrating models. The project is a multiinstitutional collaboration with more than 20 active contributors. This presentation will describe results of the first stage of the water budget analysis, whose goal was to characterize the current state of the water cycle on mean monthly, continental scales. We examine our success in closing the water budget within the expected uncertainty range and the effects of forcing budget closure as a method for refining individual flux estimates.

Rodell, Matthew↗

Numerical Simulations of Homogeneous Turbulence Using Lagrangian-Averaged Navier-Stokes Equations

The Lagrangian-averaged Navier-Stokes (LANS) equations are numerically evaluated as a turbulence closure. They are derived from a novel Lagrangian averaging procedure on the space of all volume-preserving maps and can be viewed as a numerical algorithm which removes the energy content from the small scales (smaller than some a priori fixed spatial scale alpha) using a dispersive rather than dissipative mechanism, thus maintaining the crucial features of the large scale flow. We examine the modeling capabilities of the LANS equations for decaying homogeneous turbulence, ascertain their ability to track the energy spectrum of fully resolved direct numerical simulations (DNS), compare the relative energy decay rates, and compare LANS with well-accepted large eddy simulation (LES) models.

Mohseni, Kamran↗

Integrated Design and Engineering Analysis (IDEA) Environment - Propulsion Related Module Development and Vehicle Integration

This report documents the work performed during the period from May 2011 - October 2012 on the Integrated Design and Engineering Analysis (IDEA) environment. IDEA is a collaborative environment based on an object-oriented, multidisciplinary, distributed framework using the Adaptive Modeling Language (AML). This report will focus on describing the work done in the areas of: (1) Integrating propulsion data (turbines, rockets, and scramjets) in the system, and using the data to perform trajectory analysis; (2) Developing a parametric packaging strategy for a hypersonic air breathing vehicles allowing for tank resizing when multiple fuels and/or oxidizer are part of the configuration; and (3) Vehicle scaling and closure strategies.

Kamhawi, Hilmi N.↗

Association between energetic particle bursts and Birkeland currents in the geomagnetic tail

IMP 8 observations in the geomagnetic tail of energetic particles at hundreds of keV and of the magnetic field are examined for a substorm period from 1700 to 1800 UT on Mar. 3, 1976. For this interval, it is found that bursts of energetic protons detected at a downstream distance of about 37 earth radii in the midnight tail region are accompanied by east-west magnetic field perturbations indicative of Birkeland current sheets. The magnitude of the magnetic field perturbations range from about 1 to 12 nT, corresponding to current densities of 1 to 9 mA/m when integrated over the current sheet thickness. These values are comparable to previous observations of Birkeland currents in the near-earth tail (X greater than -20 earth radii) and with the auroral-arc-associated Birkeland current densities at low altitudes scaled out to the tail region. The favorable comparison suggests that the closure of small-scale Birkeland current systems, at least in these instances, lies tailward of about 40 earth radii. The observed Birkeland current system and the magnitude of current densities are found to be consistent with ion tearing instability in the tail. The inferred wavelength of the tearing mode along the tail axis is about 2 to 17 earth radii. In the context of the tearing model, the presence of Birkeland currents and the occasional isotropic distribution of energetic electrons (of more than 0.22 MeV) suggest that the energetic particle events in this study occur on closed magnetic field lines, even though strong tailward streaming dominates the particle anisotropy in these events.

Lui, A. T. Y.↗

Global irrigation contribution to wheat and maize yield

Irrigation is the largest sector of human water use and an important option for increasing crop production and reducing drought impacts. However, the potential for irrigation to contribute to global crop yields remains uncertain. Here, we quantify this contribution for wheat and maize at global scale by developing a Bayesian framework integrating empirical estimates and gridded global crop models on new maps of the relative difference between attainable rainfed and irrigated yield (ΔY). At global scale, ΔY is 34 ± 9% for wheat and 22 ± 13% for maize, with large spatial differences driven more by patterns of precipitation than that of evaporative demand. Comparing irrigation demands with renewable water supply, we find 30–47% of contemporary rainfed agriculture of wheat and maize cannot achieve yield gap closure utilizing current river discharge, unless more water diversion projects are set in place, putting into question the potential of irrigation to mitigate climate change impacts.

Irrigation↗

Symbolic Computation of Strongly Connected Components Using Saturation

Finding strongly connected components (SCCs) in the state-space of discrete-state models is a critical task in formal verification of LTL and fair CTL properties, but the potentially huge number of reachable states and SCCs constitutes a formidable challenge. This paper is concerned with computing the sets of states in SCCs or terminal SCCs of asynchronous systems. Because of its advantages in many applications, we employ saturation on two previously proposed approaches: the Xie-Beerel algorithm and transitive closure. First, saturation speeds up state-space exploration when computing each SCC in the Xie-Beerel algorithm. Then, our main contribution is a novel algorithm to compute the transitive closure using saturation. Experimental results indicate that our improved algorithms achieve a clear speedup over previous algorithms in some cases. With the help of the new transitive closure computation algorithm, up to 10(exp 150) SCCs can be explored within a few seconds.

Zhao, Yang↗

Exploration Life Support Critical Questions for Future Human Space Missions

Exploration Life Support (ELS) is a current project under NASA's Exploration Systems Mission Directorate. The ELS Project plans, coordinates and implements the development of advanced life support technologies for human exploration missions in space. Recent work has focused on closed loop atmosphere and water systems for long duration missions, including habitats and pressurized rovers. But, what are the critical questions facing life support system developers for these and other future human missions? This paper explores those questions and how progress in the development of ELS technologies can help answer them. The ELS Project includes the following Elements: Atmosphere Revitalization Systems, Water Recovery Systems, Waste Management Systems, Habitation Engineering, Systems Integration, Modeling and Analysis, and Validation and Testing, which includes the Sub-Elements Flight Experiments and Integrated Testing. Systems engineering analysis by ELS seeks to optimize overall mission architectures by considering all the internal and external interfaces of the life support system and the potential for reduction or reuse of commodities. In particular, various sources and sinks of water and oxygen are considered along with the implications on loop closure and the resulting launch mass requirements. Systems analysis will be validated through the data gathered from integrated testing, which will demonstrate the interfaces of a closed loop life support system. By applying a systematic process for defining, sorting and answering critical life support questions, the ELS project is preparing for a variety of future human space missions

Kwert, Michael K.↗

Flight Test 4 Preliminary Results: NASA Ames SSI

Realization of the expected proliferation of Unmanned Aircraft System (UAS) operations in the National Airspace System (NAS) depends on the development and validation of performance standards for UAS Detect and Avoid (DAA) Systems. The RTCA Special Committee 228 is charged with leading the development of draft Minimum Operational Performance Standards (MOPS) for UAS DAA Systems. NASA, as a participating member of RTCA SC-228 is committed to supporting the development and validation of draft requirements as well as the safety substantiation and end-to-end assessment of DAA system performance. The Unmanned Aircraft System (UAS) Integration into the National Airspace System (NAS) Project conducted flight test program, referred to as Flight Test 4, at Armstrong Flight Research Center from April -June 2016. Part of the test flights were dedicated to the NASA Ames-developed Detect and Avoid (DAA) System referred to as JADEM (Java Architecture for DAA Extensibility and Modeling). The encounter scenarios, which involved NASA's Ikhana UAS and a manned intruder aircraft, were designed to collect data on DAA system performance in real-world conditions and uncertainties with four different surveillance sensor systems. Flight test 4 has four objectives: (1) validate DAA requirements in stressing cases that drive MOPS requirements, including: high-speed cooperative intruder, low-speed non-cooperative intruder, high vertical closure rate encounter, and Mode CS-only intruder (i.e. without ADS-B), (2) validate TCASDAA alerting and guidance interoperability concept in the presence of realistic sensor, tracking and navigational errors and in multiple-intruder encounters against both cooperative and non-cooperative intruders, (3) validate Well Clear Recovery guidance in the presence of realistic sensor, tracking and navigational errors, and (4) validate DAA alerting and guidance requirements in the presence of realistic sensor, tracking and navigational errors. The results will be presented at RTCA Special Committee 228 in support of final verification and validation of the DAA MOPS.

Detect-And-Avoid (DAA)↗

Single Field-of-View Sounding Atmospheric Products

The Single Field-of-view Sounder Atmospheric Products (SiFSAP) answer the need for a novel high spatial resolution atmospheric data product for major hyper-spectral infrared (IR) sounder missions. SiFSAP include a complete set of atmospheric vertical profiles, cloud, and surface properties, which are physically retrieved from top-of-atmosphere (TOA) spectral radiances under all-sky conditions via a rigorously defined radiative transfer relationship. By using a state-of-art fast radiative transfer model and a carefully designed optimal estimation based physical retrieval scheme, the SiFSAP algorithm ensures both an ultra-fast data processing speed needed for operational weather applications and the radiometric consistency desired by long-term climate studies. SiFSAP supplement existing operational products by providing data at the native resolution of the sounder instruments, the direct and accurate retrieval of cloud scattering properties, and the establishment of ‘radiance closure.’

Wan Wu↗