Search NASA⌕ Search

SEARCH · Search NASA

Results for “Constitutive models”

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 469 records · Page 26

High-Lift Prediction Workshop 5: Summary of Reynolds-Averaged Navier–Stokes Technology Focus Group Summary

The fifth High-Lift Prediction Workshop (HLPW-5), which involved the high-lift version of the NASA common research model in several configurations, assessed various computational fluid dynamics methods, including Reynolds-averaged Navier-Stokes (RANS) and hybrid large-eddy simulations. This paper summarizes RANS solutions computed on fixed grids. Case 1, a verification case, considered a simple wing-body configuration and focused on grid convergence of lift, drag, and pitching moment coefficients. Case 2 is a configuration buildup case that focused on predicting the effects of increasing geometric complexity. For buildup configurations with slats, flaps, and a nacelle/pylon (configurations 2.2, 2.3, and 2.4), wind-tunnel data were provided by ONERA. Case 3 focused on the reference landing configuration at four Reynolds-number conditions. For Case 1, grid-converged RANS solutions were achieved using the standard Spalart-Allmaras (SA) turbulence model and the SA model with a quadratic constitutive relation and a rotation correction. Agreement between RANS solutions was observed for the simplest configuration 2.1 of Case 2 with the standard SA model. For other configurations, agreement between RANS solutions was hampered by insufficient iterative and grid convergence, especially at high angles of attack. In comparison with the experiment, RANS solutions qualitatively showed the correct configuration buildup trend but underpredicted lift and overpredicted both drag and pitching moment at high angles of attack.

Boris Diskin↗

A Damage Model for the Simulation of Delamination in Advanced Composites under Variable-Mode Loading

A thermodynamically consistent damage model is proposed for the simulation of progressive delamination in composite materials under variable-mode ratio. The model is formulated in the context of Damage Mechanics. A novel constitutive equation is developed to model the initiation and propagation of delamination. A delamination initiation criterion is proposed to assure that the formulation can account for changes in the loading mode in a thermodynamically consistent way. The formulation accounts for crack closure effects to avoid interfacial penetration of two adjacent layers after complete decohesion. The model is implemented in a finite element formulation, and the numerical predictions are compared with experimental results obtained in both composite test specimens and structural components.

Turon, A.↗

A film thickness analysis for line contacts under pure rolling conditions with a non-Newtonian rheological model

The non-Newtonian constitutive equation of Winer and Bair (1979) is applied to a conventional isothermal film thickness analysis of line contact lubrication of rolling elements providing dimensionless film thickness equations for four lubrication methods. The equations can be applied to high viscosity and rolling speed and low limiting shear stress cases where Newtonian methods are not applicable. This analysis is based on a reasonable range of limiting shear stress which is smaller than the Newtonian values for low viscosity lubricants.

Gecim, B.↗

Transient Creep of a Composite Lower Crust: Constitutive Theory - 1

A composite model is proposed to describe the time-dependent response of the Earth's lower crust. The motivation for such it model is twofold: First, new observations of widespread postseismic deformation indicate that the deep continental crust responds viscoelastically, having both long-and short-term decay times. Second, by any number of observationally based rationales, the lower crust is compositionally and structurally heterogeneous over many length scales. For heterogeneities that have much smaller characteristic lengths than the minimum deformation wavelength of interest, the aggregate rheology can be described by composite media theory. For wavelengths of the order of the thickness of the lower crust (approx. = 25-40 km) and larger, composite theory may be applied to heterogeneities that are smaller than about several hundred meters, or equivalent to the vertical extent of a thick lower crustal mylonitic shear zone. The composite media theory developed here is constructed using both Eshelhy-Mori-Tanaka theory for aligned generalized spheroidal inclusions and a generalized self-consistent method. The inclusions and matrix are considered to be Maxwellian viscoelastic: a rheology that is consistent with past homogeneous models of postseismic stress relaxation. The composite theory presented here introduces a transient response to a suddenly imposed stress field which does not appear in homogeneous Maxwell models. Analytic expressions for the amplitude and duration of the transient and for the effective long-and short-term viscosities of the composite are given which describe the sensitivity to inclusion concentration (phi), to shape, and to ratio of inclusion-to-matrix viscosity (R).

Ivins, Erik R.↗

Modeling of the surface static displacements and fault plane slip for the 1979 Imperial Valley earthquake

