SEARCH · Engineering Papers

Results for “continuum 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 523 records · Page 29

All sources NASA NTRS DOE OSTI

NASA NTRS1991

A quick accurate model of nozzle backflow

Backflow from nozzles is a major source of contamination on spacecraft. If the craft contains any exposed high voltages, the neutral density produced by the nozzles in the vicinity of the craft needs to be known in order to assess the possibility of Paschen breakdown or the probability of sheath ionization around a region of the craft that collects electrons for the plasma. A model for backflow has been developed for incorporation into the Environment-Power System Analysis Tool (EPSAT) which quickly estimates both the magnitude of the backflow and the species makeup of the flow. By combining the backflow model with the Simons (1972) model for continuum flow it is possible to quickly estimate the density of each species from a nozzle at any position in space. The model requires only a few physical parameters of the nozzle and the gas as inputs and is therefore ideal for engineering applications.

Kuharski, R. A.↗

NASA NTRS2017

Multilayer Insulation Ascent Venting Model

The thermal and venting transient experienced by tank-applied multilayer insulation (MLI) in the Earth-to-orbit environment is very dynamic and not well characterized. This new predictive code is a first principles-based engineering model which tracks the time history of the mass and temperature (internal energy) of the gas in each MLI layer. A continuum-based model is used for early portions of the trajectory while a kinetic theory-based model is used for the later portions of the trajectory, and the models are blended based on a reference mean free path. This new capability should improve understanding of the Earth-to-orbit transient and enable better insulation system designs for in-space cryogenic propellant systems.

Tramel, R. W.↗

DOE OSTI2021

Screening of metal complexes and organic solvents using the COSMOSAC-LANL model to enhance the energy density in a non-aqueous redox flow cell: an insight into the solubility

In this paper, we have proposed a first-principles methodology to screen transition metal complexes against a particular organic solvent and organic solvents against a particular transition metal complex based on their solubility information without the knowledge of heat of fusion and melting temperature. The energy density of a non-aqueous redox flow cell directly depends on the solubility of the redox active species in the non-aqueous medium. Here, we have used the “COSMOSAC-LANL” activity coefficient model which is based on first-principles COSMO calculations where the microscopic information is passed to the macroscopic world via a dielectric continuum solvation model, followed by a post-statistical thermodynamic treatment of the self-consistent properties of the solute particle to calculate the solubility. To model the activity coefficient at infinite dilution for the binary mixtures, a 3-suffix Margules (3sM) function is introduced for the quantitative estimation of the asymmetric interactions and, for the combinatorial term, the Staverman–Guggenheim (SG) form is used. The new activity coefficient model is separately called the “LANL” activity coefficient model. The metal complex and the organic solvent have been treated as a simple binary mixture. The present model has been applied to a set of 14 different organic solvents and 16 different transition metal complexes. Using the new LANL activity coefficient model in combination with the ADF-COSMOSAC-2013 model, we have shown how one can improve the solubility of a transition metal complex in an organic solvent. We applied our model to screen 84 binary mixtures to predict the compatible pair of redox active species and organic solvent to increase the energy density. The solvation mechanism of the transition metal complexes in the organic solvents was obtained using the new model. The results have been compared with the experimental and theoretical results where they are available.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1994

A methodology to predict damage initiation, damage growth and residual strength in titanium matrix composites

In this research, a methodology to predict damage initiation, damage growth, fatigue life, and residual strength in titanium matrix composites (TMC) is outlined. Emphasis was placed on micromechanics-based engineering approaches. Damage initiation was predicted using a local effective strain approach. A finite element analysis verified the prevailing assumptions made in the formulation of this model. Damage growth, namely, fiber-bridged matrix crack growth, was evaluated using a fiber bridging (FB) model which accounts for thermal residual stresses. This model combines continuum fracture mechanics and micromechanics analyses yielding stress-intensity factor solutions for fiber-bridged matrix cracks. It is assumed in the FB model that fibers in the wake of the matrix crack are idealized as a closure pressure, and an unknown constant frictional shear stress is assumed to act along the debond length of the bridging fibers. This frictional shear stress was used as a curve fitting parameter to the available experimental data. Fatigue life and post-fatigue residual strength were predicted based on the axial stress in the first intact 0 degree fiber calculated using the FB model and a three-dimensional finite element analysis.

