SEARCH · Engineering Papers

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

SourceResearch connectionFrom yearThrough year

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

At least 1,585 records · Page 88

All sources NASA NTRS DOE OSTI

NASA NTRS2016

Modeling Fatigue Damage Onset and Progression in Composites Using an Element-Based Virtual Crack Closure Technique Combined With the Floating Node Method

A new methodology is proposed to model the onset and propagation of matrix cracks and delaminations in carbon-epoxy composites subject to fatigue loading. An extended interface element, based on the Floating Node Method, is developed to represent delaminations and matrix cracks explicitly in a mesh independent fashion. Crack propagation is determined using an element-based Virtual Crack Closure Technique approach to determine mixed-mode energy release rates, and the Paris-Law relationship to obtain crack growth rate. Crack onset is determined using a stressbased onset criterion coupled with a stress vs. cycle curve and Palmgren-Miner rule to account for fatigue damage accumulation. The approach is implemented in Abaqus/Standard® via the user subroutine functionality. Verification exercises are performed to assess the accuracy and correct implementation of the approach. Finally, it was demonstrated that this approach captured the differences in failure morphology in fatigue for two laminates of identical stiffness, but with layups containing θdeg plies that were either stacked in a single group, or distributed through the laminate thickness.

De Carvalho, Nelson V.↗

NASA NTRSDate not supplied

Status of Mars Retropropulsion Testing in the Langley Unitary Plan Wind Tunnel

Future Mars human landings will be enabled by a powered descent phase starting at supersonic conditions, something which has never been done before on a Mars mission. Significant aerosciences challenges exist due to jet interactions between the retrorocket engine plumes, freestream flow, and vehicle that will affect the aerodynamic behavior during powered descent. Historically, wind tunnel tests have been used to study the interactions with inert gas exhaust simulants in place of rocket engines. On the computational side, flowfield simulations have been completed at full-scale conditions, but the available ground and flight data are not appropriate for calibrating computational uncertainties for aerodynamic interference on proposed Mars descent vehicles, due to insufficient data, dissimilar vehicle geometries, and disparate operating conditions. A wind tunnel test has been designed to begin addressing powered descent aerodynamics risks for large-scale human Mars entry concepts and to identify gaps in computational predictive capabilities. The test will be conducted in the NASA Langley Unitary Plan Wind Tunnel and is designed with improvements in model design and data products over past tests. The test campaign will be run using sub-scale model geometries derived from NASA powered descent reference vehicles: a blunt low lift-to-drag vehicle and a more slender geometry that generates higher unpowered lift. Both models have been fabricated and are ready for testing. The blunt model is equipped with the flexibility to examine the effects of nozzle pointing direction, number, location, size, and area ratio. The main measurements are heatshield aerodynamic interference forces and moments with a custom flow-through balance, discrete and distributed heatshield pressure, and high-speed flowfield visualization. This paper covers the test objectives, facility, models and instrumentation, and planned test matrix.

Mars↗

NASA NTRSDate not supplied

Status of Mars Retropropulsion Testing in the Langley Unitary Plan Wind Tunnel

Supersonic Retropropulsion↗

NASA NTRS2021

Polarimetric sensitivity of light-absorbing carbonaceous aerosols over ocean: A theoretical assessment

