Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Closure problem”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS2013

Fracture Mechanics Analyses for Interface Crack Problems - A Review

Recent developments in fracture mechanics analyses of the interfacial crack problem are reviewed. The intent of the review is to renew the awareness of the oscillatory singularity at the crack tip of a bimaterial interface and the problems that occur when calculating mode mixity using numerical methods such as the finite element method in conjunction with the virtual crack closure technique. Established approaches to overcome the nonconvergence issue of the individual mode strain energy release rates are reviewed. In the recent literature many attempts to overcome the nonconvergence issue have been developed. Among the many approaches found only a few methods hold the promise of providing practical solutions. These are the resin interlayer method, the method that chooses the crack tip element size greater than the oscillation zone, the crack tip element method that is based on plate theory and the crack surface displacement extrapolation method. Each of the methods is validated on a very limited set of simple interface crack problems. However, their utility for a wide range of interfacial crack problems is yet to be established.

Krueger, Ronald↗
NASA NTRS1990

Navier-Stokes computation of compressible turbulent flows with a second order closure, part 1

A second order closure turbulence model for compressible flows is developed and implemented in a 2D Reynolds-averaged Navier-Stokes solver. From the beginning where a kappa-epsilon turbulence model was implemented in the bidiagonal implicit method of MACCORMACK (referred to as the MAC3 code) to the final stage of implementing a full second order closure in the efficient line Gauss-Seidel algorithm, numerous work was done, individually and collectively. Besides the collaboration itself, the final product of this work is a second order closure derived from the Launder, Reece, and Rodi model to account for near wall effects, which has been called FRAME model, which stands for FRench-AMerican-Effort. During the reporting period, two different problems were worked out. The first was to provide Ames researchers with a reliable compressible boundary layer code including a wide collection of turbulence models for quick testing of new terms, both in two equations and in second order closure (LRR and FRAME). The second topic was to complete the implementation of the FRAME model in the MAC5 code. The work related to these two different contributions is reported. dilatation in presence of stron shocks. This work, which has been conducted during a work at the Center for Turbulence Research with Zeman aimed also to cros-check earlier assumptions by Rubesin and Vandromme.

Haminh, Hieu↗
NASA NTRS1991

Advancements in engineering turbulence modeling

Some new developments in two-equation models and second order closure models are presented. Two-equation models (k-epsilon models) have been widely used in computational fluid dynamics (CFD) for engineering problems. Most of low-Reynolds number two-equation models contain some wall-distance damping functions to account for the effect of wall on turbulence. However, this often causes the confusion and difficulties in computing flows with complex geometry and also needs an ad hoc treatment near the separation and reattachment points. A set of modified two-equation models is proposed to remove the aforementioned shortcomings. The calculations using various two-equation models are compared with direct numerical simulations of channel flow and flat boundary layers. Development of a second order closure model is also discussed with emphasis on the modeling of pressure related correlation terms and dissipation rates in the second moment equations. All the existing models poorly predict the normal stresses near the wall and fail to predict the 3-D effect of mean flow on the turbulence (e.g. decrease in the shear stress caused by the cross flow in the boundary layer). The newly developed second order near-wall turbulence model is described and is capable of capturing the near-wall behavior of turbulence as well as the effect of 3-D mean flow on the turbulence.

Shih, T.-H.↗
NASA NTRS1991

Hypersonic Shock/Boundary-Layer Interaction Database

Turbulence modeling is generally recognized as the major problem obstructing further advances in computational fluid dynamics (CFD). A closed solution of the governing Navier-Stokes equations for turbulent flows of practical consequence is still far beyond grasp. At the same time, the simplified models of turbulence which are used to achieve closure of the Navier-Stokes equations are known to be rigorously incorrect. While these models serve a definite purpose, they are inadequate for the general prediction of hypersonic viscous/inviscid interactions, mixing problems, chemical nonequilibria, and a range of other phenomena which must be predicted in order to design a hypersonic vehicle computationally. Due to the complexity of turbulence, useful new turbulence models are synthesized only when great expertise is brought to bear and considerable intellectual energy is expended. Although this process is fundamentally theoretical, crucial guidance may be gained from carefully-executed basic experiments. Following the birth of a new model, its testing and validation once again demand comparisons with data of unimpeachable quality. This report concerns these issues which arise from the experimental aspects of hypersonic modeling and represents the results of the first phase of an effort to develop compressible turbulence models.

