Search NASA⌕ Search

SEARCH · Search NASA

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 145 records · Page 8

Life sciences and space research XXIV(4) - Natural and artificial ecosystems; Proceedings of the Topical Meeting of the Interdisciplinary Scientific Commission F (Meetings F10, F11, F1 and F12) of the COSPAR 28th Plenary Meeting, The Hague, Netherlands, June 25-July 6, 1990

The present conference on natural and artificial ecosystems and their application to space research encompasses both in-flight and ground-based issues of recycling and control in regenerative life support, the relationships of productivity and facility design in higher plant growth, life-support systems for manned missions to Mars, and biochemical engineering applications in space. Specific issues addressed include interface problems between material recycling systems and plants, temperature and humidity control on a lunar base, the CELSS Test-Facility Project, achieving closure in plant-growth facilities, and life-support systems for Mars transit. Also addressed are a closed equilibrated biological aquatic system, a simulated Mars outpost in the Antarctica dry valleys, analyses of human kidney-cell populations separated on the space shuttle, and the evolution of a phase-separated gravity-independent bioreactor.

Macelroy, R. D.↗

Reynolds stress closure in jet flows using wave models

Ways of implementing the turbulence closure scheme based on modeling the large scale coherent structures as instability waves were sought. The computational tools necessary to apply this scheme to jets of arbitrary geometry were developed. The model, developed earlier, was extended to the shock structure of supersonic jets of arbitrary geometry and multiple jets. It was found that though the qualititate features of the unsteady flow field could be predicted there were always difficulties with some of the quantitative features. This led to the new formation of the closure scheme. The schemes for computations tools which were developed are efficient and represent the application of the very powerful mathematical tools to the problems of practical significance.

Morris, P. J.↗

ZIP3D: An elastic and elastic-plastic finite-element analysis program for cracked bodies

ZIP3D is an elastic and an elastic-plastic finite element program to analyze cracks in three dimensional solids. The program may also be used to analyze uncracked bodies or multi-body problems involving contacting surfaces. For crack problems, the program has several unique features including the calculation of mixed-mode strain energy release rates using the three dimensional virtual crack closure technique, the calculation of the J integral using the equivalent domain integral method, the capability to extend the crack front under monotonic or cyclic loading, and the capability to close or open the crack surfaces during cyclic loading. The theories behind the various aspects of the program are explained briefly. Line-by-line data preparation is presented. Input data and results for an elastic analysis of a surface crack in a plate and for an elastic-plastic analysis of a single-edge-crack-tension specimen are also presented.

Shivakumar, K. N.↗

Coiled Brine Recovery Assembly (CoBRA): A New Approach to Recovering Water from Wastewater Brines

Brine water recovery represents a current technology gap in water recycling for human spaceflight. The role of a brine processor is to take the concentrated discharge from a primary wastewater processor, called brine, and recover most of the remaining water from it. The current state‐of‐the‐art primary processor is the ISS Urine Processor Assembly (UPA) that currently achieves 70% water recovery. Recent advancements in chemical pretreatments are expected to increase this to 85% in the near future. This is a welcome improvement, yet is still not high enough for deep space transit. Mission architecture studies indicate that at least 95% is necessary for a Mars mission, as an example. Brine water recovery is the technology that bridges the gap between 85% and 95%, and moves life support systems one step closer to full closure of the water loop. Several brine water recovery systems have been proposed for human spaceflight, most of them focused on solving two major problems: operation in a weightless environment, and management and containment of brine residual. Brine residual is the leftover byproduct of the brine recovery process, and is often a viscous, sticky paste, laden with crystallized solid particles. Due to the chemical pretreatments added to wastewater prior to distillation in a primary processor, these residuals are typically toxic, which further complicates matters. Isolation of crewmembers from these hazardous materials is paramount. The Coiled Brine Recovery Assembly (CoBRA) is a recently developed concept from the Johnson Space Center that offers solutions to these challenges. CoBRA is centered on a softgoods evaporator that enables a passive fill with brine, and regeneration by discharging liquid brine residual to a collection bag. This evaporator is meant to be lightweight, which allows it to be discarded along with the accumulated brine solids contained within it. This paper discusses design and development of a first CoBRA prototype, and reports initial test results.

Pensinger, Stuart J.↗

The influence of crack closure and elasto-plastic flow on the bending of a cracked plate

The influence of crack closure and elastoplastic flow on the bending behavior of thin cracked plates is investigated using an incremental elastoplastic plate bending finite element computer program. The finite element program was developed using assumptions consistent with Kirchhoff fourth-order plate theory while the material property treatment permits general isotropic work hardening with local elastic unloading. This technique is applied to the problem of a large centrally through cracked plate subject to remote circular bending. Comparison is drawn between two cases of the bending problem. The first neglects the possibility of crack face interference with bending, and the second includes a kinematic prohibition against the crack face from passing through the symmetry plane. Results are reported which isolate the effects of elastoplastic flow and crack closure.

