Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stress Modeling”

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 505 records · Page 28

A kinematically driven anisotropic viscoelastic constitutive model applied to tires

Aircraft tires are composite structures manufactured with viscoelastic materials such as carbon black filled rubber and nylon cords. When loaded they experience large deflections and moderately large strains. Detailed structural models of tires require the use of either nonlinear shell or nonlinear three dimensional solid finite elements. Computational predictions of the dynamic response of tires must consider the composite viscoelastic material behavior in a realistic fashion. We describe a modification to a nonlinear anisotropic shell finite element so it can be used to model viscoelastic stresses during general deformations. The model is developed by introducing internal variables of the type used to model elastic strain energy. The internal variables are strains, curvatures, and transverse shear angles which are in a one-to-one correspondence with the generalized coordinates used to model the elastic strain energy for nonlinear response. A difference-relaxation equation is used to relate changes in the observable strain field to changes in the internal strain field. The internal stress state is introduced into the equilibrium equations by converting it to nodal loads associated with the element's displacement degrees of freedom. In this form the tangent matrix in the Newton-Raphson solution algorithm is not modified from its form for the nonlinear statics problem. Only the gradient vector is modified and the modification is not computationally costly. The existing finite element model for the Space Shuttle nose gear tire is used to provide examples of the algorithm. In the first example, the tire's rim is displaced at a constant rate up to a fixed value. In the second example, the tire's rim is enforced to follow a saw tooth load and unload curve to generate hysteresis loops.

Johnson, Arthur R.↗

Computational Models of the Cardiovascular System and Its Response to Microgravity

Computational models of the cardiovascular system are powerful adjuncts to ground-based and in-flight experiments. We will provide NSBRI with a model capable of simulating the short-term effects of gravity on cardiovascular function. The model from this project will: (1) provide a rational framework which quantitatively defines interactions among complex cardiovascular parameters and which supports the critical interpretation of experimental results and testing of hypotheses. (2) permit predictions of the impact of specific countermeasures in the context of various hypothetical cardiovascular abnormalities induced by microgravity. Major progress has been made during the first 18 months of the program: (1) We have developed an operational first-order computer model capable of simulating the cardiovascular response to orthostatic stress. The model consists of a lumped parameter hemodynamic model and a complete reflex control system. The latter includes cardiopulmonary and carotid sinus reflex limbs and interactions between the two. (2) We have modeled the physiologic stress of tilt table experiments and lower body negative pressure procedures (LBNP). We have verified our model's predictions by comparing them with experimental findings from the literature. (3) We have established collaborative efforts with leading investigators interested in experimental studies of orthostatic intolerance, cardiovascular control, and physiologic responses to space flight. (4) We have established a standardized method of transferring data to our laboratory from the ongoing NSBRI bedrest studies. We use this data to estimate input parameters to our model and compare our model predictions to actual data to further verify our model. (5) We are in the process of systematically simulating current hypotheses concerning the mechanism underlying orthostatic intolerance by matching our simulations to stand test data from astronauts pre- and post-flight. (6) We are in the process of developing a JAVA version of the simulator which will be distributed amongst the cardiovascular team members. Future work on this project involves modifications of the model to represent a rodent (rat) model, further evaluation of the bedrest astronaut and animal data, and systematic investigation of specific countermeasures.

Kamm, Roger D.↗

Computation of turbulent near wake for asymmetric airfoils

A numerical procedure for studying the turbulent near wake of two dimensional airfoil sections is presented. The Reynolds Navier-Stokes equations were written for flow about bodies of arbitrary geometry and solved on an arbitrary nonuniform curvilinear computational mesh. Eddy viscosity and Reynolds stress turbulence transport models are considered. Specific examples are shown for airfoil section by using an algebraic viscosity model with streamwise relaxation and the interactive Reynolds stress model.

Deiwert, G. S.↗

Modified 316H constitutive model updated with high temperature stress relaxation test data