Bakuckas, J. G., Jr.↗

DOE OSTI2023

A sensitivity analysis of twinning crystal plasticity finite element model using single crystal and poly crystal Zircaloy

The popularity of crystal plasticity finite element method (CPFEM) models is increasing due to their ability to predict the mechanical response of crystalline materials such as metals and metal alloys more accurately than traditional continuum mechanics models. This is since the crystal plasticity models consider the effect of atomic structure, microstructural morphology, and properties of individual grains. These CPFEM models use a large number of material parameters in order to capture the mesoscale physics which comes with the downside of the tedious calibration process. In this paper, a CPFEM code was developed to include the twinning induced grain reorientation and subsequent crystallographic slip for HPC material. The developed code is incorporated in a large-scale, parallelized nonlinear solver WARP3D. Further, a sensitivity analysis with respect to 22 material parameters was then conducted using single crystal and polycrystal representative volume element (RVE) of Zircaloy material. Loading was applied along five different crystallographic orientations for single crystal RVE and along three directions namely, rolling (RD), transverse (TD), and normal (ND) direction for polycrystal RVE. Results obtained from the sensitivity analysis were used for the calibration of material parameters for Zircaloy. Finally, developed code along with calibrated material parameters was used to investigate the effect of the hydride phase formation in Zircaloy which is a typical case observed for nuclear applications. It was found that the volume fraction of the hydride phase has a significant impact on the mechanical properties of Zircaloy.

36 MATERIALS SCIENCE↗

NASA NTRS2003

Kinematics and the implementation of an elephant's trunk manipulator and other continuum style robots

Traditionally, robot manipulators have been a simple arrangement of a small number of serially connected links and actuated joints. Though these manipulators prove to be very effective for many tasks, they are not without their limitations, due mainly to their lack of maneuverability or total degrees of freedom. Continuum style (i.e., continuous "back-bone") robots, on the other hand, exhibit a wide range of maneuverability, and can have a large number of degrees of freedom. The motion of continuum style robots is generated through the bending of the robot over a given section; unlike traditional robots where the motion occurs in discrete locations, i.e., joints. The motion of continuum manipulators is often compared to that of biological manipulators such as trunks and tentacles. These continuum style robots can achieve motions that could only be obtainable by a conventionally designed robot with many more degrees of freedom. In this paper we present a detailed formulation and explanation of a novel kinematic model for continuum style robots. The design, construction, and implementation of our continuum style robot called the elephant trunk manipulator is presented. Experimental results are then provided to verify the legitimacy of our model when applied to our physical manipulator. We also provide a set of obstacle avoidance experiments that help to exhibit the practical implementation of both our manipulator and our kinematic model. c2003 Wiley Periodicals, Inc.

Non-NASA Center↗

NASA NTRS1996

Micro and Macro Segregation in Alloys Solidifying with Equiaxed Morphology

To understand macro segregation formation in Al-Cu alloys, experiments were run under terrestrial gravity (1g) and under low gravity during parabolic flights (10(exp -2) g). Alloys of two different compositions (2% and 5% Cu) were solidified at two different cooling rates. Systematic microscopic and SEM observations produced microstructural and segregation maps for all samples. These maps may be used as benchmark experiments for validation of microstructure evolution and segregation models. As expected, the macro segregation maps are very complex. When segregation was measured along the central axis of the sample, the highest macro segregation for samples solidified at 1g was obtained for the lowest cooling rate. This behavior is attributed to the longer time available for natural convection and shrinkage flow to affect solute redistribution. In samples solidified under low-g, the highest macro-segregation was obtained at the highest cooling rate. In general, low-gravity solidification resulted in less segregation. To explain the experimental findings, an analytical (Flemings-Nereo) and a numerical model were used. For the numerical model, the continuum formulation was employed to describe the macroscopic transports of mass, energy, and momentum, associated with the microscopic transport phenomena, for a two-phase system. The model proposed considers that liquid flow is driven by thermal and solutal buoyancy, and by solidification shrinkage. The Flemings-Nereo model explains well macro segregation in the initial stages of low-gravity segregation. The numerical model can describe the complex macro segregation pattern and the differences between low- and high-gravity solidification.

