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 559 records · Page 31

Third and Final Shuttle Mission of the Isothermal Dendritic Growth Experiment Conducted: Highest Supercooling Ever Recorded Achieved

Dendrites describe the treelike crystal morphology commonly assumed in metals and alloys that freeze from supercooled or supersaturated melts. There remains a high level of engineering interest in dendritic solidification because the size, shape, and orientation of the dendrites determine the final microstructure of a material. It is the microstructure that then determines the physical properties of cast or welded products. Although it is well known that dendritic growth is controlled by the transport of latent heat from the moving solid-liquid interface, an accurate and predictive model has not yet been developed. The effects of gravity-induced convection on the transfer of heat from the interface have prevented adequate testing, under terrestrial conditions, of solidification models. The Isothermal Dendritic Growth Experiment (IDGE) constituted a series of three microgravity experiments flown aboard the Space Shuttle Columbia. The apparatus was used to grow and record dendrite solidification in the absence of gravity-induced convective heat transfer, thereby producing a wealth of benchmark-quality data for testing solidification models and theories.

Glicksman, Martin E.↗

Ab initio ground states of strongly-correlated materials on quantum computers

The accurate first-principles description of strongly-correlated materials is an important and challenging problem in condensed matter physics. Ab initio downfolding has emerged as a way of deriving accurate many-body Hamiltonians including strong correlations, representing a subspace of interest of a material, using density functional theory calculations as a starting point. However, the solution of these material-specific models can scale exponentially on classical computers, constituting a challenge. Here we propose that utilizing quantum computers for obtaining the properties of downfolded Hamiltonians yields an accurate description of the ground state properties of strongly-correlated systems, while circumventing the exponential scaling problem. We benchmark the solution of Hubbard-like models obtained through downfolding by utilizing a classical tensor network implementation of variational quantum eigensolvers (VQE), and we reveal a strategy for driving the optimization through a hybrid minimization of the energy and maximization of the overlap with an approximate solution obtained through low-cost computational methods. This results in a reduction of the energy error by orders of magnitude compared to conventional VQE approaches, and allows us to reproduce long-range correlations for the first time. We demonstrate our first-principles approach for diverse strongly-correlated materials, correctly predicting the antiferromagnetic state of one-dimensional cuprate Ca 2 CuO 3 , the excitonic ground state of monolayer WTe2, and the charge-ordered state of correlated metal SrVO 3 . Our efficient computational implementation allows us to simulate large systems with up to 54 qubits and encompassing up to four correlated bands, which is indicative of the complexity that our framework can address.

Antonios M Alvertis↗

DPW-VI Results Using FUN3D with Focus on k-kL-MEAH2015 (k-kL) Turbulence Model

The Common Research Model wing-body configuration is investigated with the k-kL-MEAH2015 turbulence model implemented in FUN3D. This includes results presented at the Sixth Drag Prediction Workshop and additional results generated after the workshop with a nonlinear Quadratic Constitutive Relation (QCR) variant of the same turbulence model. The workshop provided grids are used, and a uniform grid refinement study is performed at the design condition. A large variation between results with and without a reconstruction limiter is exhibited on "medium" grid sizes, indicating that the medium grid size is too coarse for drawing conclusions in comparison with experiment. This variation is reduced with grid refinement. At a fixed angle of attack near design conditions, the QCR variant yielded decreased lift and drag compared with the linear eddy-viscosity model by an amount that was approximately constant with grid refinement. The k-kL-MEAH2015 turbulence model produced wing root junction flow behavior consistent with wind tunnel observations.

Abdol-Hamid, K. S.↗

Analytical Modeling of the High Strain Rate Deformation of Polymer Matrix Composites

The results presented here are part of an ongoing research program to develop strain rate dependent deformation and failure models for the analysis of polymer matrix composites subject to high strain rate impact loads. State variable constitutive equations originally developed for metals have been modified in order to model the nonlinear, strain rate dependent deformation of polymeric matrix materials. To account for the effects of hydrostatic stresses, which are significant in polymers, the classical 5 plasticity theory definitions of effective stress and effective plastic strain are modified by applying variations of the Drucker-Prager yield criterion. To verify the revised formulation, the shear and tensile deformation of a representative toughened epoxy is analyzed across a wide range of strain rates (from quasi-static to high strain rates) and the results are compared to experimentally obtained values. For the analyzed polymers, both the tensile and shear stress-strain curves computed using the analytical model correlate well with values obtained through experimental tests. The polymer constitutive equations are implemented within a strength of materials based micromechanics method to predict the nonlinear, strain rate dependent deformation of polymer matrix composites. In the micromechanics, the unit cell is divided up into a number of independently analyzed slices, and laminate theory is then applied to obtain the effective deformation of the unit cell. The composite mechanics are verified by analyzing the deformation of a representative polymer matrix composite (composed using the representative polymer analyzed for the correlation of the polymer constitutive equations) for several fiber orientation angles across a variety of strain rates. The computed values compare favorably to experimentally obtained results.

