Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “elastic model”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Sensitivity analysis of numerical modeling input parameters on floating offshore wind turbine loads

Abstract. Floating wind turbines must withstand a unique and challenging set of loads from the wind and ocean environment. To de-risk development, accurate predictions of these loads are necessary. Uncertainty in modeling predictions leads to larger required safety factors, increasing project costs and the levelized cost of energy. Complex aero-hydro-elastic modeling tools use many input parameters to represent the wind, waves, current, aerodynamic loads, hydrodynamic loads, and structural properties. It is helpful to understand which of these parameters ultimately drives a design. In this work, an ultimate and fatigue-proxy load sensitivity analysis was performed with 35 different input parameters, using an elementary effects approach to identify the most influential parameters for a case study involving the National Renewable Energy Laboratory (NREL) 5 MW baseline wind turbine atop the OC4-DeepCwind semisubmersible during normal operation. The importance of each parameter was evaluated using 14 response quantities of interest across three operational wind speed conditions. The study concludes that turbulent wind velocity standard deviation is the parameter with the strongest sensitivity; this value is important not just for turbine loads, but also for the global system response. The system center of mass in the wind direction is found to have the highest impact on the system rotation and tower loads. The current velocity is found to be the most dominating parameter for the system global motion and consequently the mooring loads. All tested wind turbulence parameters in addition to the standard deviation are also found to be influential. Wave characteristics are influential for some fatigue-proxy loading but do not significantly impact the extreme ultimate loads in these operational load cases. The required number of random seeds for stochastic environmental conditions is considered to ensure that the sensitivities are due to the input parameters and not due to the seed. The required number of analysis points in the parameter space is identified so that the conclusions represent a global sensitivity. The results are specific to the platform, turbine, and choice of parameter ranges, but the demonstrated approach can be applied widely to guide focus in parameter uncertainty.

17 WIND ENERGY↗
DOE OSTI2026

Cooperative effect of local active stresses on the macroscopic contractility of elastic fiber networks

The collective action of actively contractile units embedded in elastic biopolymer networks plays a crucial role in regulating the network's macroscopic mechanical response. Here, in this study, we investigate how the macroscopic boundary stress in model elastic fiber networks depends on the number and nature of embedded contractile units, each exerting an isotropic force dipole, as well as on the bending stiffness of fibers. We find that the macroscopic stress increases nonlinearly with the number of dipoles due to mutual stiffening of initially soft, bending-dominated networks. Using effective medium theory, we relate this enhanced contractility to an increase in the effective average network coordination number due to constraints imposed by the force dipoles. By comparing three distinct force dipole models that differ in their local structures, we demonstrate that the specific manner in which an active unit constrains the network strongly influences the onset and nature of the stiffening transition. Our results highlight that not only the quantity but also the local geometry of force-generating units critically determines the macroscopic mechanical behavior. This framework provides a physical basis for understanding how biological systems—such as molecular motors in the cytoskeleton, or adherent cells in the extracellular matrix—can modulate network-scale nonlinear elastic properties through local tuning of active force-generating units.

Biological and medical sciences↗
DOE OSTI2010

Connection between slow and fast dynamics of molecular liquids around the glass transition

The mean-square displacement (MSD) was measured by neutron scattering at various temperatures and pressures for a number of molecular glass-forming liquids. The MSD is invariant along the glass-transition line at the pressure studied, thus establishing an 'intrinsic' Lindemann criterion for any given liquid. A one-to-one connection between the MSD's temperature dependence and the liquid's fragility is found when the MSD is evaluated on a time scale of {approx}4 ns, but does not hold when the MSD is evaluated at shorter times. The findings are discussed in terms of the elastic model and the role of relaxations, and the correlations between slow and fast dynamics are addressed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

A quantitative assessment of the model form error of friction models across different interface representations for jointed structures