Stefanescu, Doru M.↗

NASA NTRSDate not supplied

Sensitivity of Single-Scattering Properties to Precipitation Particle Meltwater Geometry for GPM Passive Microwave and Radar Remote Sensing Applications

The objective of the current work is to describe the microwave single-scattering properties of partially melted ice-phase precipitation particles using physically-based computational methods, taking into full account the varied geometries of the initially “dry” ice particles and the distributions of liquid water that develop during the melting process. Ultimately, the bulk properties of ensembles of these partially-melted particles will be included in “scattering tables” to support radar and combined radar-radiometer remote sensing of precipitation. A meshless Lagrangian melting procedure (snowMELT) based on smoothed-particle hydrodynamics is applied to both spherical and finely-structured ice particle models to describe the full evolution of the particles from dry ice particles to liquid drops. The melting of spherical ice particles using snowMELT is compared to an alternative continuum physics model to validate melt times and internal thermodynamics. The discrete dipole approximation is then utilized to calculate the single-scattering properties of different mixed-phase particles throughout the melting process. The sensitivities of particle single-scattering properties to thermodynamic and hydrodynamic assumptions in snowMELT are explored, and the implications for combined radar-radiometer precipitation remote sensing from GPM are discussed.

William S Olson↗

NASA NTRS1999

The Role of Self-Organized Criticality in the Substorm Phenomenon and its Relation to Localized Reconnection in the Plasma Sheet

Recent observations of the magnetotail plasma sheet have shown it to be a dynamic and turbulent region. Research has found strong turbulence in the plasma sheet at approximately 20 Earth's Radius tailward of Earth; the turbulence is observed at all activity levels. The existence of strong turbulence in the plasma sheet in the region associated with substorm onset might be thought difficult to reconcile with the coherence and repeatability of the substorm cycle. We review a variety of evidence that strongly suggests the magnetotail is driven, through magnetic flux transfer, into a state of "self-organized criticality" (SOC). It is an important property of physical systems that evolve into SOC that they self-organize into a unique, global dynamic state. This global state is inevitable, and repeatable. In this state, however, small-spatiotemporal-scale system fluctuations are unpredictable and can be only described statistically. This is the basis, we think, for the global coherence and repeatability of the substorm phenomenon in the turbulent plasma sheet. At, or near, substorm onset the plasma sheet can be described by a global SOC state containing significant small scale turbulence. In several recent studies, "sandpile" models were driven into SOC and then shown to reproduce various measures of substorm activity. We discuss the plasma physical foundation of these sandpile models. The evolution of simple continuum plasma sheet models into SOC-like states of many small reconnection events in the turbulent plasma sheet under the will be demonstrated. We view the substorm phenomenon as an avalanche assumption that the plasma sheet is in a SOC state.

Klimas, A. J.↗

NASA NTRS2012

Simulating the Response of a Composite Honeycomb Energy Absorber: Dynamic Crushing of Components and Multi-Terrain Impacts - Part 1

This paper describes the experimental and analytical evaluation of an externally deployable composite honeycomb structure that is designed to attenuate impact energy during helicopter crashes. The concept, designated the Deployable Energy Absorber (DEA), utilizes an expandable Kevlar (Registered Trademark) honeycomb to dissipate kinetic energy through crushing. The DEA incorporates a unique flexible hinge design that allows the honeycomb to be packaged and stowed until needed for deployment. Experimental evaluation of the DEA included dynamic crush tests of multi-cell components and vertical drop tests of a composite fuselage section, retrofitted with DEA blocks, onto multi-terrain. Finite element models of the test articles were developed and simulations were performed using the transient dynamic code, LSDYNA (Registered Trademark). In each simulation, the DEA was represented using shell elements assigned two different material properties: Mat 24, an isotropic piecewise linear plasticity model, and Mat 58, a continuum damage mechanics model used to represent laminated composite fabrics. DEA model development and test-analysis comparisons are presented.