Settles, G. S.↗
NASA NTRS1988

Steady-state evoked potentials possibilities for mental-state estimation

The use of the human steady-state evoked potential (SSEP) as a possible measure of mental-state estimation is explored. A method for evoking a visual response to a sum-of-ten sine waves is presented. This approach provides simultaneous multiple frequency measurements of the human EEG to the evoking stimulus in terms of describing functions (gain and phase) and remnant spectra. Ways in which these quantities vary with the addition of performance tasks (manual tracking, grammatical reasoning, and decision making) are presented. Models of the describing function measures can be formulated using systems engineering technology. Relationships between model parameters and performance scores during manual tracking are discussed. Problems of unresponsiveness and lack of repeatability of subject responses are addressed in terms of a need for loop closure of the SSEP. A technique to achieve loop closure using a lock-in amplifier approach is presented. Results of a study designed to test the effectiveness of using feedback to consciously connect humans to their evoked response are presented. Findings indicate that conscious control of EEG is possible. Implications of these results in terms of secondary tasks for mental-state estimation and brain actuated control are addressed.

Junker, Andrew M.↗
NASA NTRS2001

Numerical Simulation of Delamination Growth in Composite Materials

The use of decohesion elements for the simulation of delamination in composite materials is reviewed. The test methods available to measure the interfacial fracture toughness used in the formulation of decohesion elements are described initially. After a brief presentation of the virtual crack closure technique, the technique most widely used to simulate delamination growth, the formulation of interfacial decohesion elements is described. Problems related with decohesion element constitutive equations, mixed-mode crack growth, element numerical integration and solution procedures are discussed. Based on these investigations, it is concluded that the use of interfacial decohesion elements is a promising technique that avoids the need for a pre-existing crack and pre-defined crack paths, and that these elements can be used to simulate both delamination onset and growth.

Camanho, P. P.↗
NASA NTRSDate not supplied

Composite Progressive Delamination Analysis of Thin Ply Buckling

The NASA Advanced Composites Project is a multi-year consortium of government, industry and academia to develop and validate technology that will reduce development and certification of advanced composite materials for aircraft structures. One of its tasks is to simulate and predict composite damage progression accurately through high fidelity finite element models. It included a subtask to verify an approach based on the Virtual Crack Closure Technology (VCCT) as well as Floating Node Method (FNM). This paper discusses predicted damage onset and growth in a thin-ply buckling problem. Sensitivity of the simulations to various initial crack length is studied. The compressive load-displacement curves are plotted and the crack growth critical load is predicted. Overall, the crack progression follows simulation benchmarks well. However, the critical-load overshoot phenomenon is observed in both VCCT and FNM, suggesting further investigation may be worth pursuing.

damage progression↗
NASA NTRS1972

Decay of homogeneous turbulence from a specified state

The homogeneous turbulence problem is formulated by first specifying the multipoint velocity correlations or their spectral equivalents at an initial time. Those quantities, together with the correlation or spectral equations, are then used to calculate initial time derivatives of correlations or spectra. The derivatives in turn are used in time series to calculate the evolution of turbulence quantities with time. When the problem is treated in this way, the correlation equations are closed by the initial specification of the turbulence and no closure assumption is necessary. An exponential series which is an iterative solution of the Navier stokes equations gave much better results than a Taylor power series when used with the limited available initial data. In general, the agreement between theory and experiment was good.

Deissler, R. G.↗
NASA NTRS1993

Parameterization of sub-grid scale convection

The following topics are discussed: an overview of the cumulus parameterization problem; interactions between explicit and implicit processes in mesoscale models; effects of model grid size on the cumulus parameterization problem; parameterizing convective effects on momentum fields in mesoscale models; differences between slantwise and vertical cumulus parameterization; experiments with different closure hypotheses; and coupling cumulus parameterizations to boundary layer, stable cloud, and radiation schemes.

Frank, William↗
NASA NTRS2010

Fatigue Crack Closure Analysis Using Digital Image Correlation

Fatigue crack closure during crack growth testing is analyzed in order to evaluate the critieria of ASTM Standard E647 for measurement of fatigue crack growth rates. Of specific concern is remote closure, which occurs away from the crack tip and is a product of the load history during crack-driving-force-reduction fatigue crack growth testing. Crack closure behavior is characterized using relative displacements determined from a series of high-magnification digital images acquired as the crack is loaded. Changes in the relative displacements of features on opposite sides of the crack are used to generate crack closure data as a function of crack wake position. For the results presented in this paper, remote closure did not affect fatigue crack growth rate measurements when ASTM Standard E647 was strictly followed and only became a problem when testing parameters (e.g., load shed rate, initial crack driving force, etc.) greatly exceeded the guidelines of the accepted standard.