Hysteretic models are widely used to model frictional interactions in joints to recreate experimental behavior. However, it is unclear which models are best suited for fitting or predicting the responses of structures. Here, the present study evaluates 26 friction model/interface representation combinations to quantify the model form error. A Quasi-Static Modal Analysis approach (termed Rayleigh Quotient Nonlinear Modal Analysis) is adopted to calculate the nonlinear system response, and a Multi-Objective Optimization is solved to fit experimental data of the first mode of the Brake-Reuß Beam. Optimized parameters from the first mode are applied to the second and third bending modes to quantify the predictive ability of the models. Formulations for both tracing full hysteresis loops and recreating hysteresis loops from a single loading curve (Masing assumptions) are considered. Smoothly varying models applied to a five patch representation showed the highest flexibility (for fitting mode 1) and good predictive potential (for modes 2 and 3). For a second formulation, which uses 152 frictional elements to represent the interface, the physically motivated spring in series with a Coulomb slip model (elastic dry friction) has high error for fitting mode 1 and performs near the middle for predicting higher modes. For both interface representation, the best fit models are not the most physical, but rather the ones with the most parameters (as expected); however, the more physical models perform somewhat better for predicting the higher modes.

42 ENGINEERING↗
DOE OSTI2016

Development of an empirical model relating permeability and specific stiffness for rough fractures from numerical deformation experiments: Fracture Permeability and Stiffness

The coupling between hydraulic and mechanical properties of porous and fractured geologic media are critical for many geophysical processes and practical applications. Thus, the prediction of this linkage between these properties is broadly important. Here we present a potentially predictive model that links fracture permeability and specific stiffness with empirical coefficients dependent on fracture roughness, correlation length of aperture field, and rock mechanical properties. The equation was developed empirically from results of modeling the deformation and flow through synthetic fractures with aperture fields that follow a normal distribution. The fractures were subjected to increasing normal stress and deformed following an elastic model. Specific stiffness was directly quantified from these numerical deformation experiments. Permeability was estimated through calculation of the local flow field while considering effects of local fracture roughness and tortuosity. The empirical model roughly captures the transition from effective medium to percolation flow regimes with increasing specific stiffness.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2013

A 3D Orthotropic Elastic Continuum Damage Material Model

A three dimensional orthotropic elastic constitutive model with continuum damage is implemented for polymer matrix composite lamina. Damage evolves based on a quadratic homogeneous function of thermodynamic forces in the orthotropic planes. A small strain formulation is used to assess damage. In order to account for large deformations, a Kirchhoff material formulation is implemented and coded for numerical simulation in Sandia’s Sierra Finite Element code suite. The theoretical formulation is described in detail. An example of material parameter determination is given and an example is presented.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Predictive modelling of hydration and mechanical performance of low Ca composite cements: Possibilities and limitations from industrial perspective

We quantified hydration and performance development of composites cements containing metakaolin in the rage of 0–40% and investigated applicability of thermodynamic and micromechanical models to predict the evolution of hydrates assemblage and of performance, respectively. The reaction of metakaolin results in a formation of significant amounts of C-S-H phase for lower metakaolin content and of strätlingite for higher contents. This leads to very dense matrix with finer porosity allowing co-existence of thermodynamically incompatible phases, namely portlandite and strätlingite. Hence, thermodynamic modelling cannot be applied to correctly predict the phase assemblage for the systems with high metakaolin content. By introducing constraints into the thermodynamic model in a form of minimum phase content, the XRD-amorphous phases can be quantified. However, such constraints are only valid for the cases studied and limit the general applicability and predictive power of the model. Elastic properties and strength of mortars are influenced by the metakaolin additions in an opposite trend. It indicates that they are sensitive to different microstructural features/mechanisms. Continuum micromechanic-based model is capable to accurately predict the performance evolution based on the volume of hydrates predicted by the thermodynamic modelling. An additional step required for the accurate prediction - porosity is partitioned based on its saturation at 95% RH – is introduced and substantiated.

36 MATERIALS SCIENCE↗
NASA NTRS1972

Computer program for investigating effects of nonlinear suspension-system elastic properties on parachute inflation loads and motions

A computer program is presented by which the effects of nonlinear suspension-system elastic characteristics on parachute inflation loads and motions can be investigated. A mathematical elastic model of suspension-system geometry is coupled to the planar equations of motion of a general vehicle and canopy. Canopy geometry and aerodynamic drag characteristics and suspension-system elastic properties are tabular inputs. The equations of motion are numerically integrated by use of an equivalent fifth-order Runge-Kutta technique.

Poole, L. R.↗
DOE OSTI2016

Brownian relaxation of an inelastic sphere in air