Jackson, Karen E.↗

NASA NTRS2014

Simulating the Impact Response of Composite Airframe Components

In 2010, NASA Langley Research Center obtained residual hardware from the US Army's Survivable Affordable Repairable Airframe Program (SARAP). The hardware consisted of a composite fuselage section that was representative of the center section of a Black Hawk helicopter. The section was fabricated by Sikorsky Aircraft Corporation and designated the Test Validation Article (TVA). The TVA was subjected to a vertical drop test in 2008 to evaluate a tilting roof concept to limit the intrusion of overhead mass items, such as the rotor transmission, into the fuselage cabin. As a result of the 2008 test, damage to the hardware was limited primarily to the roof. Consequently, when the post-test article was obtained in 2010, the roof area was removed and the remaining structure was cut into six different types of test specimens including: (1) tension and compression coupons for material property characterization, (2) I-beam sections, (3) T-sections, (4) cruciform sections, (5) a large subfloor section, and (6) a forward framed fuselage section. In 2011, NASA and Sikorsky entered into a cooperative research agreement to study the impact responses of composite airframe structures and to evaluate the capabilities of the explicit transient dynamic finite element code, LS-DYNA®, to simulate these responses including damage initiation and progressive failure. Finite element models of the composite specimens were developed and impact simulations were performed. The properties of the composite material were represented using both a progressive in-plane damage model (Mat 54) and a continuum damage mechanics model (Mat 58) in LS-DYNA. This paper provides test-analysis comparisons of time history responses and the location and type of damage for representative I-beam, T-section, and cruciform section components.

Jackson, Karen E.↗

DOE OSTI2020

Modeling Single and Two-Phase Transport in Thin Porous Layers Using a Composite Continuum-Pore Network Formulation

In this work, a composite continuum-pore network formulation is presented to model single and two-phase transport thin porous layers, such as gas diffusion layers in polymer electrolyte fuel cells (PEFCs) and active electrodes in redox flow batteries (RFBs). The formulation can be integrated into CFD codes, thus combining the ease of implementation of continuum-based modeling and the computational power of pore-network modeling. The composite model includes a control volume (CV) mesh at the layer scale, which embeds a cubic pore network [1,2]. The pore-network model is used to determine analytically local anisotropic effective transport properties (local effective diffusivity and permeability), which are mapped onto the CV mesh to simulate transport in the porous transport layer [3,4]. Good agreement is found between the predicted global effective transport properties (global effective diffusivity and permeability) under dry and wet conditions and previous experimental data reported in the literature for Toray TGP-H series carbon paper. Water saturation distributions are also compared with results obtained using X-ray computed tomography [4]. [1] P.A. García-Salaberri, I.V. Zenyuk, J.T. Gostick, A.Z. Weber, Modeling Gas Diffusion Layer in Polymer Electrolyte Fuel Cells Using a Continuum-Based Pore-Network Formulation, ECS Trans. 97 (2020) 615. [2] P.A. García-Salaberri, Modeling diffusion and convection in thin porous transport layers using a composite continuum-network model: Application to gas diffusion layers in polymer electrolyte fuel cells, Int J. Heat Mass Transf. (2020), submitted. [3] P.A. García-Salaberri, J.T. Gostick, G. Hwang, A.Z. Weber, M. Vera, Effective diffusivity in partially-saturated carbon-fiber gas diffusion layers: Effect of local saturation and application to macroscopic continuum models, J. Power Sources 296 (2015) 440–453. [4] P.A. García-Salaberri, G. Hwang, M. Vera, A.Z. Weber, J.T. Gostick, Effective diffusivity in partially-saturated carbon-fiber gas diffusion layers: Effect of through-plane saturation distribution, Int. J. Heat Mass Transf. 86 (2015) 319–333.

Garcia-Salaberri, Pablo↗

NASA NTRS2010

Study of Plume Impingement Effects in the Lunar Lander Environment