Jones, D. P.↗

Accuracy Assessment of Response Surface Approximations for Supersonic Turbine Design

There is a growing trend to employ CFD tools to supply the necessary information for design optimization of fluid dynamics components/systems. Such results are prone to uncertainties due to reasons including discretization. errors, incomplete convergence of computational procedures, and errors associated with physical models such as turbulence closures. Based on this type of information, gradient-based optimization algorithms often suffer from the noisy calculations, which can seriously compromise the outcome. Similar problems arise from the experimental measurements. Global optimization techniques, such as those based on the response surface (RS) concept are becoming popular in part because they can overcome some of these barriers. However, there are also fundamental issues related to such global optimization technique such as RS. For example, in high dimensional design spaces, typically only a small number of function evaluations are available due to computational and experimental costs. On the other hand, complex features of the design variables do not allow one to model the global characteristics of the design space with simple quadratic polynomials. Consequently a main challenge is to reduce the size of the region where we fit the RS, or make it more accurate in the regions where the optimum is likely to reside. Response Surface techniques using either polynomials or and Neural Network (NN) methods offer designers alternatives to conduct design optimization. The RS technique employs statistical and numerical techniques to establish the relationship between design variables and objective/constraint functions, typically using polynomials. In this study, we aim at addressing issues related to the following questions: (1) How to identify outliers associated with a given RS representation and improve the RS model via appropriate treatments? (2) How to focus on selected design data so that RS can give better performance in regions critical to design optimization? (3) How to combine NN and polynomial techniques for improving the accuracy of the RS model?

Papila, Nilay↗

Design Study for a Global Magnetospheric Dynamics Mission

A successful design was developed, one with many advantages over the original mission. The time spent in orbit was more evenly spread over the region being investigated. The radiation close was significantly lower and the mission did not rely on gravity assist at the moon and thus did not have to make measurements that far out in the tail. A spacecraft design was developed that keeps interference from the engines to a minimum. The design however was quite specific for four spacecraft. It could not be easily scaled to five spacecraft for example. One problem was discovered that is a concern for all similar missions. Inter- spacecraft communication can determine the spacing of the vehicles easily and to the accuracy required. However, the orientation of the polyhedron with the spacecraft at its vertices is not well known for small separations. Ground station range measurements give the line of sight location well but not the angle around that vector. This is a problem any such mission needs to solve. Neither the navigation teams at Goddard nor at Lewis were willing to attempt to solve this problem. At the completion of the study a report was made to the AGU meeting in San Francisco and a paper published in the volume "Science Closure and Enabling Technologies for Constellation Class Missions". This paper is attached.

Russell, C. T.↗

Analytical and phenomenological studies of rotating turbulence

A framework, which combines mathematical analysis, closure theory, and phenomenological treatment, is developed to study the spectral transfer process and reduction of dimensionality in turbulent flows that are subject to rotation. First, we outline a mathematical procedure that is particularly appropriate for problems with two disparate time scales. The approach which is based on the Green's method leads to the Poincare velocity variables and the Poincare transformation when applied to rotating turbulence. The effects of the rotation are now reflected in the modifications to the convolution of a nonlinear term. The Poincare transformed equations are used to obtain a time-dependent analog of the Taylor-Proudman theorem valid in the asymptotic limit when the non-dimensional parameter mu is identical to Omega(t) approaches infinity (Omega is the rotation rate and t is the time). The 'split' of the energy transfer in both direct and inverse directions is established. Secondly, we apply the Eddy-Damped-Quasinormal-Markovian (EDQNM) closure to the Poincare transformed Euler/Navier-Stokes equations. This closure leads to expressions for the spectral energy transfer. In particular, an unique triple velocity decorrelation time is derived with an explicit dependence on the rotation rate. This provides an important input for applying the phenomenological treatment of Zhou. In order to characterize the relative strength of rotation, another non-dimensional number, a spectral Rossby number, which is defined as the ratio of rotation and turbulence time scales, is introduced. Finally, the energy spectrum and the spectral eddy viscosity are deduced.

Mahalov, Alex↗

Elevated temperature crack growth

A study was performed to examine the applicability of path-independent (P-I) integrals to crack growth problems in hot section components of gas turbine aircraft engines. Alloy 718 was used and the experimental parameters included combined temperature and strain cycling, thermal gradients, elastic-plastic strain levels, and mean strains. A literature review was conducted of proposed P-I integrals, and those capable of analyzing hot section component problems were selected and programmed into the postprocessor of a finite element code. Detailed elastic-plastic finite element analyses were conducted to simulate crack growth and crack closure of the test specimen, and to evaluate the P-I integrals. It was shown that the selected P-I integrals are very effective for predicting crack growth for isothermal conditions.