The procedures that are used to calculate the forces and moments on an aerodynamic body in the rarefied gas of the upper atmosphere are applied to a small sphere of the size of an aerosol particle at sea level. While the gas-surface interaction model that provides accurate results for macroscopic bodies may not be appropriate for bodies that are comprised of only about a thousand atoms, it provides a limiting case that is more realistic than the elastic model. The paper concentrates on the transfer of energy from the air to an initially stationary sphere as it acquires Brownian motion. Individual particle trajectories vary wildly, but a clear relaxation process emerges from an ensemble average over tens of thousands of trajectories. The translational and rotational energies in equilibrium Brownian motion are determined. Empirical relationships are obtained for the mean translational and rotational relaxation times, the mean initial power input to the particle, the mean rates of energy transfer between the particle and air, and the diffusivity. These relationships are functions of the ratio of the particle mass to an average air molecule mass and the Knudsen number, which is the ratio of the mean free path in the air to the particle diameter. The ratio of the molecular radius to the particle radius also enters as a correction factor. The implications of Brownian relaxation for the second law of thermodynamics are discussed.

42 ENGINEERING↗
DOE OSTI2014

Definition of the Semisubmersible Floating System for Phase II of OC4

Phase II of the Offshore Code Comparison Collaboration Continuation (OC4) project involved modeling of a semisubmersible floating offshore wind system as shown below. This report documents the specifications of the floating system, which were needed by the OC4 participants for building aero-hydro-servo-elastic models.

17 WIND ENERGY↗
DOE OSTI2016

Distributed parameter model and experimental validation of a compressive-mode energy harvester under harmonic excitations

This paper presents the modeling and parametric analysis of the recently proposed nonlinear compressive-mode energy harvester (HC-PEH) under harmonic excitation. Both theoretical and experimental investigations are performed in this study over a range of excitation frequencies. Specially, a distributed parameter electro-elastic model is analytically developed by means of the energy-based method and the extended Hamilton’s principle. An analytical formulation of bending and stretching forces are derived to gain insight on the source of nonlinearity. Furthermore, the analytical model is validated against with experimental data and a good agreement is achieved. Both numerical simulations and experiment illustrate that the harvester exhibits a hardening nonlinearity and hence a broad frequency bandwidth, multiple coexisting solutions and a large-amplitude voltage response. Using the derived model, a parametric study is carried out to examine the effect of various parameters on the harvester voltage response. It is also shown from parametric analysis that the harvester’s performance can be further improved by selecting the proper length of elastic beams, proof mass and reducing the mechanical damping.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1976

Correlation of full-scale helicopter rotor performance in air with model-scale Freon data

An investigation was conducted in a transonic dynamics tunnel to measure the performance of a 1/5 scale model helicopter rotor in a Freon atmosphere. Comparisons were made between these data and full scale data obtained in air. Both the model and full scale tests were conducted at advance ratios between 0.30 and 0.40 and advancing tip Mach numbers between 0.79 and 0.95. Results show that correlation of model scale rotor performance data obtained in Freon with full scale rotor performance data in air is good with regard to data trends. Mach number effects were found to be essentially the same for the model rotor performance data obtained in Freon and the full scale rotor performance data obtained in air. It was determined that Reynolds number effects may be of the same magnitude or smaller than rotor solidity effects or blade elastic modeling in rotor aerodynamic performance testing.

Yeager, W. T., Jr.↗
NASA NTRS1993

Constraints on interseismic deformation at Japan trench from VLBI data

Space geodetic data from very long baseline interferometry (VLBI) was used to estimate velocity relative to the plate interiors of two sites on the deforming leading edge at the Japan trench. Elastic models of interseismic deformation and results obtained were used to put constraints on the slip rate along the main thrust of the Japan subduction zone. Observed velocities reflect the sum of permanent west-northwest shortening in Honshu, elastic deformation due to locking of the main thrust fault at the Japan trench, and deformation associated with the subducting Phillipine plate. These velocities limit the locked segment of the main thrust at the Japan trench to 27 km vertically and 100 km along the dip. This indicates that the main Pacific plate thrust fault is not strongly coupled and probably does not generate strong earthquakes.

Argus, Donald F.↗
DOE OSTI2021

Efficient Failure Probability Calculations and Modeling Interface Debonding in TRISO Particles with BISON