Visible-light-absorbing carbonaceous aerosols within the boundary layer affect the radiance and polarization states of the radiation at the top of the atmosphere. Remote sensing from suborbital and satellite-based platforms utilizes these radiance and polarization signals to retrieve the key properties of these aerosols. Recent retrieval algorithms have shown a progressive trend toward including multi-angular and multi-spectral polarimetric measurements to produce better retrieval accuracy in comparison to those using measurements based on a single viewing angle. Here, we perform a theoretical investigation of the top of atmosphere (TOA) radiance-related reflectance factor (bidirectional reflectance factor (BRF)) and the two types of polarimetry-related factors (polarized bidirectional reflectance factor (pBRF) and the degree of linear polarization (DoLP)) for different types of atmospheric light-absorbing carbonaceous aerosols as a function of particle size distribution. We selected three polarimetric bands corresponding to those utilized by NASA's Airborne Multiangle SpectroPolarimetric Imager (AirMSPI)—near-UV (470 nm), visible (660 nm), and near-infrared (865 nm)—for our simulations which were performed over ocean surface using the successive order of scattering (SOS) algorithm coupled to a Lorenz-Mie aerosol optics model. The analysis of particle phase matrix elements indicates a close relationship between the angular dependencies of DoLP and associated phase matrix components at the shortest polarimetric band (470 nm). Using Jacobian analysis, we find that the radiance- and polarimetry-related reflectance factors of weakly light-absorbing aerosols, such as brown carbon, are more sensitive to changes in particle size and imaginary refractive index in comparison with those of black carbon, which is strongly light-absorbing. Our results suggest that the DoLP data could be used by future retrieval algorithms for reliably estimating microphysical properties of absorbing carbonaceous aerosols with imaginary refractive index less than 0.4.

refractive index↗

DOE OSTI2016

JASMINE: a fuel rod thermal-mechanical performance's code

JASMINE, which is developed by China General Nuclear Power Corporation (CGNPC), is a code for calculating thermal-mechanical response of PWR fuel rods. The code models many physical phenomena, including heat conduction through the fuel and cladding to the coolant, cladding elastic and plastic deformation, fuel elastic deformation, fuel swelling and densification,fuel-cladding mechanical interaction, fission gas release from the fuel, rod internal pressure, cladding oxidation and hydriding. Just as most fuel rod codes, JASMINE, which is a quasi-two dimensional code, includes a thermal module, a mechanical module, a fission gas release (FGR) module, and an internal pressure module. For thermal module, the time dependent thermal conduction equation of axisymmetric cylindrical geometry, Halden's fuel thermal conductivity, thermal conductivity integration method and multi-dimension Newton-Raphson procedure are used; for the mechanical module, generalized plane strain hypothesis, cladding-fuel axial lockup model, cladding thermal stress and transfer matrix method are applied; for fission gas release module, the time dependent diffusion equation for a spherical grain, Speight-Turnbull grain surface boundary condition, and Forsberg-Massih's method are used; for the internal pressure module, the ideal gas state equation is used. The code, which is on validation stage, will be used for fuel rod design of nuclear power plants in China in a near future. So far, the preliminary validation process for pellet centerline temperature and fission gas release, which is based on fuel rod test data in FUMEX-III project, has finished. According to the validation results, JASMINE has the ability to reasonably predicts fuel behavior and performance in-pile. JASMINE has good prediction in fuel temperature, and good prediction on FGR as a whole. However, for FGR in the range 1% to 15%, JASMINE over-predicts FGR; for other range, JASMINE underestimates FGR. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2023

Complementarity and the unitarity of the black hole S -matrix

Recently, Akers et al. proposed a non-isometric holographic map from the interior of a black hole to its exterior. Within this model, we study properties of the black hole S-matrix, which are in principle accessible to observers who stay outside the black hole. Specifically, we investigate a scenario in which an infalling agent interacts with radiation both outside and inside the black hole. Because the holographic map involves postselection, the unitarity of the S-matrix is not guaranteed in this scenario, but we find that unitarity is satisfied to very high precision if suitable conditions are met. If the internal black hole dynamics is described by a pseudorandom unitary transformation, and if the operations performed by the infaller have computational complexity scaling polynomially with the black hole entropy, then the S-matrix is unitary up to corrections that are superpolynomially small in the black hole entropy. Furthermore, while in principle quantum computation assisted by postselection can be very powerful, we find under similar assumptions that the S-matrix of an evaporating black hole has polynomial computational complexity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Neutrinoless Double-Beta Decay Investigations of 82Se Using Three Shell Model Hamiltonians