Leser, William P.↗
NASA NTRS1984

Simulations of the stratocumulus-topped boundary layer with a third-order closure model

A third order closure model is proposed by Andre et al. (1982), in which the time rate of change terms, the relaxation and rapid effects for the pressure related terms, and the clipping approximation are included along with the quasi-normal closure, to study turbulence in a cloudy layer which is cooled radiatively from above. A spurious oscillation which is strongest near the inversion occurs. An analysis of the problem shows that the oscillation arises from the mean gradient and buoyancy terms of the triple moment equations; these terms are largest near the cloud top. The oscillation is physical, rather than computational. In nature the oscillation is effectively damped, by a mechanism which apparently is not included in our model. In the stably stratified layer just above the mixed layer top, turbulence can excite gravity waves, whose energy is radiated away. Because the closure assumption for the pressure terms does not take into account the transport of wave energy, the model generates spurious oscillations. Damping of the oscillations is possible by introducing diffusion terms into the triple moment equations. With a large enough choice for the diffusion coefficient, the oscillation is effectively eliminated. The results are quite sensitive to the ad hoc eddy coefficient.

Moeng, C. H.↗
NASA NTRS1983

The mechanics of delamination in fiber-reinforced composite materials. Part 1: Stress singularities and solution structure

The fundamental mechanics of delamination in fiber composite laminates is studied. Mathematical formulation of the problem is based on laminate anisotropic elasticity theory and interlaminar fracture mechanics concepts. Stress singularities and complete solution structures associated with general composite delaminations are determined. For a fully open delamination with traction-free surfaces, oscillatory stress singularities always appear, leading to physically inadmissible field solutions. A refined model is introduced by considering a partially closed delamination with crack surfaces in finite-length contact. Stress singularities associated with a partially closed delamination having frictional crack-surface contact are determined, and are found to be diferent from the inverse square-root one of the frictionless-contact case. In the case of a delamination with very small area of crack closure, a simplified model having a square-root stress singularity is employed by taking the limit of the partially closed delamination. The possible presence of logarithmic-type stress singularity is examined; no logarithmic singularity of any kind is found in the composite delamination problem. Numerical examples of dominant stress singularities are shown for delaminations having crack-tip closure with different frictional coefficients between general (1) and (2) graphite-epoxy composites.

Wang, S. S.↗
NASA NTRS1983

The mechanics of delamination in fiber-reinforced composite materials. I - Stress singularities and solution structure

The fundamental mechanics of delamination in fiber composite laminates is studied. Mathematical formulation of the problem is based on laminate anisotropic elasticity theory and interlaminar fracture mechanics concepts. Stress singularities and complete solution structures associated with general composite delaminations are determined. For a fully open delamination with traction-free surfaces, oscillatory stress singularities always appear, leading to physically inadmissible field solutions. A refined model is introduced by considering a partially closed delamination with crack surfaces in finite-length contact. Stress singularities associated with a partially closed delamination having frictional crack-surface contact are determined, and are found to be different from the inverse square-root one of the frictionless-contact case. In the case of a delamination with very small area of crack closure, a simplified model having a square-root stress singularity is employed by taking the limit of the partially closed delamination. The possible presence of logarithmic-type stress singularity is examined; no logarithmic singularity of any kind is found in the composite delamination problem. Numerical examples of dominant stress singularities are shown for delaminations having crack-tip closure with different frictional coefficients between general (1) and (2) graphite-epoxy composites. Previously announced in STAR as N84-13221

Wang, S. S.↗
NASA NTRS2011

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.↗
NASA NTRS1986

Delamination growth in composite materials