This document constitutes completion of the NEAMS milestone, which is titled Demonstrate efficiency improvements in TRISO failure probability calculations and the effect of debonding on fission product transport in TRISO fuel particles in Bison . In this report we present the development of the following: (1) the anisotropic elasticity model of pyrolytic carbon and extended material models using the local coordinate system for aspherical particle geometry; (2) a capability of modeling interface debonding in TRISO particles; and (3) an efficient high-fidelity approach to calculate the failure probability of TRISO particles. Benchmark problems and AGR-2 tests are used to verify the models and demonstrate these new capabilities.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

Modeling the Effective Elasticity of Anisotropic Porous Materials

The development and optimization of composite materials designed for thermal protection of NASA’s spacecraft require understanding their physical response to high-enthalpy environments. To predict their macro-scale properties and behavior, high-fidelity 3D simulations are performed at the microscale on realistic representations of these composites. The digital microstructures are generated either synthetically or through X-ray micro-computed tomography reconstructions. One of the main challenges in the prediction of the structural response of heatshield materials is the computation of the effective elasticity of the fibrous composite, as well as the understanding of the deformation and stresses generated at the microscale. These are driven by the fiber layout within the microstructure and the distribution of the infused matrix. In this effort, the micro-mechanical linear elastic behavior of fibrous ablators is modeled using a numerical method based on the Multi-Point Stress Approximation (MPSA) finite volume scheme, a generalization of the more commonly used Multi-Point Flux Approximation (MPFA) that was presented at the 10th Ablation Workshop. To predict the behavior of fibrous and woven architectures, algorithms that compute the local fiber orientation are used. The implementation of the MPSA was verified using analytical solutions, engineering test cases, and compared against legacy Finite Element Analysis (FEA) software. The stress analysis models were then applied to real geometries used by NASA in thermal protection systems such as fibrous preforms and woven materials and the results were compared to experimental data.

Elasticity↗
DOE OSTI2015

An Elastic Plastic Contact Model with Strain Hardening for the LAMMPS Granular Package

The following details the implementation of an analytical elastic plastic contact model with strain hardening for normal im pacts into the LAMMPS granular package. The model assumes that, upon impact, the co llision has a period of elastic loading followed by a period of mixed elastic plas tic loading, with contributions to each mechanism estimated by a hyperbolic seca nt weight function. This function is implemented in the LAMMPS source code as the pair style gran/ep/history. Preliminary tests, simulating the pouring of pure nickel spheres, showed the elastic/plastic model took 1.66x as long as similar runs using gran/hertz/history.

36 MATERIALS SCIENCE↗
NASA NTRS1976

Computer simulation of earthquakes

Two computer simulation models of earthquakes were studied for the dependence of the pattern of events on the model assumptions and input parameters. Both models represent the seismically active region by mechanical blocks which are connected to one another and to a driving plate. The blocks slide on a friction surface. In the first model elastic forces were employed and time independent friction to simulate main shock events. The size, length, and time and place of event occurrence were influenced strongly by the magnitude and degree of homogeniety in the elastic and friction parameters of the fault region. Periodically reoccurring similar events were frequently observed in simulations with near homogeneous parameters along the fault, whereas, seismic gaps were a common feature of simulations employing large variations in the fault parameters. The second model incorporated viscoelastic forces and time-dependent friction to account for aftershock sequences. The periods between aftershock events increased with time and the aftershock region was confined to that which moved in the main event.

Cohen, S. C.↗
DOE OSTI2022

Equation of state and strength of diamond in high-pressure ramp loading

Diamond is used extensively as a component in high-energy-density experiments, but existing equation of state (EOS) models do not capture its observed response to dynamic loading. In particular, in contrast with first-principles theoretical models of equilibrium multiphase EOS, no solid-solid phase changes have been detected, and no general-purpose EOS models match the measured ambient isotherm. We have performed density functional theory (DFT) calculations of the diamond phase to ~10 TPa, well beyond its predicted range of thermodynamic stability, and used these results as the basis of a Mie-Grüneisen EOS. We also performed DFT calculations of the elastic moduli, and calibrated an algebraic elasticity model for use in simulations. We then estimated the flow stress of diamond by comparison with the stress-density relation measured experimentally in ramp-loading experiments. Furthermore, the resulting constitutive model allows us to place a constraint on the Taylor-Quinney factor (the fraction of plastic work converted to heat) from the observation that diamond does not melt on ramp compression.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