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 865 records · Page 48

Waterhammer Modeling for the Ares I Upper Stage Reaction Control System Cold Flow Development Test Article

The Upper Stage Reaction Control System provides three-axis attitude control for the Ares I launch vehicle during active Upper Stage flight. The system design must accommodate rapid thruster firing to maintain the proper launch trajectory and thus allow for the possibility to pulse multiple thrusters simultaneously. Rapid thruster valve closure creates an increase in static pressure, known as waterhammer, which propagates throughout the propellant system at pressures exceeding nominal design values. A series of development tests conducted in the fall of 2009 at Marshall Space Flight Center were performed using a water-flow test article to better understand fluid performance characteristics of the Upper Stage Reaction Control System. A subset of the tests examined waterhammer along with the subsequent pressure and frequency response in the flight-representative system and provided data to anchor numerical models. This thesis presents a comparison of waterhammer test results with numerical model and analytical results. An overview of the flight system, test article, modeling and analysis are also provided.

Williams, Jonathan H.↗

Evaluation of Water and Energy Cycles in Contemporary Atmospheric Reanalyses in US Basins

Global atmospheric reanalyses have become an important tool for climate studies, owing to the inclusion of observation driven model quantities that are not easily observed, if at all. This includes the study of regional water and energy budgets, despite the fact that the physical terms of the budgets do not balance, owing in large part to the influence of assimilated observations. Here, we focus on the regional water and energy budgets for key river basins in the US, related to the Humans and Hydroclimate in the US (H2US, a Global Energy and Watercycle EXchanges, GEWEX, regional hydroclimate project). Regional budgets and their residuals (related to the assimilated observations) from contemporary reanalyses will be evaluated and intercompared to assess the improvements over previous reanalyses and identify their strengths and weaknesses as contributions to H2US. When available, the observational analysis increments will be quantified directly, and as residuals when that information is missing. Early results show improvements over the Mississippi River basin water cycle closure, and reduced analysis increments. However, the energy analysis increments show non negligible seasonality. Results for the Colorado River basin and National Climate Assessment (NCA) regions will also be presented. The results will emphasize the need for reanalyses and models to produce complete budget equation terms.

Michael Bosilovich↗

Navier-Stokes turbine heat transfer predictions using two-equation turbulence closures

Navier-Stokes calculations were carried out in order to predict the heat-transfer rates on turbine blades. The calculations were performed using TRAF2D which is a k-epsilon, explicit, finite volume mass-averaged Navier-Stokes solver. Turbulence was modeled using Coakley's q-omega and Chien's k-epsilon two-equation models and the Baldwin-Lomax algebraic model. The model equations along with the flow equations were solved explicitly on a nonperiodic C grid. Implicit residual smoothing (IRS) or a combination of multigrid technique and IRS was applied to enhance convergence rates. Calculations were performed to predict the Stanton number distributions on the first stage vane and blade row as well as the second stage vane row of the SSME high-pressure fuel turbine. The comparison serves to highlight the weaknesses of the turbulence models for use in turbomachinery heat-transfer calculations.

Ameri, Ali A.↗

Theory of inhomogeneous fluids

A theory of nonuniform liquids is presented which is based on the Yvon-Born-Green equation for the one-particle density and the Ornstein-Zernike equation of an inhomogeneous system. The necessary closure is effected by exploiting the solution of a modified hypernetted-chain equation and making a local approximation on the (highly universal) bridge function. The theory is successfully applied to model problems that have also been studied by direct-simulation methods. A brief generalization to quantum liquids is also given.

Nieminen, R. M.↗

The compositional structure and the effects of exothermicity in an unpremixed planar jet flame

Results are presented of direct numerical simulations (DNS) of a randomly perturbed compressible, spatially developing planar jet under the influence of a finite rate Arrhenius chemical reaction of the type F + O yields Product + Heat. The objectives of the simulations are to assess the compositional structure of the flame and to determine the influence of reaction exothermicity by means of statistical sampling of the data generated by DNS. It is shown that even with this idealized kinetics model, the simulated results exhibit features in accordance with experimental data. These results indicate that the Damkoehler number is an important parameter in determining the statistical composition of the reacting field, and that the results are insensitive to the mechanism by which this parameter is varied. It is demonstrated that as the intensity of mixing is increased and the effect of finite rate chemistry is more pronounced, the magnitudes of the ensemble mean and rms of the product mass fraction decrease, and those of the reactants mass fraction increase. Also, at higher mixing rates the joint probability density functions of the reactants' mass fractions shift towards higher values within the composition domain indicating a lower reactedness. These trends are consistent with those observed experimentally and are very useful in portraying the statistical structure of non-equilibrium diffusion flames. The DNS generated data are also utilized to examine the applicability of the 'laminar diffusion flamelet model' in predicting the rate of the reactant conversion with finite rate chemistry. This examination indicates that the bounds of the product formation scale reasonably well with those obtained by the flamelet model. However, in the range of the relatively low values of the Damkoehler numbers considered, the scatter of the results is substantially more than that to be completely determined by the closure.