Goldberg, Robert K.↗

Equations of state and constitutive equations

Equations-of state information in the otherwise undeformed state as a starting point for the development of constitutive equations is considered in addition to thermodynamics, free energy functions, and conceptual difficulties including the definition of reference states for strain. V-T effects in the form of the Simha-Somcynsky (1969) equation of state are explicitly discussed, and it is shown how this model can be modified to produce a constitutive equation. Continuum mechanics approaches are considered, and examples are given of developments based on linear viscoelastic theory which directly incorporate stress-induced volume changes, and on large-strain elastic theory.

Landel, Robert F.↗

Advances in Aero-Propulsive Modeling for Fixed-Wing and eVTOL Aircraft Using Experimental Data

Small unmanned aircraft and electric vertical takeoff and landing (eVTOL) aircraft have recently emerged as vehicles able to perform new missions and stimulate future air transportation methods. This dissertation presents several system identification research advancements for these modern aircraft configurations enabling accurate mathematical model development for flight dynamics simulations based on wind-tunnel and flight-test data. The first part of the dissertation focuses on advances in flight-test system identification methods using small, fixed-wing, remotely-piloted, electric, propeller-driven aircraft. A generalized approach for flight dynamics model development for small fixed-wing aircraft from flight data is described and is followed by presentation of novel flight-test system identification applications, including: aero-propulsive model development for propeller aircraft and nonlinear dynamic model identification without mass properties. The second part of the dissertation builds on established fixed-wing and rotary-wing aircraft system identification methods to develop modeling strategies for transitioning, distributed propulsion, eVTOL aircraft. Novel wind-tunnel experiment designs and aero-propulsive modeling approaches are developed using a subscale, tandem tilt-wing, eVTOL aircraft, leveraging design of experiments and response surface methodology techniques. Additionally, a method applying orthogonal phase-optimized multisine input excitations to aircraft control effectors in wind-tunnel testing is developed to improve test efficiency and identified model utility. Finally, the culmination of this dissertation is synthesis of the techniques described throughout the document to form a flight-test system identification approach for eVTOL aircraft that is demonstrated using a high-fidelity flight dynamics simulation. The research findings highlighted throughout the dissertation constitute substantial progress in efficient empirical aircraft modeling strategies that are applicable to many current and future aeronautical vehicles enabling accurate flight simulation development, which can subsequently be used to foster advancement in many other pertinent technology areas.

system identification↗

A structural model of the adaptive human pilot

A compensatory tracking model of the human pilot is offered which attempts to provide a more realistic representation of the human's signal processing structure than that which is exhibited by pilot models currently in use. Two features of the model distinguish it from other representations of the human pilot. First, proprioceptive information from the control stick or manipulator constitutes one of the major feedback paths in the model, providing feedback of vehicle output rate due to control activity. Implicit in this feedback loop is a model of the vehicle dynamics which is valid in and beyond the region of crossover. Second, error-rate information is continuously derived and independently but intermittently controlled. An output injected remnant model is offered and qualitatively justified on the basis of providing a measure of the effect of inaccuracies such as time variations in the pilot's internal model of the controlled element dynamics. The data from experimental tracking tasks involving five different controlled element dynamics and one nonideal viewing condition were matched with model generated describing functions and remnant power spectral densities.

Hess, R. A.↗

Thermal Effects Modeling Developed for Smart Structures

Applying smart materials in aeropropulsion systems may improve the performance of aircraft engines through a variety of vibration, noise, and shape-control applications. To facilitate the experimental characterization of these smart structures, researchers have been focusing on developing analytical models to account for the coupled mechanical, electrical, and thermal response of these materials. One focus of current research efforts has been directed toward incorporating a comprehensive thermal analysis modeling capability. Typically, temperature affects the behavior of smart materials by three distinct mechanisms: Induction of thermal strains because of coefficient of thermal expansion mismatch 1. Pyroelectric effects on the piezoelectric elements; 2. Temperature-dependent changes in material properties; and 3. Previous analytical models only investigated the first two thermal effects mechanisms. However, since the material properties of piezoelectric materials generally vary greatly with temperature (see the graph), incorporating temperature-dependent material properties will significantly affect the structural deflections, sensory voltages, and stresses. Thus, the current analytical model captures thermal effects arising from all three mechanisms through thermopiezoelectric constitutive equations. These constitutive equations were incorporated into a layerwise laminate theory with the inherent capability to model both the active and sensory response of smart structures in thermal environments. Corresponding finite element equations were formulated and implemented for both the beam and plate elements to provide a comprehensive thermal effects modeling capability.