Three-dimensional finite element modeling techniques are used to synthesize geodetic and seismological results for 1979 Imperial Valley earthquake. The strategy pursued consists of two principal steps. In the first step, the seismologically-derived coseismic fault slip is taken as a function of position in the fault plane and is applied directly to the three-dimensional dislocation model. In the second step, a physical model of stresses and constitutive parameters is perturbed so as to reproduce the observed fault slip. Hence, the principal features of the coseismic slip pattern are explained by a stress-driven fault model in which: (1) a spatially unresolved asperity is found equivalent to a stress drop of 18 MPa averaged over an area of 15 sq km, and (2) driving stress is essentially absent on the fault segment overlapping the 1940 earthquake rupture zone.

Slade, M. A.↗

A TRMM Microwave Radiometer Rain Rate Estimation Method with Convective and Stratiform Discrimination

Tropical Rainfall Measuring Mission (TRMM) Microwave Imager (TMI) radiometer brightness temperature data in the 85 GHz channel (T85) reveal distinct local minima (T85min) in a regional map containing a Mesoscale Convective System (MCS). This is because of relatively small footprint size (approximately 5.5 km) and strong extinction properties in this channel of the TMI. A map of surface rain rate for that region, deduced from simultaneous measurements made by the Precipitation Radar (PR) on board the TRMM satellite, reveals that these T85(sub min), produced by scattering, correspond to local PR rain maxima. Utilizing the PR rain rate map as a guide, we infer empirically from TMI data the presence of three different kinds of thunderstorms or Cbs. These Cbs are classified as young, mature, and decaying types, and are assumed to have a scale of about 20 km on the average. Two parameters are used to classify these three kinds of Cbs based on the T85 data: a) the magnitude of scattering depression deduced from local T85(sub min) and b) the mean horizontal gradient of T85 around such minima. Knowing the category of a given Cb, we can estimate the rain rate associated with it. Such estimation is done with the help of relationships linking T85min to rain rate in each Cb type. Similarly, a weak background rain rate in all the areas where T85 is less than 260 K is deduced with another relationship linking T85 to rain rate. In our rain retrieval model, this background rain constitutes stratiform rain where the Cbs are absent. Initially, these relationships are optimized or tuned utilizing the PR and TMI data of a few MCS events. After such tuning, the model is applied to independent MCS cases. The areal distribution of light (1-10 mm/hr), moderate (10-20 mm/hr), and intense (>= 20 mm/hr) rain rates are retrieved satisfactorily. Accuracy in the estimates of the light, moderate, and intense rain areas and the mean rain rates associated with such areas in these independent MCS cases is on the average about 15 %. Taking advantage of this ability of our retrieval method, one could derive the latent heat input into the atmosphere over the 760 km wide swath of the TMI radiometer in the tropics.

Prabhakara, C.↗

TMI Rain Rate Estimation Over Land and Ocean Utilizing Convective and Stratiform Discrimination

Tropical Rainfall Measuring Mission (TRMM) Microwave Imager (TMI) radiometer brightness temperature data in the 85 GHz channel (T85) reveal distinct local minima in a regional map containing a Mesoscale Convective System (MCS). This is because of relatively small footprint size (approximately 5.5 km) and strong extinction properties in this channel of the TMI. A map of rain rate for that region, deduced from simultaneous measurements made by the Precipitation Radar (PR) on board the TRMM satellite, reveals that these T85 minima, produced by scattering, correspond to local PR rain maxima. Utilizing the PR rain rate map as a guide, we infer from TMI data the presence of three different kinds of thunderstorms or Cbs. They are young, mature, and decaying Cbs that have a scale of about 20 km on the average. Two parameters enable us to infer these three kinds of Cbs objectively: a) the magnitude of scattering depression deduced from local T85 minima and b) the mean horizontal gradient of T85 around such minima. Knowing the category of a given Cb, we can estimate the rain rate associated with it. Such estimation is done with the help of relationships linking T85 minimum to rain rate in each Cb type. Similarly, a weak background rain rate in all the areas where T85 is less than 260 K is deduced with another relationship linking T85 to rain rate. In our rain retrieval model, this background rain constitutes the stratiform rain where the Cbs are absent. Initially, these relationships are optimized or tuned utilizing the PR and TMI data of a few MCS events. After such tuning, the model is applied to independent MCS cases. The areal distribution of light (1-10 mm/hr), moderate (10-20 mm/hr), and intense (> 20 mm/hr) rain rates are retrieved satisfactorally. Accuracy in the estimates of the light, moderate and intense rain areas and the mean rain rates associated with such areas in these independent MCS cases is on the average about 15%. Taking advantage of this ability of our retrieval method, one could derive the latent heat input into the atmosphere over the 760 km wide swath of the TMI radiometer in the tropics.

Prabhakara, C.↗

Arctic Sea ice by passive microwave observations from the Nimbus-5 Satellite