Steinberger, C. J.↗

Inclusive inelastic scattering of heavy ions in the independent particle model

We consider the inclusive inelastic scattering of heavy ions using the Glauber (1959) model and the independent particle approximation. Inclusive inelastic distributions for projectile excitation of the target and total inelastic scattering, where all projectile and target excited states are summed, are discussed using closure. The total inelastic distribution, when integrated, is shown to be equivalent to the absorption cross section, found from applying the optical theorem to the elastic scattering amplitude in the coherent approximation. Calculations are presented for several heavy-ion pairs, using realistic nuclear densities in a large mass number approximation.

Cucinotta, Francis A.↗

Turbulence modeling for separated flow

Two projects are described in this report. The first involves assessing turbulence models in separated flow. The second addresses the anomalous behavior of certain turbulence models in stagnation point flow. The primary motivation for developing turbulent transport models is to provide tools for computing non-equilibrium, or complex, turbulent flows. Simple flows can be analyzed using data correlations or algebraic eddy viscosities, but in more complicated flows such as a massively separated boundary layer, a more elaborate level of modeling is required. It is widely believed that at least a two-equation transport model is required in such cases. The transport equations determine the evolution of suitable velocity and time-scales of the turbulence. The present study included assessment of second-moment closures in several separated flows, including sharp edge separation; smooth wall, pressure driven separation; and unsteady vortex shedding. Flows with mean swirl are of interest for their role in enhancing mixing both by turbulent and mean motion. The swirl can have a stabilizing effect on the turbulence. An axi-symmetric extension to the INS-2D computer program was written adding the capability of computing swirling flow. High swirl can produce vortex breakdown on the centerline of the jet and it occurs in various combustors.

Durbin, Paul A.↗

A Semi-Analytical Method for Determining the Energy Release Rate of Cracks in Adhesively-Bonded Single-Lap Composite Joints

A semi-analytical method for determining the strain energy release rate due to a prescribed interface crack in an adhesively-bonded, single-lap composite joint subjected to axial tension is presented. The field equations in terms of displacements within the joint are formulated by using first-order shear deformable, laminated plate theory together with kinematic relations and force equilibrium conditions. The stress distributions for the adherends and adhesive are determined after the appropriate boundary and loading conditions are applied and the equations for the field displacements are solved. Based on the adhesive stress distributions, the forces at the crack tip are obtained and the strain energy release rate of the crack is determined by using the virtual crack closure technique (VCCT). Additionally, the test specimen geometry from both the ASTM D3165 and D1002 test standards are utilized during the derivation of the field equations in order to correlate analytical models with future test results. The system of second-order differential field equations is solved to provide the adherend and adhesive stress response using the symbolic computation tool, Maple 9. Finite element analyses using J-integral as well as VCCT were performed to verify the developed analytical model. The finite element analyses were conducted using the commercial finite element analysis software ABAQUS. The results determined using the analytical method correlated well with the results from the finite element analyses.

Yang, Charles↗

Development of Benchmark Examples for Delamination Onset and Fatigue Growth Prediction

An approach for assessing the delamination propagation and growth capabilities in commercial finite element codes was developed and demonstrated for the Virtual Crack Closure Technique (VCCT) implementations in ABAQUS. The Double Cantilever Beam (DCB) specimen was chosen as an example. First, benchmark results to assess delamination propagation capabilities under static loading were created using models simulating specimens with different delamination lengths. For each delamination length modeled, the load and displacement at the load point were monitored. The mixed-mode strain energy release rate components were calculated along the delamination front across the width of the specimen. A failure index was calculated by correlating the results with the mixed-mode failure criterion of the graphite/epoxy material. The calculated critical loads and critical displacements for delamination onset for each delamination length modeled were used as a benchmark. The load/displacement relationship computed during automatic propagation should closely match the benchmark case. Second, starting from an initially straight front, the delamination was allowed to propagate based on the algorithms implemented in the commercial finite element software. The load-displacement relationship obtained from the propagation analysis results and the benchmark results were compared. Good agreements could be achieved by selecting the appropriate input parameters, which were determined in an iterative procedure.