Lee, Ho-Jun↗

Preliminary Work for Modeling the Propellers of an Aircraft as a Noise Source in an Acoustic Boundary Element Analysis

An algorithm for generating appropriate velocity boundary conditions for an acoustic boundary element analysis from the kinematics of an operating propeller is presented. It constitutes the initial phase of Integrating sophisticated rotorcraft models into a conventional boundary element analysis. Currently, the pressure field is computed by a linear approximation. An initial validation of the developed process was performed by comparing numerical results to test data for the external acoustic pressure on the surface of a tilt-rotor aircraft for one flight condition.

Vlahopoulos, Nickolas↗

Analysis of Compressible Mixing Layers Using Dilatational Covariances Model

Compressible mixing layers are analyzed using a dilatational covariances model based on a pseudo-sound constitutive relation. The calculations are used to evaluate the different physical phenomena affecting compressible mixing layers. The rate of growth of the mixing layer is retarded by both the compressible dissipation and the pressure-dilatational covariances. The pressure-dilatational, essentially a nonequilibrium effect, reduces the amount of excess production over dissipation available for the turbulence energy growth. The pseudo-sound model also includes a history dependent portion: this is also investigated. All constants in the model and used in these computations are predicted by the theory.

Thangam, S.↗

Nonlinearity and Strain-Rate Dependence in the Deformation Response of Polymer Matrix Composites Modeled

There has been no accurate procedure for modeling the high-speed impact of composite materials, but such an analytical capability will be required in designing reliable lightweight engine-containment systems. The majority of the models in use assume a linear elastic material response that does not vary with strain rate. However, for containment systems, polymer matrix composites incorporating ductile polymers are likely to be used. For such a material, the deformation response is likely to be nonlinear and to vary with strain rate. An analytical model has been developed at the NASA Glenn Research Center at Lewis Field that incorporates both of these features. A set of constitutive equations that was originally developed to analyze the viscoplastic deformation of metals (Ramaswamy-Stouffer equations) was modified to simulate the nonlinear, rate-dependent deformation of polymers. Specifically, the effects of hydrostatic stresses on the inelastic response, which can be significant in polymers, were accounted for by a modification of the definition of the effective stress. The constitutive equations were then incorporated into a composite micromechanics model based on the mechanics of materials theory. This theory predicts the deformation response of a composite material from the properties and behavior of the individual constituents. In this manner, the nonlinear, rate-dependent deformation response of a polymer matrix composite can be predicted.

Goldberg, Robert K.↗

Does The Earth Have an Adaptive Infrared Iris?

Observations and analyses of water vapor and clouds in the tropics over the past decade suggest a different approach to radiative climate feedbacks: namely, that high clouds and high free-tropospheric relative humidity are largely tied to each other, and that the main feedback consists in changing the relative areas of cloudy/moist regions vis a vis clear/dry regions in response to the surface temperature of the cloudy/moist regions - as opposed to altering the humidity in either of the regions. This is an intrinsically 2-dimensional (horizontal and vertical) effect which does not readily enter simple 1-dimensional (vertical) radiative-convective schemes which emphasize average humidity, etc. Preliminary analyses of cloud data for the eastern part of the Western Pacific from the Japanese GMS-5(Geostationary Meteorological Satellite), are supportive of this suggestion - pointing to a 15% reduction in cloudy/moist area for a 1C increase of the sea surface temperature as measured by the cloud-weighted SST (sea surface temperature). The implication of this result is examined using a simple 2-dimensional radiative-convective model. The calculations show that such a change in the tropics would lead to a strong negative feedback in the global climate, with a feedback factor of about -1.7, which, if correct, would easily dominate the positive water vapor feedback found in current models. This new feedback mechanism, in effect, constitutes an adaptive infrared iris that opens and closes in order to control the OLR (outgoing longwave radiation) in response to changes in surface temperature in a manner similar to the way in which an eye's iris opens and closes in response to changing light levels. The climate sensitivity resulting from this thermostatic mechanism is consistent with the independent determination by Lindzen and Giannitisis (1998). Preliminary attempts to replicate observations with GCMs (General Circulation Models) suggest that models lack such a negative cloud/moist areal feedback.

Lindzen, Richard S.↗

Mathmatical modeling for diffractive optics