Kim, K. S.↗

Investigations of turbulent scalar fields using probability density function approach

Scalar fields undergoing random advection have attracted much attention from researchers in both the theoretical and practical sectors. Research interest spans from the study of the small scale structures of turbulent scalar fields to the modeling and simulations of turbulent reacting flows. The probability density function (PDF) method is an effective tool in the study of turbulent scalar fields, especially for those which involve chemical reactions. It has been argued that a one-point, joint PDF approach is the one to choose from among many simulation and closure methods for turbulent combustion and chemically reacting flows based on its practical feasibility in the foreseeable future for multiple reactants. Instead of the multi-point PDF, the joint PDF of a scalar and its gradient which represents the roles of both scalar and scalar diffusion is introduced. A proper closure model for the molecular diffusion term in the PDF equation is investigated. Another direction in this research is to study the mapping closure method that has been recently proposed to deal with the PDF's in turbulent fields. This method seems to have captured the physics correctly when applied to diffusion problems. However, if the turbulent stretching is included, the amplitude mapping has to be supplemented by either adjusting the parameters representing turbulent stretching at each time step or by introducing the coordinate mapping. This technique is still under development and seems to be quite promising. The final objective of this project is to understand some fundamental properties of the turbulent scalar fields and to develop practical numerical schemes that are capable of handling turbulent reacting flows.

Gao, Feng↗

Program Promotion Can Distort Space Systems Engineering and Deny Risk

NASA's spectacular success in the Apollo moon landings was achieved against the odds by an obsessive dedication to reducing the great risk. But risk analysis predicted so many astronaut fatalities that it was thought to be unreasonably pessimistic and potentially damaging to the Apollo program. Risk analysis was discontinued, risk was neglected in space shuttle engineering, and so the space shuttle design was unnecessarily dangerous. Since the Apollo era it has been understood that long human space missions would recycle oxygen and water to avoid the very high launch cost of directly supplying them. The development of recycling systems was justified by the need to increase material closure and reduce launch mass. When it was recognized that increasing closure leads to rapidly diminishing returns, the program goal was changed to reducing launch mass and reliability, cost, and risk were considered irrelevant. Systems engineering and especially the discouraging problems of risk and cost have been deliberately ignored because they detract from program promotion, with unfortunate results. Current human launch system design does account for risk and the result strongly resembles Apollo. Current life support design continues to assume recycling, even though the recent great reduction in launch cost now allows direct supply of oxygen and water with significantly better quality, reliability, cost, and risk.

Advocacy↗

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO). PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Material Response↗

Recent Advancements in the PATO Material Response Code

Introduction: Predicting the complicated multiphysics phenomena during atmospheric entry requires high-fidelity modeling tools to refine estimates of mission risks during entry. To this end, new capabilities are being added to the Porous-material Analysis Toolbox based on OpenFOAM (PATO) [1,2,3]. PATO is an open-source software for Computational Material Response (CMR) of reactive porous materials submitted to high-temperature environments. The objective of this work is to highlight current efforts to add to and improve upon the modeling capabilities of PATO. These include efforts to loosely couple PATO with other discipline specialized codes including hypersonic Computational Fluid Dynamics (CFD), to assess the interaction effects between pyrolysis gas blowing and the boundary layer, and Computational Solid Mechanics (CSM), to address modeling of mechanical erosion. Other refinements include surface phenomena modeling capabilities to address the effects of silicone-based coatings applied to the TPS during flight preparation, and a unified multiphase solver for a mixed porous-material and plain-fluid domain. Coupling CMR with CFD (CMR/CFD): A loose coupling between PATO and the Data Parallel Line Relaxation (DPLR) [4] CFD code has been achieved by making use of a blowing boundary condition at the heatshield surface available in DPLR. Starting with heat flux estimates with no pyrolysis gas blowing at the surface, blowing gases are computed by the CMR and passed to the CFD such that aerothermal properties of the environment can be recomputed for a new CMR computation. This leads to an iterative process which is supplemented with an estimate of the radiative heat flux using the Nonequilibrium air radiation (NEQAIR) [5] program. The entire iterative process is illustrated in Figure 1. This coupling strategy has been utilized in computing the MSL material response. The goal is to compare the coupled CMR/CFD results with material response results obtained using traditional blowing corrections. Coupling CMS with CMR: A mechanical erosion model is currently being implemented in PATO to account for the additional mass removal induced by high shear conditions. The modeling process at each timestep consists of updating the mechanical properties as a function of temperature and computing the stress tensor and displacement fields of the material. Then, a failure criteria model determines the regions in which the stress exceeds the ultimate strength values resulting in mesh movement to account for mass removal. This model allows the material response simulation to compute the recession due to both oxidation and shear-induced erosion. The model is demonstrated by computing material response of sphere-cone arc jet samples. Surface Modeling Capabilities: NuSil, a silicone-based coating, was sprayed onto the MSL and Mars 2020 heatshields to mitigate shedding of phenolic dust. To better understand the effects of the NuSil coating on the material response, a novel model has been implemented in PATO. In this model, the equilibrium of the charred NuSil surface is modeled as pure silica, and a constant offset, inspired by the classical spallation model, is added to the the char blowing rate and wall enthalpy to reproduce HyMETS experimental results. The model has also been used to estimate the 3D material response of the MSL heatshield [6]. Unified Solver: In addition to the iterative loose coupling approach mentioned above, a multiphase unified solver is being developed to couple the environment (plain-fluid phase) and the porous-material phase. The solver is based on the volume averaged conservation of mass, momentum, and energy for the macroscale with closure models which include microscale effects through effective physicochemical properties. The unified solver has been used to compute flow through a porous plug and solve the Beavers and Joseph problem [7]. Since the strong coupling between phases is inherent to this solver, modeling assumptions present in other coupling methods of material response are mitigated. This strategy also makes it feasible to capture the competition between surface and volume ablation in the same computational domain, which is usually not possible with other coupling approaches.