The Double Cantilever Beam (DCB) and the End Notched Flexure (ENF) specimens are employed to characterize MODE I and MODE II interlaminar fracture resistance of graphite/epoxy (CYCOM 982) and graphite/PEEK (APC2) composites. Sizing of test specimen geometries to achieve crack growth in the linear elastic regime is presented. Data reduction schemes based upon beam theory are derived for the ENF specimen and include the effects of shear deformation and friction between crack surfaces on compliance, C, and strain energy release rate, G sub II. Finite element (FE) analyses of the ENF geometry including the contact problem with friction are presented to assess the accuracy of beam theory expressions for C and G sub II. Virtual crack closure techniques verify that the ENF specimen is a pure Mode II test. Beam theory expressions are shown to be conservative by 20 to 40 percent for typical unidirectional test specimen geometries. A FE parametric study investigating the influence of delamination length and depth, span, thickness and material properties on G sub II is presented. Mode I and II interlaminar fracture test results are presented. Important experimental parameters are isolated, such as precracking techniques, rate effects, and nonlinear load-deflection response. It is found that subcritical crack growth and inelastic materials behavior, responsible for the observed nonlinearities, are highly rate-dependent phenomena with high rates generally leading to linear elastic response.

Gillespie, J. W., Jr.↗
NASA NTRS1992

Turbulence modeling and experiments

The best way of verifying turbulence is to do a direct comparison between the various terms and their models. The success of this approach depends upon the availability of the data for the exact correlations (both experimental and DNS). The other approach involves numerically solving the differential equations and then comparing the results with the data. The results of such a computation will depend upon the accuracy of all the modeled terms and constants. Because of this it is sometimes difficult to find the cause of a poor performance by a model. However, such a calculation is still meaningful in other ways as it shows how a complete Reynolds stress model performs. Thirteen homogeneous flows are numerically computed using the second order closure models. We concentrate only on those models which use a linear (or quasi-linear) model for the rapid term. This, therefore, includes the Launder, Reece and Rodi (LRR) model; the isotropization of production (IP) model; and the Speziale, Sarkar, and Gatski (SSG) model. Which of the three models performs better is examined along with what are their weaknesses, if any. The other work reported deal with the experimental balances of the second moment equations for a buoyant plume. Despite the tremendous amount of activity toward the second order closure modeling of turbulence, very little experimental information is available about the budgets of the second moment equations. Part of the problem stems from our inability to measure the pressure correlations. However, if everything else appearing in these equations is known from the experiment, pressure correlations can be obtained as the closing terms. This is the closest we can come to in obtaining these terms from experiment, and despite the measurement errors which might be present in such balances, the resulting information will be extremely useful for the turbulence modelers. The purpose of this part of the work was to provide such balances of the Reynolds stress and heat flux equations for the buoyant plume.

Shabbir, Aamir↗
NASA NTRS1993

Uncertainty analysis of diffuse-gray radiation enclosure problems: A hypersensitive case study

An uncertainty analysis of diffuse-gray enclosure problems is presented. The genesis was a diffuse-gray enclosure problem which proved to be hypersensitive to the specification of view factors. This genesis is discussed in some detail. The uncertainty analysis is presented for the general diffuse-gray enclosure problem and applied to the hypersensitive case study. It was found that the hypersensitivity could be greatly reduced by enforcing both closure and reciprocity for the view factors. The effects of uncertainties in the surface emissivities and temperatures are also investigated.

Taylor, Robert P.↗
NASA NTRS2002

Solid-State Compressor for Space Station Oxygen Recovery

At present, the life support system on the International Space Station Alpha vents overboard the carbon dioxide (CO2) produced by the crew members. Recovering the oxygen contained in the CO2 has the potential to reduce resupply mass by 2000 pounds per year or more, a significant weight that could be used for experimental payloads and other valuable items. The technologies used to remove CO2 from the air and to recover O2 from CO2 are flight-ready; however, the interface between the devices is a problem for the Space Station system. Ames Research Center has developed a new technology that solves the interface issue, possibly allowing closure of the oxygen loop in a spacecraft for the first time. CO2 produced by the crew is removed in the Carbon Dioxide Removal Assembly (CDRA). This device effectively produces a pure CO2 stream, but at a very low pressure. Elsewhere, the oxygen generation system which makes O2 by electrolyzing water produces a hydrogen stream. In principle the CO2 and H2 can react to form methane and water over a suitable catalyst. Water produced in this methane-formation reactor can be returned to the water electrolyzer, where the O2 can be returned to the cabin; however, the methane-formation reactor requires CO2 at a much higher pressure than that produced by the CDRA. Furthermore, the CO2 and H2 are often not available at the same time, due to power management and scheduling on the space station. In order to get the CO2 to the reactor at the right pressure and at the right time, a device or assembly that functions as a vacuum pump, compressor, and storage tank is required.

Finn, John E.↗