Krueger, Ronald↗

TPSAS-NF1676L-13122-DND

An approach for assessing the delamination propagation and growth capabilities in commercial finite element codes was developed and demonstrated for the Virtual Crack Closure Technique (VCCT) implementations in ABAQUS?. The Double Cantilever Beam (DCB) specimen was chosen as an example. First, benchmark results to assess delamination propagation capabilities under static loading were created using models simulating specimens with different delamination lengths. Second, starting from an initially straight front, the delamination was allowed to propagate based on the algorithms implemented in the commercial finite element software. The load-displacement relationship obtained from the propagation analysis results and the benchmark results were compared. Good agreements could be achieved by selecting the appropriate input parameters, which were determined in an iterative procedure. A benchmark case to assess the fatigue growth prediction capabilities was also created. Input control parameters were varied to study the effect on the computed delamination growth. The benchmarks enabled the selection of the appropriate input parameters that yielded good agreements between the results obtained from the automated analysis and the benchmark results. Overall, the results are encouraging but further assessment for mixed-mode delamination is required.

Ronald Krueger↗

A comparative study of some mathematical models of the mean wind structure and aerodynamic drag of plant canopies

A semianalytical method for describing the mean wind profile and shear stress within plant canopies and for estimating the roughness length and the displacement height is presented. This method incorporates density and vertical structure of the canopy and includes simple parameterizations of the roughness sublayer and shelter factor. Some of the wind profiles examined are consistent with first-order closure techniques while others are consistent with second-order closure techniques. Some profiles show a shearless region near the base of the canopy; however, none displays a secondary maximum there. Comparing several different analytical expressions for the canopy wind profile against observations suggests that one particular type of profile (an Airy function which is associated with the triangular foliage surface area density distribution) is superior to the others. Because of the numerical simplicity of the methods outlined, it is suggested that they may be profitably used in large-scale models of plant-atmosphere exchanges.

Massman, William↗

Separation of crack extension modes in orthotropic delamination models

In the analysis of an interface crack between dissimilar elastic materials, the mode of crack extension is typically not unique, due to oscillatory behavior of near-tip stresses and displacements. This behavior currently limits the applicability of interfacial fracture mechanics as a means to predict composite delamination. The Virtual Crack Closure Technique (VCCT) is a method used to extract mode 1 and mode 2 energy release rates from numerical fracture solutions. The mode of crack extension extracted from an oscillatory solution using the VCCT is not unique due to the dependence of mode on the virtual crack extension length, Delta. In this work, a method is presented for using the VCCT to extract Delta-independent crack extension modes for the case of an interface crack between two in-plane orthotropic materials. The method does not involve altering the analysis to eliminate its oscillatory behavior. Instead, it is argued that physically reasonable, Delta-independent modes of crack extension can be extracted from oscillatory solutions. Knowledge of near-tip fields is used to determine the explicit Delta dependence of energy release rate parameters. Energy release rates are then defined that are separated from the oscillatory dependence on Delta. A modified VCCT using these energy release rate definitions is applied to results from finite element analyses, showing that Delta-independent modes of crack extension result. The modified technique has potential as a consistent method for extracting crack extension modes from numerical solutions. The Delta-independent modes extracted using this technique can also serve as guides for testing the convergence of finite element models. Direct applications of this work include the analysis of planar composite delamination problems, where plies or debonded laminates are modeled as in-plane orthotropic materials.

Beuth, Jack L.↗

Conjecture on superrotation in planetary atmospheres - A diffusion model with mixing length theory

The component of rigid shell superrotation on Venus is discussed in the context of comparative planetary atmospheres. A simplified, heuristic analysis, utilizing mixed length theory to describe the small scale nonlinear advections of energy and angular momentum, thereby providing a closure of the dynamic system, is presented, on the basis of which a crude estimate of zonal velocity is made. The rigid shell (global average) component on Venus was calculated to be 105 m/sec.

Mayr, H. G.↗

Improved Earth Oblateness Rate Reveals Increased Ice Sheet Losses and Mass-Driven Sea Level Rise