This report describes the calibration of a new high temperature constitutive model for 316H stainless steel, suitable for use with the ASME Boiler & Pressure Vessel Section III, Division 5, Class A rules for design by inelastic analysis. The model retains the same mathematical form used by the reference model included in Nonmandatory Appendix Z of the Code, but refits the model to an expanded dataset including all the data used to fit the original model plus seven new stress relaxation tests. The addition of these high temperature stress relaxation tests improves the model's accuracy in predicting relaxation at temperatures greater than 700 ⁰C, without compromising the accuracy of the model versus the original calibration data.

36 MATERIALS SCIENCE↗

NASA Trapezoidal Wing Simulation Using Stress-w and One- and Two-Equation Turbulence Models

The Wilcox 2006 stress-omega model (also referred to as WilcoxRSM-w2006) has been implemented in the NASA Langley code CFL3D and used to study a variety of 2-D and 3-D configurations. It predicted a variety of basic cases reasonably well, including secondary flow in a supersonic rectangular duct. One- and two-equation turbulence models that employ the Boussinesq constitutive relation were unable to predict this secondary flow accurately because it is driven by normal turbulent stress differences. For the NASA trapezoidal wing at high angles of attack, the WilcoxRSM-w2006 model predicted lower maximum lift than experiment, similar to results of a two-equation model.

Rodio, J. J.↗

Modification of the MTS model for high strain-rate behavior of TI-6AL-4V

The Mechanical Threshold Stress (MTS) model provides excellent predictive capabilities for the material constitutive response for a wide range of temperatures and strain rates. However, the MTS model fails to capture the rapidly increasing yield stress at high strain rate behavior as the deformation controlling mechanism transitions from thermal activation to drag mechanisms, only capturing the linear behavior. Further, the model typically over predicts the flow stress behavior at yield and post yield due to its use of a constant work hardening rate parameter derived from the stress–strain response at constant saturation stress. An alternative approach to fitting portions of the MTS model is investigated and mathematical models are developed to address these issues. The results show that with appropriate experimental data, the mechanical threshold stress and work hardening rate parameters within the MTS model can quite easily and accurately be modified to extend applicability to high strain rate behavior and more accurately model the initial flow stress behavior at early work hardening rates without modification of the functions core to the MTS model itself.

36 MATERIALS SCIENCE↗

Turbulence Modeling for Shock Wave/Turbulent Boundary Layer Interactions

Accurate aerodynamic computational predictions are essential for the safety of space vehicles, but these computations are of limited accuracy when large pressure gradients are present in the flow. The goal of the current project is to improve the state of compressible turbulence modeling for high speed flows with shock wave / turbulent boundary layer interactions (SWTBLI). Emphasis will be placed on models that can accurately predict the separated region caused by the SWTBLI. These flows are classified as nonequilibrium boundary layers because of the very large and variable adverse pressure gradients caused by the shock waves. The lag model was designed to model these nonequilibrium flows by incorporating history effects. Standard one- and two-equation models (Spalart Allmaras and SST) and the lag model will be run and compared to a new lag model. This new model, the Reynolds stress tensor lag model (lagRST), will be assessed against multiple wind tunnel tests and correlations. The basis of the lag and lagRST models are to preserve the accuracy of the standard turbulence models in equilibrium turbulence, when the Reynolds stresses are linearly related to the mean strain rates, but create a lag between mean strain rate effects and turbulence when nonequilibrium effects become important, such as in large pressure gradients. The affect this lag has on the results for SWBLI and massively separated flows will be determined. These computations will be done with a modified version of the OVERFLOW code. This code solves the RANS equations on overset grids. It was used for this study for its ability to input very complex geometries into the flow solver, such as the Space Shuttle in the full stack configuration. The model was successfully implemented within two versions of the OVERFLOW code. Results show a substantial improvement over the baseline models for transonic separated flows. The results are mixed for the SWBLI assessed. Separation predictions are not as good as the baseline models, but the over prediction of the peak heat flux downstream of the reattachment shock that plagues many models is reduced.

Lillard, Randolph P.↗

Early warning and crop condition assessment research