Neutrinoless double-beta decay is considered one of the most promising processes that would help clarify some of the symmetry-breaking problems in our understanding of the observable universe. Recent studies of neutrinoless double-beta decay matrix elements have employed statistical approaches based on modified shell model effective Hamiltonians for 48Ca (Phys. Rev. C 106, 054302 (2022)) and 136Xe (Phys. Rev. C 107, 045501 (2023)). The analyses rely on inducing perturbations in the starting effective Hamiltonians to observe the behavior of a wide range of observables, besides the 0νββ) NME, that are compared with experimental data. Following a Bayesian model averaging approach, the range of probable values for the neutrinoless double-beta decay matrix elements is presented. In this paper, we present a similar study for 82Se, which is described in the same model space as 76Ge that is under experimental observation. Due to its faster calculation time compared to 76Ge, 82Se can be used as an appropriate substitute in our complex statistical study. Using the calculations performed for the statistical analysis of the neutrinoless double-beta decay matrix elements, we also search for the correlations between the observables that we can compare to experimental data.

Science & Technology - Other Topics↗

DOE OSTI2014

Trigonometrical sums connected with the chiral Potts model, Verlinde dimension formula, two-dimensional resistor network, and number theory

We have recently developed methods for obtaining exact two-point resistance of the complete graph minus N edges. We use these methods to obtain closed formulas of certain trigonometrical sums that arise in connection with one-dimensional lattice, in proving Scott’s conjecture on permanent of Cauchy matrix, and in the perturbative chiral Potts model. The generalized trigonometrical sums of the chiral Potts model are shown to satisfy recursion formulas that are transparent and direct, and differ from those of Gervois and Mehta. By making a change of variables in these recursion formulas, the dimension of the space of conformal blocks of SU(2) and SO(3) WZW models may be computed recursively. Our methods are then extended to compute the corner-to-corner resistance, and the Kirchhoff index of the first non-trivial two-dimensional resistor network, 2×N. Finally, we obtain new closed formulas for variant of trigonometrical sums, some of which appear in connection with number theory. -- Highlights: • Alternative derivation of certain trigonometrical sums of the chiral Potts model are given. • Generalization of these trigonometrical sums satisfy recursion formulas. • The dimension of the space of conformal blocks may be computed from these recursions. • Exact corner-to-corner resistance, the Kirchhoff index of 2×N are given.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Neutrinoless double beta decay via light neutralinos in R-parity violating supersymmetry

We perform a study of neutrinoless double beta (0νββ) decay mediated by the lightest neutralino of arbitrary mass in the minimal supersymmetric Standard Model (MSSM) under the presence of R-parity violating trilinear interactions. In this scenario, the exchange of the lightest neutralino can result in 0νββ decay of either long-range or short-range behaviour, depending on the neutralino mass. Using nuclear matrix elements calculated in the interacting boson model, we use an interpolation between the long- and short-range behaviours with an approximate formula. The non-observation of 0νββ decay is then used to place constraints on the supersymmetry parameter space, compatible with constraints from collider experiments. We compare these constraints to bounds from pion decays, CKM unitarity and big bang nucleosynthesis.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1991

Matrix cracking in laminated composites under monotonic and cyclic loadings

An analytical model based on the internal state variable (ISV) concept and the strain energy method is proposed for characterizing the monotonic and cyclic response of laminated composites containing matrix cracks. A modified constitution is formulated for angle-ply laminates under general in-plane mechanical loading and constant temperature change. A monotonic matrix cracking criterion is developed for predicting the crack density in cross-ply laminates as a function of the applied laminate axial stress. An initial formulation for a cyclic matrix cracking criterion for cross-ply laminates is also discussed. For the monotonic loading case, a number of experimental data and well-known models are compared with the present study for validating the practical applicability of the ISV approach.