Satellite laser ranging (SLR) observations are routinely applied toward the estimation of dynamic oblateness, C(sub 20), which is the largest globally integrated component of Earth's time-variable gravity field. Since 2002, GRACE and GRACE Follow-On have revolutionized the recovery of higher spatial resolution features of global time-variable gravity, with SLR continuing to provide the most reliable estimates of C (sub 20).We quantify the effect of various SLR processing strategies on estimating C(sub 20) and demonstrate better signal recovery with the inclusion of GRACE-derived low-degree gravity information in the forward model. This improved SLR product modifies the Antarctic and Greenland Ice Sheet mass trends by -15.4 and -3.5 Gt/year, respectively, as compared to CSR TN11, and improves global mean sea level budget closure by modifying sea level rise by +0.08 mm/year. We recommend that this new C(sub 20) product be applied to RL06 GRACE data products for enhanced accuracy and scientific interpretation.

Loomis, B. D.↗

Turbulence energy and diffusion transport in a separating and reattaching flow

For accurate prediction of the turbulent flow in separated and reattaching regions, it is necessary to incorporate the second- and third-moments of turbulent fluctuations. The turbulence energy and the energy dissipation rate equations are modified by incorporating the second-order closure. Moreover, the third-order closure with near-wall correction is developed for the evaluation of the diffusive action of the second-moments. After comparison of the results with experimental data, it is shown that the models developed here improve the prediction of triple-velocity correlations in both recirculating and developing flow regions.

Amano, R. S.↗

Environmental fatigue in aluminum-lithium alloys

Aluminum-lithium alloys exhibit similar environmental fatigue crack growth characteristics compared to conventional 2000 series alloys and are more resistant to environmental fatigue compared to 7000 series alloys. The superior fatigue crack growth behavior of Al-Li alloys 2090, 2091, 8090, and 8091 is due to crack closure caused by tortuous crack path morphology and crack surface corrosion products. At high R and reduced closure, chemical environment effects are pronounced resulting in accelerated near threshold da/dN. The beneficial effects of crack closure are minimized for small cracks resulting in rapid growth rates. Limited data suggest that the 'chemically small crack' effect, observed in other alloy system, is not pronounced in Al-Li alloys. Modeling of environmental fatigue in Al-Li-Cu alloys related accelerated fatigue crack growth in moist air and salt water to hydrogen embrittlement.

Piascik, Robert S.↗

A long-term space astrophysics research program: The evolution of the quasar continuum

The research program supported by this grant now has great momentum. Numerous papers are in progress, and a strong multi-wavelength observing program is rapidly accumulating data on samples of high redshift quasars across the spectrum. ROSAT spectra of quasars continue to yield surprises. Of four z = 3 quasars with X-ray spectra, three show strong absorption. This contrasts strongly with the situation for luminous AGN at low redshifts where fewer than 1 in 20 show X-ray absorption. A new site for this absorption is probably needed, either around the quasar (e.g. in a cluster cooling flow) or along the line of sight (e.g. in a Damped Lyman-alpha system). The unabsorbed quasar allows limits on the physical conditions in a damped Lyman-alpha cloud to be calculated, and will allow a X-ray Gunn-Peterson test to be applied that will limit the fraction of the closure mass in an intergalactic medium. The X-ray spectral indices of these z = 3 quasars show no change from those of similar objects at low z, suggesting that 'short-lifetime' models apply. Eight other z = 3-4 quasars have been detected and their energy distributions from X-rays to Infrared (using new infrared spectrographs) have been compiled. These are now being compared with the low z continua from the 'Atlas of Quasar Energy Distributions' to search for evolutionary changes. The discovery of a likely warm absorber in 3C351 made recognition of another example simple. Also, modeling of the conditions in the absorber in 3C351 using the OVI absorption line from HST and the high ionization emission lines, suggests that the broad line region is indeed the origin of the warm absorber in this quasar, and by extension, others. Warm absorbers can now be used as a new diagnostic of this region. The X-ray spectrum of a 'Red Quasar', 3C212, has a cut-off spectrum, which could be fitted by an absorbed power-law, or more remarkably, by an unabsorbed black body. Using our quasi-simultaneous optical data and photoionization modeling we are able to rule out this latter possibility. The long-sought key to understanding red quasars has now been found. (Reddening was earlier ruled out because no 2175A absorption feature was seen. This is now seen to have been misleading.)

Elvis, M.↗

Model-based automatic generation of grasping regions

The problem of automatically generating stable regions for a robotic end effector on a target object, given a model of the end effector and the object is discussed. In order to generate grasping regions, an initial valid grasp transformation from the end effector to the object is obtained based on form closure requirements, and appropriate rotational and translational symmetries are associated with that transformation in order to construct a valid, continuous grasping region. The main result of this algorithm is a list of specific, valid grasp transformations of the end effector to the target object, and the appropriate combinations of translational and rotational symmetries associated with each specific transformation in order to produce a continuous grasp region.

Bloss, David A.↗