The Early Warning Crop Condition Assessment Project of AgRISTARS was a multiagency and multidisciplinary effort. Its mission and objectives were centered around development and testing of remote-sensing techniques that enhance operational methodologies for global crop-condition assessments. The project developed crop stress indicators models that provide data filter and alert capabilities for monitoring global agricultural conditions. The project developed a technique for using NOAA-n satellite advanced very-high-resolution radiometer (AVHRR) data for operational crop-condition assessments. This technology was transferred to the Foreign Agricultural Service of the USDA. The project developed a U.S. Great Plains data base that contains various meteorological parameters and vegetative index numbers (VIN) derived from AVHRR satellite data. It developed cloud screening techniques and scan angle correction models for AVHRR data. It also developed technology for using remotely acquired thermal data for crop water stress indicator modeling. The project provided basic technology including spectral characteristics of soils, water, stressed and nonstressed crop and range vegetation, solar zenith angle, and atmospheric and canopy structure effects.

Boatwright, G. O.↗

The Effect of Creep on the Residual Stresses Generated During Silicon Sheet Growth

The modeling of stresses generated during the growth of thin silicon sheets at high speeds is an important part of the EFG technique since the experimental measurement of the stresses is difficult and prohibitive. The residual stresses which arise in such a growth process lead to serious problems which make thin Si ribbons unsuitable for fabrication. The constitutive behavior is unrealistic because at high temperature (close to the melting point) Si exhibits considerable creep which significantly relaxes the residual stresses. The effect of creep on the residual stresses generated during the growth of Si sheets at high speeds was addressed and the basic qualitative effect of creep are reported.

Hutchinson, J. W.↗

Machine and Deep Learning: Artificial Intelligence Application in Biotic and Abiotic Stress Management in Plants

Biotic and abiotic stresses significantly affect plant fitness, resulting in a serious loss in food production. Biotic and abiotic stresses predominantly affect metabolite biosynthesis, gene and protein expression, and genome variations. However, light doses of stress result in the production of positive attributes in crops, like tolerance to stress and biosynthesis of metabolites, called hormesis. Advancement in artificial intelligence (AI) has enabled the development of high-throughput gadgets such as high-resolution imagery sensors and robotic aerial vehicles, i.e., satellites and unmanned aerial vehicles (UAV), to overcome biotic and abiotic stresses. These High throughput (HTP) gadgets produce accurate but big amounts of data. Significant datasets such as transportable array for remotely sensed agriculture and phenotyping reference platform (TERRA-REF) have been developed to forecast abiotic stresses and early detection of biotic stresses. For accurately measuring the model plant stress, tools like Deep Learning (DL) and Machine Learning (ML) have enabled early detection of desirable traits in a large population of breeding material and mitigate plant stresses. In this review, advanced applications of ML and DL in plant biotic and abiotic stress management have been summarized.

59 BASIC BIOLOGICAL SCIENCES↗

Nonlinear Spring Finite Elements for Predicting Mode I-Dominated Delamination Growth in Laminated Structure with Through-Thickness reinforcement