Allen, David H.↗

DOE OSTI2015

Updated constraints on the light-neutrino exchange mechanisms of the 0νββ-decay

The neutrinoless double-beta (0νββ) decay associated with light neutrino exchange mechanisms, which are due to both left-handed V-A and right-handed V+A leptonic and hadronic currents, is discussed by using the recent progress achieved by the GERDA, EXO and KamlandZen experiments. The upper limits for effective neutrino mass m{sub ββ} and the parameters 〈λ〉 and 〈η〉 characterizing the right handed current mechanisms are deduced from the data on the 0νββ-decay of {sup 76}Ge and {sup 136}Xe using nuclear matrix elements calculated within the nuclear shell model and quasiparticle random phase approximation and phase-space factors calculated with exact Dirac wave functions with finite nuclear size and electron screening. The careful analysis of upper constraints on effective lepton number violating parameters assumes a competition of the above mechanisms and arbitrary values of involved CP violating phases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS2006

Equivalent Plate Analysis of Aircraft Wing with Discrete Source Damage

An equivalent plate procedure is developed to provide a computationally efficient means of matching the stiffness and frequencies of flight vehicle wing structures for prescribed loading conditions. First, the equivalent plate is used to match the stiffness of a stiffened panel without damage and the stiffness of a stiffened panel with damage. For both stiffened panels, the equivalent plate models reproduce the deformation of a corresponding detailed model exactly for the given loading conditions. Once the stiffness was matched, the equivalent plate models were then used to predict the frequencies of the panels. Two analytical procedures using the lumped-mass matrix were used to match the first five frequencies of the corresponding detailed model. In both the procedures, the lumped-mass matrix for the equivalent plate is constructed by multiplying the diagonal terms of the consistent-mass matrix by a proportionality constant. In the first procedure, the proportionality constant is selected such that the total mass of the equivalent plate is equal to that of the detailed model. In the second method, the proportionality constant is selected to minimize the sum of the squares of the errors in a set of pre-selected frequencies between the equivalent plate model and the detailed model. The equivalent plate models reproduced the fundamental first frequency accurately in both the methods. It is observed that changing only the mass distribution in the equivalent plate model did not provide enough flexibility to match all of the frequencies.

Krishnamurthy, T.↗

NASA NTRS2000

Constituent-Based Life Models Being Developed for SiC/SiC Composites

For the successful utilization of ceramic matrix composites (CMC) as hot-section components in advanced aeropropulsion engines, the CMC constituents will need to be tailored and optimized to meet all the critical property demands of each component. Under the High-Speed Research (HSR) and Advanced High-Temperature Engine Materials Technology (HITEMP) Programs, the NASA Glenn Research Center at Lewis Field initiated research to develop mechanistic models for key CMC thermostructural properties. These models would describe the effects of different constituent factors (composition, geometry, and volume fraction) and of potential application conditions (stress, time, temperature, and environment) on these properties. Particular focus was placed on both analytical and numerical modeling of state-of-the-art SiC/SiC composites where the primary load-bearing constituents are stoichiometric SiC fibers in a complex multiphase SiC matrix produced by chemical vapor infiltration and melt infiltration. Recent studies have resulted in computer-based numerical models for the elastic modulus, thermal expansion, and thermal conductivity properties of the SiC/SiC system. Additional studies have generated analytical and empirical models for the time dependence of composite rupture strength at temperatures above 2200 F (1200 C), where CMC's have an important thermostructural advantage over current nickel-based superalloys. These life models utilize thermal activation theory and fiber stress-rupture results measured at Glenn to generate Larson-Miller (L M) plots of fiber rupture strength versus q, a single time- and temperature-dependent parameter. Assuming a worse case in which the SiC matrix is cracked, rupture is then controlled by the time-dependent fracture characteristics of the fiber bundles bridging the matrix cracks. With this as the controlling mechanism, one can then use simple composite theory and the fiber L M plots to predict CMC rupture strength versus the q parameter. The dashed line shows the predicted rupture strength of a SiC/SiC composite that is reinforced by a state-of-the-art stoichiometric SiC fiber. For the q parameter, time is in hours and temperature in degrees kelvin. To generate these predictions, a two-dimensional 0 /90 composite with approximately 16 percent fiber in the applied stress direction and an air test environment were assumed. As such, it is possible to compare the model predictions against limited stress-rupture data for this CMC as shown by the data points. The good agreement confirms the rupture model at least for the selected CMC and test conditions. Thus for this particular SiC/SiC composite, one can estimate a 1000-hr rupture strength of approximately 12 ksi at 2400 F (T = 1588 K and q = 39 700 K). At lower CMC application stresses, the SiC matrix is typically uncracked, so both the fiber and matrix constituents share the composite load. In this case, CMC rupture is controlled by the constituent with the longest rupture time based on the creep rate of the composite. Measured Monkman-Grant plots of rupture time versus creep rate for the two SiC constituents have been used to develop CMC life models for this important application condition. NASA and DOD are currently using this information to establish application and material goals for more advanced CMC's that can be used at even higher temperatures.