We consider a 'diffractive optic' to be a biperiodic surface separating two half-spaces, each having constant constitutive parameters; within a unit cell of the periodic surface and across the transition zone between the two half-spaces, the constitutive parameters can be a continuous, complex-valued function. Mathematical models for diffractive optics have been developed, and implemented as numerical codes, both for the 'direct' problem and for the 'inverse' problem. In problems of the 'direct' class, the diffractive optic is specified, and the full set of Maxwell's equations is cast in a variational form and solved numerically by a finite element approach. This approach is well-posed in the sense that existence and uniqueness of the solution can be proved and specific convergence conditions can be derived. An example of a metallic grating at a Wood anomaly is presented as a case where other approaches are known to have convergence problems. In problems of the 'inverse' class, some information about the diffracted field (e.g., the far-field intensity) is given, and the problem is to find the periodic structure in some optimal sense. Two approaches are described: phase reconstruction in the far-field approximation; and relaxed optimal design based on the Helmholtz equation. Practical examples are discussed for each approach to the inverse problem, including array generators in the far-field case and antireflective structures for the relaxed optimal design.

Dobson, David↗

Contributions to HiLiftPW-3 Using Structured, Overset Grid Methods

The High-Lift Common Research Model (HL-CRM) and the JAXA Standard Model (JSM) were analyzed computationally using both the OVERFLOW and LAVA codes for the third AIAA High-Lift Prediction Workshop. Geometry descriptions and the test cases simulated are described. With the HL-CRM, the effects of surface smoothness during grid projection and the effect of partially sealing a flap gap were studied. Grid refinement studies were performed at two angles of attack using both codes. For the JSM, simulations were performed with and without the nacelle/pylon. Without the nacelle/pylon, evidence of multiple solutions was observed when a quadratic constitutive relation is used in the turbulence modeling; however, using time-accurate simulation seemed to alleviate this issue. With the nacelle/pylon, no evidence of multiple solutions was observed. Laminar-turbulent transition modeling was applied to both JSM configuration, and had an overall favorable impact on the lift predictions.

HiLiftPW-3↗

User-Defined Material Model for Progressive Failure Analysis

An overview of different types of composite material system architectures and a brief review of progressive failure material modeling methods used for structural analysis including failure initiation and material degradation are presented. Different failure initiation criteria and material degradation models are described that define progressive failure formulations. These progressive failure formulations are implemented in a user-defined material model (or UMAT) for use with the ABAQUS/Standard1 nonlinear finite element analysis tool. The failure initiation criteria include the maximum stress criteria, maximum strain criteria, the Tsai-Wu failure polynomial, and the Hashin criteria. The material degradation model is based on the ply-discounting approach where the local material constitutive coefficients are degraded. Applications and extensions of the progressive failure analysis material model address two-dimensional plate and shell finite elements and three-dimensional solid finite elements. Implementation details and use of the UMAT subroutine are described in the present paper. Parametric studies for composite structures are discussed to illustrate the features of the progressive failure modeling methods that have been implemented.

Knight, Norman F. Jr.↗

Factors Influencing Progressive Failure Analysis Predictions for Laminated Composite Structure

Progressive failure material modeling methods used for structural analysis including failure initiation and material degradation are presented. Different failure initiation criteria and material degradation models are described that define progressive failure formulations. These progressive failure formulations are implemented in a user-defined material model for use with a nonlinear finite element analysis tool. The failure initiation criteria include the maximum stress criteria, maximum strain criteria, the Tsai-Wu failure polynomial, and the Hashin criteria. The material degradation model is based on the ply-discounting approach where the local material constitutive coefficients are degraded. Applications and extensions of the progressive failure analysis material model address two-dimensional plate and shell finite elements and three-dimensional solid finite elements. Implementation details are described in the present paper. Parametric studies for laminated composite structures are discussed to illustrate the features of the progressive failure modeling methods that have been implemented and to demonstrate their influence on progressive failure analysis predictions.

Knight, Norman F., Jr.↗

Two-Fluid Models and Interfacial Area Transport in Microgravity Condition

The objective of the present study is to develop a two-fluid model formulation with interfacial area transport equation applicable for microgravity conditions. The new model is expected to make a leapfrog improvement by furnishing the constitutive relations for the interfacial interaction terms with the interfacial area transport equation, which can dynamically model the changes of the interfacial structures. In the first year of this three-year project supported by the U.S. NASA, Office of Biological and Physics Research, the primary focus is to design and construct a ground-based, microgravity two-phase flow simulation facility, in which two immiscible fluids with close density will be used. In predicting the two-phase flow behaviors in any two-phase flow system, the interfacial transfer terms are among the most essential factors in the modeling. These interfacial transfer terms in a two-fluid model specify the rate of phase change, momentum exchange, and energy transfer at the interface between the two phases. For the two-phase flow under the microgravity condition, the stability of the fluid particle interface and the interfacial structures are quite different from those under normal gravity condition. The flow structure may not reach an equilibrium condition and the two fluids may be loosely coupled such that the inertia terms of each fluid should be considered separately by use of the two-fluid model. Previous studies indicated that, unless phase-interaction terms are accurately modeled in the two-fluid model, the complex modeling does not necessarily warrant an accurate solution.

Ishii, Mamoru↗