The results of a dynamic/thermodynamic numerical model of Arctic sea ice are compared with satellite images from the Nimbus 5 electrically scanning microwave radiometer. The model combines aspects of two previous sea ice models those of Parkinson and Washington and Ling, Rasmussen, and Campbell. A solid/fluid model basically follows the formulation of the Parkinson and Washington model with the addition of the constitutive equation and equation of state from the Ling model. The Parkinson and Washington model simulates the seasonal cycle of sea ice thicknesses and concentrations with a horizontal resolution of roughly 200 km and a timestep of 8 hours. The thermodynamics are calculated through energy balances at the interfaces between ice and air, water and ice, and water and air. The ice dynamics are calculated through a momentum equation balancing air stress, water stress, dynamic topography, and Coriolis force, with an adjustment for internal ice resistance.

Campbell, W. J.↗

Prediction of stress-strain response of SCS-6/Timetal-21S subjected to a hypersonic flight profile

Thermomechanical response of a cross-ply SCS-6/Timetal-21S composite subjected to a generic hypersonic flight profile with the temperature ranging from -130 C to 816 C was evaluated experimentally and analytically. A two dimensional micromechanical anlaysis, VISCOPLY, was used to predict the stress-strain response of the laminate and of the constituents in each ply during thermomechanical loading conditions. In the analysis, the fiber was modeled as elastic with transverse orthotropic and temperature dependent properties and the matrix was modeled using a thermoviscoplastic constitutive relation. The fiber transverse modulus was reduced in the analysis to simulate fiber-matrix interface failure. Reasonable agreement was found between measured and predicted laminate stress-strain response when fiber-matrix debonding was modeled.

Mirdamadi, Massoud↗

Prediction of thermal and mechanical stress-strain responses of TMC's subjected to complex TMF histories

This paper presents an experimental and analytical evaluation of cross-plied laminates of Ti-15V-3Cr-3Al-3Sn (Ti-15-3) matrix reinforced with continuous silicon-carbide fibers (SCS-6) subjected to a complex TMF loading profile. Thermomechanical fatigue test techniques were developed to conduct a simulation of a generic hypersonic flight profile. A micromechanical analysis was used. The analysis predicts the stress-strain response of the laminate and of the constituents in each ply during thermal and mechanical cycling by using only constituent properties as input. The fiber was modeled as elastic with transverse orthotropic and temperature-dependent properties. The matrix was modeled using a thermoviscoplastic constitutive relation. The fiber transverse modulus was reduced in the analysis to simulate the fiber-matrix interface failures. Excellent correlation was found between measured and predicted laminate stress-strain response due to generic hypersonic flight profile when fiber debonding was modeled.

Johnson, W. S.↗

The Analysis of Adhesively Bonded Advanced Composite Joints Using Joint Finite Elements

The design and sizing of adhesively bonded joints has always been a major bottleneck in the design of composite vehicles. Dense finite element (FE) meshes are required to capture the full behavior of a joint numerically, but these dense meshes are impractical in vehicle-scale models where a course mesh is more desirable to make quick assessments and comparisons of different joint geometries. Analytical models are often helpful in sizing, but difficulties arise in coupling these models with full-vehicle FE models. Therefore, a joint FE was created which can be used within structural FE models to make quick assessments of bonded composite joints. The shape functions of the joint FE were found by solving the governing equations for a structural model for a joint. By analytically determining the shape functions of the joint FE, the complex joint behavior can be captured with very few elements. This joint FE was modified and used to consider adhesives with functionally graded material properties to reduce the peel stress concentrations located near adherend discontinuities. Several practical concerns impede the actual use of such adhesives. These include increased manufacturing complications, alterations to the grading due to adhesive flow during manufacturing, and whether changing the loading conditions significantly impact the effectiveness of the grading. An analytical study is conducted to address these three concerns. Furthermore, proof-of-concept testing is conducted to show the potential advantages of functionally graded adhesives. In this study, grading is achieved by strategically placing glass beads within the adhesive layer at different densities along the joint. Furthermore, the capability to model non-linear adhesive constitutive behavior with large rotations was developed, and progressive failure of the adhesive was modeled by re-meshing the joint as the adhesive fails. Results predicted using the joint FE was compared with experimental results for various joint configurations, including double cantilever beam and single lap joints.

Stapleton, Scott E.↗

Analysis of Sample Acquisition Dynamics Using Discrete Element Method