DiCarlo, James A.↗

DOE OSTI2022

The effect of differential mineral shrinkage on crack formation and network geometry

Rock, concrete, and other engineered materials are often composed of several minerals that change volumetrically in response to variations in the moisture content of the local environment. Such differential shrinkage is caused by varying shrinkage rates between mineral compositions during dehydration. Using both 3D X-ray imaging of geo-architected samples and peridynamic (PD) numerical simulations, we show that the spatial distribution of the clay affects the crack network geometry with distributed clay particles yielding the most complex crack networks and percent damage (99.56%), along with a 60% reduction in material strength. We also demonstrate that crack formation, growth, coalescence, and distribution during dehydration, are controlled by the differential shrinkage rates between a highly shrinkable clay and a homogeneous mortar matrix. Sensitivity tests performed with the PD models show a clay shrinkage parameter of 0.4 yields considerable damage, and reductions in the parameter can result in a significant reduction in fracturing and an increase in material strength. Additionally, isolated clay inclusions induced localized fracturing predominantly due to debonding between the clay and matrix. These insights indicate differential shrinkage is a source of potential failure in natural and engineered barriers used to sequester anthropogenic waste.

58 GEOSCIENCES↗

DOE OSTI2012

Enhanced Generic Phase-field Model of Irradiation Materials: Fission Gas Bubble Growth Kinetics in Polycrystalline UO 2

Experiments show that inter-granular and intra-granular gas bubbles have different growth kinetics which results in heterogeneous gas bubble microstructures in irradiated nuclear fuels. A science-based model predicting the heterogeneous microstructure evolution kinetics is desired, which enables one to study the effect of thermodynamic and kinetic properties of the system on gas bubble microstructure evolution kinetics and morphology, improve the understanding of the formation mechanisms of heterogeneous gas bubble microstructure, and provide the microstructure to macroscale approaches to study their impact on thermo-mechanical properties such as thermo-conductivity, gas release, volume swelling, and cracking. In our previous report 'Mesoscale Benchmark Demonstration, Problem 1: Mesoscale Simulations of Intra-granular Fission Gas Bubbles in UO 2 under Post-irradiation Thermal Annealing', we developed a phase-field model to simulate the intra-granular gas bubble evolution in a single crystal during post-irradiation thermal annealing. In this work, we enhanced the model by incorporating thermodynamic and kinetic properties at grain boundaries, which can be obtained from atomistic simulations, to simulate fission gas bubble growth kinetics in polycrystalline UO 2 fuels. The model takes into account of gas atom and vacancy diffusion, vacancy trapping and emission at defects, gas atom absorption and resolution at gas bubbles, internal pressure in gas bubbles, elastic interaction between defects and gas bubbles, and the difference of thermodynamic and kinetic properties in matrix and grain boundaries. We applied the model to simulate gas atom segregation at grain boundaries and the effect of interfacial energy and gas mobility on gas bubble morphology and growth kinetics in a bi-crystal UO 2 during post-irradiation thermal annealing. The preliminary results demonstrate that the model can produce the equilibrium thermodynamic properties and the morphology of gas bubbles at grain boundaries for given grain boundary properties. More validation of the model capability in polycrystalline is underway.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