Thermal Protection Systems↗

Fatigue crack propagation in aerospace aluminum alloys

This paper reviews fracture mechanics based, damage tolerant characterizations and predictions of fatigue crack growth in aerospace aluminum alloys. The results of laboratory experimentation and modeling are summarized in the areas of: (1) fatigue crack closure, (2) the wide range crack growth rate response of conventional aluminum alloys, (3) the fatigue behavior of advanced monolithic aluminum alloys and metal matrix composites, (4) the short crack problem, (5) environmental fatigue, and (6) variable amplitude loading. Remaining uncertainties and necessary research are identified. This work provides a foundation for the development of fatigue resistant alloys and composites, next generation life prediction codes for new structural designs and extreme environments, and to counter the problem of aging components.

Gangloff, R. P.↗

Turbulent transport modelling of separating and reattaching shear flows

The improvement of capabilities for computer simulation of turbulent recirculating flows was investigated. Attention has been limited to two dimensional flows and principally to statistically stationary motion. Improvement of turbulence modeling explored the treatment of the near wall sublayer and of the exterior fully turbulent region, working within the framework of turbulence closures requiring the solution of transport equations for the turbulence energy and its dissipation rate. The work on the numerical procedure, based on the Gosman-Pun program TEACH, addressed the problems of incorporating the turbulence model as well as the extension to time dependent flows, the incorporation of a third order approximation of convective transport, and the treatment of non-orthogonal boundaries.

Launder, B. E.↗

A study of hypersonic swept shock wave/turbulent boundary layer interactions using a conical Navier-Stokes code

A computational study has been performed of sharp fin-induced swept shock wave/turbulent boundary layer interactions at low hypersonic Mach numbers. The objective was to determine if results obtained using a conical Navier-Stokes code, particularly the peak heating and pressure, are adequate for engineering predictions. The advantage of the conical approach is that the problem becomes two-dimensional and requires much less computational effort than a fully three-dimensional calculation. In this code, the standard Baldwin-Lomax model is used for turbulent closure and its performance is studied in some detail. To assess the approach interactive flowfields generated by unswept sharp fins at two angles of attack at each of three hyprsonic freestream Mach numbers (5, 6, 11) have been calculated and the results compared with experimental wall pressure and heat transfer data. Although the conical Navier-Stokes/Baldwin-Lomax approach is reasonably successful at Mach numbers up to 5, the performance deteriorates as the Mach number is increased. Nevertheless, the approach could be a valuable tool in preliminary parametric design studies.

Rodi, Patrick E.↗

Fast parallel algorithms that compute transitive closure of a fuzzy relation

The notion of a transitive closure of a fuzzy relation is very useful for clustering in pattern recognition, for fuzzy databases, etc. The original algorithm proposed by L. Zadeh (1971) requires the computation time O(n(sup 4)), where n is the number of elements in the relation. In 1974, J. C. Dunn proposed a O(n(sup 2)) algorithm. Since we must compute n(n-1)/2 different values s(a, b) (a not equal to b) that represent the fuzzy relation, and we need at least one computational step to compute each of these values, we cannot compute all of them in less than O(n(sup 2)) steps. So, Dunn's algorithm is in this sense optimal. For small n, it is ok. However, for big n (e.g., for big databases), it is still a lot, so it would be desirable to decrease the computation time (this problem was formulated by J. Bezdek). Since this decrease cannot be done on a sequential computer, the only way to do it is to use a computer with several processors working in parallel. We show that on a parallel computer, transitive closure can be computed in time O((log(sub 2)(n))2).

Kreinovich, Vladik YA.↗

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