One particular concern of polymer matrix composite laminates is the relatively low resistance to delamination cracking, in particular when the dominant type of failure is mode I opening. One method proposed for alleviating this problem involves the insertion pultruded carbon pins through the laminate thickness. The pins, known as z-pins, are inserted into the prepreg laminate using an ultrasonic hammer prior to the curing process, resulting in a field of pins embedded normal to the laminate plane as illustrated in Figure. 1. Pin diameters range between 0.28-mm to 0.5-mm and standard areal densities range from 0.5% to 4%. The z-pins are provided by the manufacturer, Aztex(Registered TradeMark) , in a low-density foam preform, which acts to stabilize orientation of the pins during the insertion process [1-3]. Typical pin materials include boron and carbon fibers embedded in a polymer matrix. A number of methods have been developed for predicting delamination growth in laminates reinforced with z-pins. During a study on the effect of z-pin reinforcement on mode I delamination resistance, finite element analyses of z-pin reinforced double cantilever beam (DCB) specimens were performed by Cartie and Partridge [4]. The z-pin bridging stresses were modeled by applying equivalent forces at the pin locations. Single z-pin pull-out tests were performed to characterize the traction law of the pins under mode I loading conditions. Analytical solutions for delamination growth in z-pin reinforced DCB specimens were independently derived by Robinson and Das [5] and Ratcliffe and O'Brien [6]. In the former case, pin bridging stresses were modeled using a distributed load and in the latter example the bridging stresses were discretely modeled by way of grounded springs. Additionally, Robinson and Das developed a data reduction strategy for calculating mode I fracture toughness, G(sub Ic), from a z-pin reinforced DCB specimen test [5]. In both cases a traction law similar to that adopted by Cartie and Partridge was used to represent z-pin failure under mode I loading conditions. In the current work spring elements available in most commercial finite element codes were used to model z-pins. The traction law used in previous analyses [4-6] was employed to represent z-pin damage. This method is intended for and is limited to simulating z-pins in composite laminate structure containing mode I-dominated delamination cracking. The current technique differs from previous analyses in that spring finite elements (available in commercial codes) are employed for simulating zpins, reducing the complexity of the analysis construction process. Furthermore, the analysis method can be applied to general structure that experiences mode I-dominated delamination cracking, in contrast to existing analytical solutions that are only applicable to coupon DCB specimens.

Ratcliffe, James G.↗

Mixing in Low Reynolds Number Reacting Impinging Jets in Crossflow

Previous efforts to model uranyl fluoride formation in an impinging jet gas reactor underpredicted spatial mixing and overpredicted chemical conversion into particulates. The previous fluid dynamics model was based on the solution of the Reynolds Averaged Navier Stokes equations. After simulating fluid dynamics, aerosol dynamics were superimposed onto CFD-simulated gas reactant species concentrations. The current work explores the influence of complex unsteady flow features on the overall flow physics and chemistry for a low Reynolds number, opposed flow, impinging jet gas reactor where there is a low Reynolds number cross flow. The objective of this study was to assess the impact of model formulation on scalar mixing and transport. Here, transient flow simulations were performed using Scale Resolving Simulations. Large-Eddy Simulations with the dynamic Smagorinsky turbulence model were performed along with simulations which directly resolved the flow. Average and root-mean-square (RMS) velocities and species concentrations were computed along with modeled and resolved turbulence kinetic energy (TKE), modeled turbulence dissipation, and modeled turbulent viscosity. Lagrangian flow tracers were also used to quantify species concentrations along path lines emanating from the jet tips. Transient simulation data were compared to results from RANS simulations using the k-ω shear stress transport (SST) model and Reynolds Stress Model (RSM). Transient simulations showed spatial mixing patterns which were more consistent with experimental data and helped elucidate the process of particle formation observed in experiments.

42 ENGINEERING↗

Size effects in models for mechanically-stressed protein crystals and aggregates

As protein aggregates increase in size, they become easier to disrupt mechanically. Using the scaling properties of models proposed to govern protein aggregation, the effect of thermal vibrations and gravity are investigated as deforming forces. For typical protein assemblies made of 30 A proteins, the assembled diameter must remain less than 100-10,000 times the molecular radius to survive in finite thermal and gravity fields. The analysis predicts the following experimental outcomes: (1) reductions in gravitational strain should favor larger protein aggregates; (2) in comparing the aggregate stability of different proteins, the addition of peptide chains should stabilize against thermal strain, but should not affect gravitational strain; (3) critical aggregate sizes should show significant (exponential) sensitivity to cluster geometry, solution preparation and growth conditions. The analysis is extended to consider qualitative size effects in crystal damage during X-ray exposure.

Noever, David A.↗

Computation of confined coflow jets with three turbulence models