The analysis presented in this paper is conducted in the framework of the Ocean Worlds Autonomy Testbed for Exploration Research and Simulation (OceanWATERS) project, currently under development at NASA Ames Research Center. OceanWATERS aims at designing a simulation environment which allows for testing autonomy of scientific lander missions to the icy moons of our solar system. Mainly focused on reproducing the end effector interaction with the inherent terrain, this paper introduces a novel discrete element method (DEM)-based approach to determine forces and torques acting on the lander’s scoop during the sample acquisition process. An accurate force feedback from the terrain on the scoop is required by fault-detection and autonomous decision-making algorithms to identify when the requested torque on the robotic arm’s joints exceeds the maximum available torque. Knowledge of the terrain force feedback significantly helps evaluating the arm’s links structural properties and properly selecting actuators for the joints. Models available in literature constitute a partial representation of the dynamics of the interaction. As an example, Balovnev derived an analytical expression of the vertical and horizontal force acting on a bucket while collecting a sample as a function of its geometry and velocity, soil parameters and reached depth. Although the model represents an adequate approximation of the two force components, it ignores the direction orthogonal to the scoop motion and neglects the torque. This work relies on DEM analysis to compensate for analytical models’ deficiencies and inaccuracies, i. e. provide force and torque 3D vectors, defined in the moving reference (body) frame attached to the scoop, at each instant of the sample collection process. Results from the first presented analysis relate to the specific OceanWATERS sampling strategy, which consists of collecting the sample through five consecutive passes with increasing depth, each pass following the same circularlinear- circular trajectory. Data is collected given a specific scoop design interacting with two types of bulk materials, which may characterize the surface of icy planetary bodies: snow and ice. Although specifically concerned with the OceanWATERS design, this first analysis provides the expected force trends for similar sampling strategies and allows to deduce phenomenological information about the general scooping process. In order to further instruct the community on the use of DEM tools as a solution to the sampling collection problem, two more analyses have been carried out, mainly focused on reducing the DEM computation time, which increases with a decrease in particle size. After running a set of identical simulations, where the only changing parameter is the size of the spherical particle, it is observed that the resulting force trajectories, starting from a given particle size, converge to the true trend. It is deducible that a further decrease in size yields negligible improvements in the accuracy, while it sensibly increases computation time. A final analysis aims at discussing limitations of approximating bulk material particles having a complex shape, e. g. ice fragments, with spheres, by comparing force trends resulting in the two cases for the same simulation scenario.

Catanoso, Damiana↗

Notched fatigue of single crystal PWA 1480 at turbine attachment temperatures

The focus is on the lower temperature, uncoated and notched features of gas turbine blades. Constitutive and fatigue life prediction models applicable to these regions are being developed. Fatigue results are presented which were obtained thus far. Fatigue tests are being conducted on PWA 1480 single crystal material using smooth strain controlled specimens and three different notched specimens. Isothermal fatigue tests were conducted at 1200, 1400, and 1600 F. The bulk of the tests were conducted at 1200 F. The strain controlled tests were conducted at 0.4 percent per second strain rate and the notched tests were cycled at 1.0 cycle per second. A clear orientation dependence is observed in the smooth strain controlled fatigue results. The fatigue lifes of the thin, mild notched specimens agree fairly well with this smooth data when elastic stress range is used as a correlating parameter. Finite element analyses were used to calculate notch stresses. Fatigue testing will continue to further explore the trends observed thus far. Constitutive and life prediction models are being developed.

Meyer, T. G.↗

Multiobjective optimization of hybrid regenerative life support technologies. Topic D: Technology Assessment

We developed simple mathematical models for many of the technologies constituting the water reclamation system in a space station. These models were employed for subsystem optimization and for the evaluation of the performance of individual water reclamation technologies, by quantifying their operational 'cost' as a linear function of weight, volume, and power consumption. Then we performed preliminary investigations on the performance improvements attainable by simple hybrid systems involving parallel combinations of technologies. We are developing a software tool for synthesizing a hybrid water recovery system (WRS) for long term space missions. As conceptual framework, we are employing the state space approach. Given a number of available technologies and the mission specifications, the state space approach would help design flowsheets featuring optimal process configurations, including those that feature stream connections in parallel, series, or recycles. We visualize this software tool to function as follows: given the mission duration, the crew size, water quality specifications, and the cost coefficients, the software will synthesize a water recovery system for the space station. It should require minimal user intervention. The following tasks need to be solved for achieving this goal: (1) formulate a problem statement that will be used to evaluate the advantages of a hybrid WRS over a single technology WBS; (2) model several WRS technologies that can be employed in the space station; (3) propose a recycling network design methodology (since the WRS synthesis task is a recycling network design problem, it is essential to employ a systematic method in synthesizing this network); (4) develop a software implementation for this design methodology, design a hybrid system using this software, and compare the resulting WRS with a base-case WRS; and (5) create a user-friendly interface for this software tool.