NASA NTRS2004

Modeling the Nonlinear, Strain Rate Dependent Deformation of Shuttle Leading Edge Materials with Hydrostatic Stress Effects Included

An analysis method based on a deformation (as opposed to damage) approach has been developed to model the strain rate dependent, nonlinear deformation of woven ceramic matrix composites, such as the Reinforced Carbon Carbon (RCC) material used on the leading edges of the Space Shuttle. In the developed model, the differences in the tension and compression deformation behaviors have also been accounted for. State variable viscoplastic equations originally developed for metals have been modified to analyze the ceramic matrix composites. To account for the tension/compression asymmetry in the material, the effective stress and effective inelastic strain definitions have been modified. The equations have also been modified to account for the fact that in an orthotropic composite the in-plane shear response is independent of the stiffness in the normal directions. The developed equations have been implemented into LS-DYNA through the use of user defined subroutines (UMATs). Several sample qualitative calculations have been conducted, which demonstrate the ability of the model to qualitatively capture the features of the deformation response present in woven ceramic matrix composites.

Goldberg, Robert K.↗

DOE OSTI2017

In silico investigation of blast-induced intracranial fluid cavitation as it potentially leads to traumatic brain injury

In this paper, we conducted computational macroscale simulations predicting blast-induced intracranial fluid cavitation possibly leading to brain injury. To further understanding of this problem, we developed microscale models investigating the effects of blast-induced cavitation bubble collapse within white matter axonal fiber bundles of the brain. We model fiber tracks of myelinated axons whose diameters are statistically representative of white matter. Nodes of Ranvier are modeled as unmyelinated sections of axon. Extracellular matrix envelops the axon fiber bundle, and gray matter is placed adjacent to the bundle. Cavitation bubbles are initially placed assuming an intracranial wave has already produced them. Pressure pulses, of varied strengths, are applied to the upper boundary of the gray matter and propagate through the model, inducing bubble collapse. Simulations, conducted using the shock wave physics code CTH, predict an increase in pressure and von Mises stress in axons downstream of the bubbles after collapse. This appears to be the result of hydrodynamic jetting produced during bubble collapse. Interestingly, results predict axon cores suffer significantly lower shear stresses from proximal bubble collapse than does their myelin sheathing. Finally, simulations also predict damage to myelin sheathing, which, if true, degrades axonal electrical transmissibility and general health of the white matter structures in the brain.

60 APPLIED LIFE SCIENCES↗

NASA NTRS1999

Bone Proteoglycan Changes During Skeletal Unloading

Skeletal adaptability to mechanical loads is well known since the last century. Disuse osteopenia due to the microgravity environment is one of the major concerns for space travelers. Several studies have indicated that a retardation of the mineralization process and a delay in matrix maturation occur during the space flight. Mineralizing fibrillar type I collagen possesses distinct cross-linking chemistries and their dynamic changes during mineralization correlate well with its function as a mineral organizer. Our previous studies suggested that a certain group of matrix proteoglycans in bone play an inhibitory role in the mineralization process through their interaction with collagen. Based on these studies, we hypothesized that the altered mineralization during spaceflight is due in part to changes in matrix components secreted by cells in response to microgravity. In this study, we employed hindlimb elevation (tail suspension) rat model to study the effects of skeletal unloading on matrix proteoglycans in bone.

Yamauchi, M.↗