Plume impingement effects from the descent and ascent engine firings of the Lunar Lander were analyzed in support of the Lunar Architecture Team under the Constellation Program. The descent stage analysis was performed to obtain shear and pressure forces on the lunar surface as well as velocity and density profiles in the flow field in an effort to understand lunar soil erosion and ejected soil impact damage which was analyzed as part of a separate study. A CFD/DSMC decoupled methodology was used with the Bird continuum breakdown parameter to distinguish the continuum flow from the rarefied flow. The ascent stage analysis was performed to ascertain the forces and moments acting on the Lunar Lander Ascent Module due to the firing of the main engine on take-off. The Reacting and Multiphase Program (RAMP) method of characteristics (MOC) code was used to model the continuum region of the nozzle plume, and the Direct Simulation Monte Carlo (DSMC) Analysis Code (DAC) was used to model the impingement results in the rarefied region. The ascent module (AM) was analyzed for various pitch and yaw rotations and for various heights in relation to the descent module (DM). For the ascent stage analysis, the plume inflow boundary was located near the nozzle exit plane in a region where the flow number density was large enough to make the DSMC solution computationally expensive. Therefore, a scaling coefficient was used to make the DSMC solution more computationally manageable. An analysis of the effectiveness of this scaling technique was performed by investigating various scaling parameters for a single height and rotation of the AM. Because the inflow boundary was near the nozzle exit plane, another analysis was performed investigating three different inflow contours to determine the effects of the flow expansion around the nozzle lip on the final plume impingement results.

Marichalar, Jeremiah↗

NASA NTRS1992

Decoupled predictions of radiative heating in air using a particle simulation method

A particle simulation technique with decoupled radiation is used to estimate the radiative emission along the stagnation streamline and the radiative heating at the stagnation point of a blunt-nosed vehicle entering the earth's atmosphere at hypersonic speed. The direct simulation Monte Carlo (DSMC) method is used to compute the fluid mechanics of the weakly ionized flow. The radiative emission is computed using the NEQAIR computer code with the radiation decoupled from the flowfield solution. It is concluded that the new models make it possible to decrease the predicted total radiative heating at the stagnation point of the vehicle by a factor of 15. A comparison of the DSMC approach with a continuum flow model shows that the total predicted radiative heating estimates agree within a factor of 2.

Boyd, Iain D.↗

DOE OSTI2020

A continuum damage mechanics (CDM) based Wilshire model for creep deformation, damage, and rupture prediction

In this report, a new constitutive model is derived that combines the Wilshire equations with continuum damage mechanics (CDM) to enable the long-term prediction of creep deformation, damage, and rupture. The Wilshire equations have been demonstrated to accurately extrapolate the stress-rupture, minimum-creep-strain-rate, and time-to-strain of various alloys across decade of life. Recently, the time-to-creep-strain equation has been rearranged to predict creep deformation curves; however, the resulting equation is difficult to calibrate and does not track damage. The CDM-based Sinh constitutive model accurately predicts creep deformation, damage, and stress-rupture for a variety of alloys; however, it lacks an explicit description of stress and temperature co-dependence which limits interpolation and extrapolation ability. In this study, the Wilshire and the CDM-based Sinh equations are combined to create the “WCS” model. The WCS model combines the best features of each model while eliminating their deficiencies. Experimental data for alloy P91 is gathered and the WCS model is calibrated to the data. The WCS model accurately predicts the stress rupture, minimum-creep-strain-rate, creep deformation, and damage of alloy P91 across the experimental stress and temperature range as well as consistency in parametric simulations.

36 MATERIALS SCIENCE↗

NASA NTRS2001

The Constitutive Modeling of Thin Films with Randon Material Wrinkles