A numerical study of confined jets in a cylindrical duct is carried out to examine the performance of two recently proposed turbulence models: an RNG-based K-epsilon model and a realizable Reynolds stress algebraic equation model. The former is of the same form as the standard K-epsilon model but has different model coefficients. The latter uses an explicit quadratic stress-strain relationship to model the turbulent stresses and is capable of ensuring the positivity of each turbulent normal stress. The flow considered involves recirculation with unfixed separation and reattachment points and severe adverse pressure gradients, thereby providing a valuable test of the predictive capability of the models for complex flows. Calculations are performed with a finite-volume procedure. Numerical credibility of the solutions is ensured by using second-order accurate differencing schemes and sufficiently fine grids. Calculations with the standard K-epsilon model are also made for comparison. Detailed comparisons with experiments show that the realizable Reynolds stress algebraic equation model consistently works better than does the standard K-epsilon model in capturing the essential flow features, while the RNG-based K-epsilon model does not seem to give improvements over the standard K-epsilon model under the flow conditions considered.

Zhu, J.↗

Torque Limits for Fasteners in Composites

The two major classes of laminate joints are bonded and bolted. Often the two classes are combined as bonded-bolted joints. Several characteristics of fiber reinforced composite materials render them more susceptible to joint problems than conventional metals. These characteristics include weakness in in-plane shear, transverse tension/compression, interlaminar shear, and bearing strength relative to the strength and stiffness in the fiber direction. Studies on bolted joints of composite materials have been focused on joining assembly subject to in-plane loads. Modes of failure under these loading conditions are net-tension failure, cleavage tension failure, shear-out failure, bearing failure, etc. Although the studies of torque load can be found in literature, they mainly discussed the effect of the torque load on in-plane strength. Existing methods for calculating torque limit for a mechanical fastener do not consider connecting members. The concern that a composite member could be crushed by a preload inspired the initiation of this study. The purpose is to develop a fundamental knowledge base on how to determine a torque limit when a composite member is taken into account. Two simplified analytical models were used: a stress failure analysis model based on maximum stress criterion, and a strain failure analysis model based on maximum strain criterion.

Yi Zhao↗

Research activities at the Center for Modeling of Turbulence and Transition

The main research activities at the Center for Modeling of Turbulence and Transition (CMOTT) are described. The research objective of CMOTT is to improve and/or develop turbulence and transition models for propulsion systems. The flows of interest in propulsion systems can be both compressible and incompressible, three dimensional, bounded by complex wall geometries, chemically reacting, and involve 'bypass' transition. The most relevant turbulence and transition models for the above flows are one- and two-equation eddy viscosity models, Reynolds stress algebraic- and transport-equation models, pdf models, and multiple-scale models. All these models are classified as one-point closure schemes since only one-point (in time and space) turbulent correlations, such as second moments (Reynolds stresses and turbulent heat fluxes) and third moments, are involved. In computational fluid dynamics, all turbulent quantities are one-point correlations. Therefore, the study of one-point turbulent closure schemes is the focus of our turbulence research. However, other research, such as the renormalization group theory, the direct interaction approximation method, and numerical simulations are also pursued to support the development of turbulence modeling.

Shih, Tsan-Hsing↗

Evaluation of Industry Standard Turbulence Models on an Axisymmetric Supersonic Compression Corner

Reynolds-averaged Navier-Stokes computations of a shock-wave/boundary-layer interaction (SWBLI) created by a Mach 2.85 flow over an axisymmetric 30-degree compression corner were carried out. The objectives were to evaluate four turbulence models commonly used in industry, for SWBLIs, and to evaluate the suitability of this test case for use in further turbulence model benchmarking. The Spalart-Allmaras model, Menter's Baseline and Shear Stress Transport models, and a low-Reynolds number k- model were evaluated. Results indicate that the models do not accurately predict the separation location; with the SST model predicting the separation onset too early and the other models predicting the onset too late. Overall the Spalart-Allmaras model did the best job in matching the experimental data. However there is significant room for improvement, most notably in the prediction of the turbulent shear stress. Density data showed that the simulations did not accurately predict the thermal boundary layer upstream of the SWBLI. The effect of turbulent Prandtl number and wall temperature were studied in an attempt to improve this prediction and understand their effects on the interaction. The data showed that both parameters can significantly affect the separation size and location, but did not improve the agreement with the experiment. This case proved challenging to compute and should provide a good test for future turbulence modeling work.

Turbulent Boundary Layer↗