Manousiouthakis, Vasilios↗

Computational Analysis of the X-57 Maxwell Airplane, the Landing Configuration with High-Lift Blowing and Aileron Deflections (Preliminary Fuselage)

The X-57 Maxwell is an all-electric airplane with a distributed electric propulsion system used for a high-lift system at takeoff and landing conditions. The Kestrel and USM3D flow solvers were used at NASA Langley to investigate the performance of the X-57 Maxwell in the development of an aerodynamic database. The configuration investigated in this paper had a 30° flap deflection, and in addition the pilot’s right aileron were deflected and undeflected in different cases. The solutions were computed at an airspeed of 58 KEAS, for an altitude of 6000 feet, and a flight Reynolds number of 0.588e+06 per foot. To evaluate the high-lift distributed electric propulsion system without aileron deflections, the solutions were computed for an angle-of-attack sweep from −2° to 20°, with the high-lift propellers blowing and the cruise propellers excluded from the simulation. To investigate the aileron effectiveness, the solutions were computed at angles of attack of −2° and 14°, for aileron deflections from −25° to 18° with the high-lift propellers blowing and the cruise propellers windmilling at idle-power. The high-lift propellers and the idle-power cruise propellers were modeled with an actuator disk. Results show negligible differences in lift, drag and pitching moment whether the idle-powered cruise propellers were included or excluded from the simulation. In general, the Kestrel and USM3D codes compared well for lift, drag and pitching moment for the landing configuration with no aileron control. The Kestrel code, using the Spalart-Allmaras turbulence model with the rotation corrections terms, predicted an increasing lift coefficient up to maximum lift coefficient of 4.65 at a 15° angle of attack. The USM3D code, using the Spalart-Allmaras turbulence model turbulence model with the Quadratic Constitutive Relation predicted an increasing lift coefficient up to maximum lift coefficient of 4.5 at a 12° angle of attack, with the lift remaining constant through a 14◦ angle of attack. A possible difference in lift coefficient between the codes for high angles of attack may result from the different available options used with the standard Spalart-Allmaras turbulence model. The Kestrel code predicted more drag across the range of angle of attack than the USM3D code. The codes compared well for pitching moment coefficient across the range of angle of attack. The Kestrel and USM3D codes compared well for aileron effectiveness at a 2° angle of attack. The Kestrel code predicts better aileron effectiveness than USM3D at a 14° angle of attack. There is more flow separation in the region outboard of the last high-lift nacelle for the USM3D solutions at a 14° angle of attack, which diminishes the ability of the aileron to be effective.

X-57 Maxwell Airplane↗

A complete model of turbulence

A set of constitutive equations suitable for a priori computation of turbulent shear flows has been developed. Since no properties of a given turbulent flow need be known in advance in order to obtain a solution, the equations comprise a complete model of turbulence. Perturbation analysis shows that the model predicts a composite five-layer structure for an incompressible turbulent boundary layer, viz, a defect layer, a law-of-the-wall layer, a viscous sublayer, a near-surface roughness layer, and a viscous superlayer at the boundary-layer edge. Analysis of the defect layer demonstrates the key improvement of the model over its predecessor, the Saffman-Wilcox two-equation model of turbulence. Examination of model-predicted sublayer structure yields model-parameter boundary conditions appropriate for surfaces with roughness and mass injection. Results of numerical computations of compressible and incompressible equilibrium boundary layers show that, for such flows, the model is as accurate as mixing-length theory. Applications to transitional boundary layers and to nonequilibrium relaxation of a boundary layer passing from a rough to a smooth surface indicate that the model's applicability extends far beyond that of mixing-length theory's.

Wilcox, D. C.↗

Modeling Progressive Failure of Bonded Joints Using a Single Joint Finite Element

Enhanced finite elements are elements with an embedded analytical solution which can capture detailed local fields, enabling more efficient, mesh-independent finite element analysis. In the present study, an enhanced finite element is applied to generate a general framework capable of modeling an array of joint types. The joint field equations are derived using the principle of minimum potential energy, and the resulting solutions for the displacement fields are used to generate shape functions and a stiffness matrix for a single joint finite element. This single finite element thus captures the detailed stress and strain fields within the bonded joint, but it can function within a broader structural finite element model. The costs associated with a fine mesh of the joint can thus be avoided while still obtaining a detailed solution for the joint. Additionally, the capability to model non-linear adhesive constitutive behavior has been included within the method, and progressive failure of the adhesive can be modeled by using a strain-based failure criteria and re-sizing the joint as the adhesive fails. Results of the model compare favorably with experimental and finite element results.

Stapleton, Scott E.↗