Material wrinkles drastically alter the structural constitutive properties of thin films. Normally linear elastic materials, when wrinkled, become highly nonlinear and initially inelastic. Stiffness' reduced by 99% and negative Poisson's ratios are typically observed. This paper presents an effective continuum constitutive model for the elastic effects of material wrinkles in thin films. The model considers general two-dimensional stress and strain states (simultaneous bi-axial and shear stress/strain) and neglects out of plane bending. The constitutive model is derived from a traditional mechanics analysis of an idealized physical model of random material wrinkles. Model parameters are the directly measurable wrinkle characteristics of amplitude and wavelength. For these reasons, the equations are mechanistic and deterministic. The model is compared with bi-axial tensile test data for wrinkled Kaptong(Registered Trademark) HN and is shown to deterministically predict strain as a function of stress with an average RMS error of 22%. On average, fitting the model to test data yields an RMS error of 1.2%

Murphey, Thomas W.↗

NASA NTRS1995

Multiwavelength monitoring of the BL Lacertae object PKS 2155-304. 4: Multiwavelength analysis

Simultaneous X-ray, ultraviolet, optical, infrared, and radio monitoring data were used to test and constrain models of continuum emission from the BL Lacertae object PKS 2155-304. Intensively sampled ultraviolet and soft X-ray light curves showed a clear temporal correlation with the X-rays leading the ultraviolet by 2-3 hr. This lag was found to be significantly different from zero after an exhaustive comparison of four different techniques for measuring temporal correlations. Variations in the ultraviolet trough optical wave bands were also strongly correlated, with no measurable lag down to limiting timescales of approximately less than 1-2 hr. This strong correlation extends to the near-infrared, but the less intensive sampling precludes measurement of any lag beyomnd an upper limit of approximately less than 1 day. These lags and limits of the order of hours are much shorter than most rapid observed single-band variations. Because of the very sparse radio sampling, it was not possible to measure quantitatively the correlation and lag with shorter wavelengths, but the data do suggest that the radio may lag the optical/ultraviolet by approximately 1 week, with longer delays and weaker variations to longer radio wavelengths. The epoch-folding Q(exp 2) statistic was used to test for periodicity, and no evidence for strict or quasi-periodicity was found in any of the light curves. Because they lead the lower frequencies, the soft X-rays (approximately less than 1 keV) cannot arise from synchrotron self-Compton scattering. These results also rule out the accretion disk model, which predicts a measurable lag between ultraviolet/optical wavelength bands and a correlation between hardness and brightness, neither of which were seen. They are consistent with the entire radio through X-ray continuum arising from direct synchrotron emission from a relativistic jet. However, the tapered jet model, in which the X-ray emission is produced closer in, has problems explaining the magnitude of the ultraviolet/X-ray lag, because the X-ray-emitting electrons have very short lifetimes (t(sub 1/2) much less than 1 s). The result that the lag is much smaller than the variability timescale suggests instead that the radiation may be produced in a flattened region such as a shock front.

Edelson, R.↗

DOE OSTI2024

A B-spline based gradient-enhanced micropolar implicit material point method for large localized inelastic deformations

The quasi-brittle response of cohesive-frictional materials in numerical simulations is commonly represented by softening plasticity or continuum damage models, either individually or in combination. However, classical models, particularly when coupled with non-associated plasticity, often suffer from ill-posedness and a lack of objectivity in numerical simulations. Moreover, the performance of the finite element method significantly degrades in simulations involving finite strains when mesh distortion reaches excessive levels. This represents a challenge for modeling cohesive-frictional materials, given their tendency to experience strongly localized deformations, such as those occurring during shear band dominated failure. Hence, accurate modeling of the response of cohesive-frictional solids is a demanding task. To address these challenges, we present an extension of the material point method (MPM) for the unified gradient-enhanced micropolar continuum, aiming at the analysis of finite localized inelastic deformations in cohesive-frictional materials. The generalized gradient-enhanced micropolar continuum formulation is employed to tackle challenges related to localization and softening material behavior, while the MPM addresses issues arising from excessive deformations. The method utilizes a B-spline formulation for the rigid background mesh to mitigate the well-known cell crossing errors of the MPM. To demonstrate the performance of the method, 2D and 3D numerical studies on localized failure in sandstone in plane strain compression and triaxial extension tests are presented. A comparison with finite element results confirms the suitability of the formulation. Moreover, an efficient numerical implementation of the formulation is presented, and it is demonstrated that the additional MPM specific overhead is negligible.

B-spline↗