Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “composite 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 523 records · Page 29

On fracture of fibrous composites

The problem of an existing through-the-thickness center notch in a unidirectional fibrous composite is investigated using a homogeneous anisotropic elasticity analysis in conjunction with various crack growth criteria. Tensor polynomial, strain energy density, and normal stress ratio criteria are applied as crack growth direction criteria and compared to experimental observations. For the cases and criteria studied, only the normal stress ratio criterion exhibits an ability to correctly predict crack growth direction. The normal stress ratio criterion is subsequently applied to the same cases as a failure criterion. Comparison with experiment indicates that these predictions exhibit reasonable correlation with observed critical stresses. Finally, the macroscopic approach to the problem used in this study is compared to a microscopic shear lag analysis of experimental results in the literature. The comparison indicates that both approaches yield essentially identical predictions. The comparisons made with other criteria and approaches indicate that application of the normal stress ratio within a macroscopic level analysis has great potential as an approach to modeling composite crack growth.

Beuth, J. L., Jr.↗

Seismic parameters of hcp‐Fe alloyed with Ni and Si in the Earth's inner core

Abstract Iron alloyed with Ni and Si has been suggested to be a major component of the Earth's inner core. High‐pressure results on the combined alloying effects of Ni and Si on seismic parameters of iron are thus essential for establishing satisfactory geophysical and geochemical models of the region. Here we have investigated the compressional ( V P ) and shear ( Vs ) wave velocity‐density ( ρ ) relations, Poisson's ratio ( ν ), and seismic heterogeneity ratios ( dlnρ / dlnV P , dlnρ / dlnV S , and dlnV P / dlnV S ) of hcp‐Fe and hcp‐Fe 86.8 Ni 8.6 Si 4.6 alloy up to 206 GPa and 136 GPa, respectively, using multiple complementary techniques. Compared with the literature velocity values for hcp‐Fe and Fe‐Ni‐Si alloys, our results show that the combined addition of 9.0 wt % Ni and 2.3 wt % Si slightly increases the V P but significantly decreases the V S of hcp‐Fe at a given density relevant to the inner core. Such distinct alloying effects on velocities of hcp‐Fe produce a high ν of about 0.40 for the alloy at inner core densities, which is approximately 20% higher than that for hcp‐Fe. Analysis of the literature high P‐T results on V P and V S of Fe alloyed with light elements shows that high temperature can further enhance the ν of hcp‐Fe alloyed with Ni and Si. Most significantly, the derived seismic heterogeneity ratios of this hcp alloy present a better match with global seismic observations. Our results provide a multifactored geophysical constraint on the compositional model of the inner core which is consistent with silicon being a major light element alloyed with Fe and 5 wt % Ni.

Liu, Jin↗

Permittivity effects of particle agglomeration in ferroelectric ceramic-epoxy composites using finite element modeling

Here, the size dependence of the dielectric constants of barium titanate or other ferroelectric particles can be explored by embedding particles into an epoxy matrix whose dielectric constant can be measured directly. However, to extract the particle dielectric constant requires a model of the composite medium. We compare a finite element model for various volume fractions and particle arrangements to several effective medium approximations, which do not consider particle arrangement explicitly. For a fixed number of particles, the composite dielectric constant increases with the degree of agglomeration, and we relate this increase to the number of regions of enhanced electric field along the applied field between particles in an agglomerate. Additionally, even for dispersed particles, we find that the composite method of assessing the particle dielectric constant may not be effective if the particle dielectric constant is too high compared to the background medium dielectric constant.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Parallelized Quadrupole Simulations of Thermographic Responses of Composites

Thermography has been shown to be a viable technique for inspection of composites. Model inversion of the thermography data requires a fast method for performing the forward problem. Viable numerical methods for the thermal response forward problem are finite element, finite difference and the quadrupole method. Normally both the finite element and finite difference methods solve for the thermal response in the time domain which limits one’s ability to increase the speed of the simulation by parallelization. In contrast, the quadrupole method solves for the Laplace transform of the thermal response. One of the features of the Laplace transform methodology is the solution at any discrete time is independent of the solution at all other times. Therefore, it is easy to separate into a set of independent calculations with each of the times of interest being performed in parallel. Additionally, the numeric inversion of the Laplace transform typically involves numerically solving for the Laplace transform at multiple Laplace frequencies. Each of those solutions are also independent of solutions at other frequencies and can be calculated in parallel. By parallelization of this method, it is possible to perform the simulations of three-dimensional configurations in seconds. When the input stimulus for thermal response is a delta function heat flux (a reasonable approximation for flash heating), the thermal response is smooth. For this case, it is possible to accurately estimate the thermal response at any time within a given time interval from a set of simulations separated by exponentially increasing time steps. From these simulations, it is possible to accurately interpolate to find the response at intermediate times by a spline interpolation of the logarithm of time versus logarithm of temperature. The thermal response with exponential time stepping is shown to produce values for the thermal response which are within 1% of values within the time interval. The simulations are compared to finite element simulations of the same inspection configurations. The simulations are also compared to the thermographic measurements on composites where shape and depth of the delaminations are obtained from other inspection methods.

Thermography↗

Parallelized Quadrupole Simulations of Thermographic Responses of Composites

Thermography has been shown to be a viable technique for inspection of composites. Model inversion of the thermography data requires a fast method for performing the forward problem. Viable numerical methods for the thermal response forward problem are finite element, finite difference and the quadrupole method. Normally both the finite element and finite difference methods solve for the thermal response in the time domain which limits one’s ability to increase the speed of the simulation by parallelization. In contrast, the quadrupole method solves for the Laplace transform of the thermal response. One of the features of the Laplace transform methodology is the solution at any discrete time is independent of the solution at all other times. Therefore, it is easy to separate into a set of independent calculations with each of the times of interest being performed in parallel. Additionally, the numeric inversion of the Laplace transform typically involves numerically solving for the Laplace transform at multiple Laplace frequencies. Each of those solutions are also independent of solutions at other frequencies and can be calculated in parallel. By parallelization of this method, it is possible to perform the simulations of three-dimensional configurations in seconds. When the input stimulus for thermal response is a delta function heat flux (a reasonable approximation for flash heating), the thermal response is smooth. For this case, it is possible to accurately estimate the thermal response at any time within a given time interval from a set of simulations separated by exponentially increasing time steps. From these simulations, it is possible to accurately interpolate to find the response at intermediate times by a spline interpolation of the logarithm of time versus logarithm of temperature. The thermal response with exponential time stepping is shown to produce values for the thermal response which are within 1% of values within the time interval. The simulations are compared to finite element simulations of the same inspection configurations. The simulations are also compared to the thermographic measurements on composites where shape and depth of the delaminations are obtained from other inspection methods.

Thermography↗

Micromechanical Modeling of High-Strain Thin-Ply Composites

This paper presents a micromechanical analysis for the elastic and viscoelastic behavior of high-strain thin-ply composites. The modeling approach is based on unit cell homogenization. The geometry of the internal weave architecture and ply configuration is characterized via micrographic analysis and explicitly modeled in the unit cell. The composites are modeled as Kirchhoff plates and the homogenization analysis computes the effective relaxation ABD matrix represented by Prony series using the elastic and viscoelastic properties of the constituent fiber and matrix. The formulation of the micromechanical model and numerical implementation d. Composite laminates with 3-ply and 4-ply configurations are studied.

Hamillage, Milinda Yapa↗

Thermal equation of state of hcp-iron: Constraint on the density deficit of Earth's solid inner core: THERMAL EQUATION OF STATE OF HCP-IRON

We conducted high-pressure experiments on hexagonal close packed iron (hcp-Fe) in MgO, NaCl, and Ne pressure-transmitting media and found general agreement among the experimental data at 300 K that yield the best fitted values of the bulk modulus K 0 = 172.7(±1.4) GPa and its pressure derivative K 0 '= 4.79(±0.05) for hcp-Fe, using the third-order Birch-Murnaghan equation of state. Using the derived thermal pressures for hcp-Fe up to 100 GPa and 1800 K and previous shockwave Hugoniot data, we developed a thermal equation of state of hcp-Fe. The thermal equation of state of hcp-Fe is further used to calculate the densities of iron along adiabatic geotherms to define the density deficit of the inner core, which serves as the basis for developing quantitative composition models of the Earth's inner core. We determine the density deficit at the inner core boundary to be 3.6%, assuming an inner core boundary temperature of 6000 K.

58 GEOSCIENCES↗

Reactor fuel fraction information on the antineutrino anomaly

We analyzed the evolution data of the Daya Bay reactor neutrino experiment in terms of short-baseline active-sterile neutrino oscillations taking into account the theoretical uncertainties of the reactor antineutrino fluxes. We found that oscillations are disfavored at 2.6σ with respect to a suppression of the 235 U reactor antineutrino flux and at 2.5σ with respect to variations of the 235 U and 239 Pu fluxes. On the other hand, the analysis of the rates of the short-baseline reactor neutrino experiments favor active-sterile neutrino oscillations and disfavor the suppression of the 235 U flux at 3.1σ and variations of the 235 U and 239 Pu fluxes at 2.8σ. We also found that both the Daya Bay evolution data and the global rate data are well-fitted with composite hypotheses including variations of the 235 U or 239 Pu fluxes in addition to active-sterile neutrino oscillations. A combined analysis of the Daya Bay evolution data and the global rate data shows a slight preference for oscillations with respect to variations of the 235 U and 239 Pu fluxes. However, the best fits of the combined data are given by the composite models, with a preference for the model with an enhancement of the 239 Pu flux and relatively large oscillations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Learning from Higgs physics at future Higgs factories

Future Higgs factories can reach impressive precision on Higgs property measurements. In this paper, instead of conventional focus of Higgs precision in certain interaction bases, we explored its sensitivity to new physics models at the electron-positron colliders. In particular, we studied two categories of new physics models, Standard Model (SM) with a real scalar singlet extension, and Two Higgs Double Model (2HDM) as examples of weakly-interacting models, Minimal Composite Higgs Model (MCHM) and three typical patterns of the more general operator counting for strong interacting models as examples of strong dynamics. We performed a global fit to various Higgs search channels to obtain the 95% C.L. constraints on the model parameter space. In the SM with a singlet extension, we obtained the limits on the singlet-doublet mixing angle sin(theta), as well as the more general Wilson coefficients of the induced higher dimensional operators. In the 2HDM, we analyzed tree level effects in tan(beta) vs. cos(beta-alpha) plane, as well as the one-loop contributions from the heavy Higgs bosons in the alignment limit to obtain the constraints on heavy Higgs masses for different types of 2HDM. In strong dynamics models, we obtained lower limits on the strong dynamics scale. In addition, once deviations of Higgs couplings are observed, they can be used to distinguish different models. Here, we also compared the sensitivity of various future Higgs factories, namely Circular Electron Positron Collider (CEPC), Future Circular Collider (FCC)-ee and International Linear Collider (ILC).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Drivers and potential predictability of summer time North Atlantic polar front jet variability

The variability of the North Atlantic polar front jet stream is crucial in determining summer weather around the North Atlantic basin. Recent extreme summers in western Europe and North America have highlighted the need for greater understanding of this variability, in order to aid seasonal forecasting and mitigate societal, environmental and economic impacts. Here we find that simple linear regression and composite models based on a few predictable factors are able to explain up to 35 % of summertime jet stream speed and latitude variability from 1955 onwards. Sea surface temperature forcings impact predominantly on jet speed, whereas solar and cryospheric forcings appear to influence jet latitude. The cryospheric associations come from the previous autumn, suggesting the survival of an ice-induced signal through the winter season, whereas solar influences lead jet variability by a few years. Regression models covering the earlier part of the twentieth century are much less effective, presumably due to decreased availability of data, and increased uncertainty in observational reanalyses. Wavelet coherence analysis identifies that associations fluctuate over the study period but it is not clear whether this is just internal variability or genuine non-stationarity. Finally we identify areas for future research.

Meteorology & Atmospheric Sciences↗

Porosity Prediction in A356 Wheel Casting

A FLUENT™-based thermal-fluid-compositional model has been developed and applied to a low-pressure, die-cast, A356 aluminum alloy wheel to explore the effects of macro-segregation of hydrogen and late-stage liquid encapsulation on pore formation during solidification. Based on the evolution of temperature, pressure, and hydrogen concentration in the liquid output from FLUENT, the porosity size distribution has been predicted at a number of locations in a commercially produced wheel using an in-house micro-porosity model. The porosity model incorporates pore nucleation kinetics based on hydrogen supersaturation, and two modes of pore growth: the first, occurring at low solid fractions is based on hydrogen diffusion; and the second, occurring at high solid fractions is based on conservation of volume. Samples from a commercially cast wheel have been analyzed using X-ray micro-tomography to provide basic validation of the micro-porosity model. The predicted results showed that the evolution in pressure has the dominant effect on pore growth, but only under conditions where pores have nucleated prior to late-stage liquid encapsulation. Otherwise, the cooling rate appears to have the dominant effect. The accuracy of the model is discussed in terms of the predicted pore size distribution and number density of pores. Areas for continued development are presented.

36 MATERIALS SCIENCE↗

Synchrotron X-ray microbeam diffraction measurements of full elastic long range internal strain and stress tensors in commercial-purity aluminum processed by multiple passes of equal-channel angular pressing

Synchrotron X-ray microbeam diffraction was used to measure the full elastic long range internal strain and stress tensors of low dislocation density regions within the submicrometer grain/subgrain structure of equal-channel angular pressed (ECAP) aluminum alloy AA1050 after 1, 2, and 8 passes using route B C . This is the first time that full tensors were measured in plastically deformed metals at this length scale. The maximum (most tensile or least compressive) principal elastic strain directions for the unloaded 1 pass sample for the grain/subgrain interiors align well with the pressing direction, and are more random for the 2 and 8 pass samples. The measurements reported here indicate that the local stresses and strains become increasingly isotropic (homogenized) with increasing ECAP passes using route BC. The average maximum (in magnitude) LRISs are -0.43 σ a for 1 pass, -0.44 σ a for 2 pass, and 0.14 σ a for the 8 pass sample. Furthermore, these LRISs are larger than those reported previously because those earlier measurements were unable to measure the full stress tensor. Significantly, the measured stresses are inconsistent with the two-component composite model.

36 MATERIALS SCIENCE↗

Effect of ion interactions on the Raman spectrum of NO 3 − : Toward monitoring of low-activity nuclear waste at Hanford

Raman spectroscopy is a valuable in-situ technique for many applications. The concentrations of species in complex ionic mixtures of nuclear waste can be estimated using Raman measurements. However, it has been experimentally observed that ion interactions can cause a modification of Raman peak intensities and positions. The present work explores the nonlinear behavior associated with the nitrate anion, NO 3 − , which is present in abundance in low activity nuclear waste. We examine changes in the main Raman peak of the nitrate anion in the presence of other ions. A wide range of concentrations are covered, including those expected during direct-feed low-activity waste (DFLAW) processing at the Hanford site in the State of Washington. The experiments showed that the ions interact and associate to form ion pairs, which results in a blue shift (i.e., a shift towards higher wavenumbers) in the main Raman peak of nitrate. The results indicate that cation concentration is a better predictor of the peak shift, compared to ionic strength, both for binary and multicomponent mixtures. These findings have direct implications on the development of spectra-to-composition models for the DFLAW system, since they show deviations from the linearity assumptions used in common chemometric models.

42 ENGINEERING↗

Corrosion resistant hot melt adhesive to bind metals

Hot melt adhesives (HMAs) play an important role in many industries, and their demand is expected to grow. HMAs don't require any solvents, and their application results in the formation of strong bonds with the substrate upon cooling within seconds. These properties differentiate them from liquid glues and make them preferable for practical application. Currently, commercial HMAs are used in bonding lightweight materials such as paper, polymers, and cartons, and have limited usage in areas necessitating the bonding of heavier objects like metals. Here, in this study, we report a design and testing of versatile platform comprising an ion-coordinating polymer and ionic fillers for performance optimization and understanding of structure–property relationships, enabling the rational design of HMAs with improved adhesion to metal surfaces. All-atom Molecular dynamics (MD) simulations and various characterization methods are used to elucidate the adhesion mechanism in model composite system containing polyethylene oxide mixed with chemically diverse salt particles. The maximum adhesion strength is found in composites with Al(OH) 3 and FeCl 3 fillers. Interestingly, the presence of Al(OH) 3 also provides a multifunctional anticorrosion property as measured electrochemically using the Tafel method. The discovered path to formulations with improved adhesion to metal surfaces constitutes an important step toward advancing HMAs for use in the structural and semi-structural metal work domain.

Adhesive and cohesive forces↗

Stability of Fe-bearing hydrous phases and element partitioning in the system MgO–Al 2 O 3 –Fe 2 O 3 –SiO 2 –H 2 O in Earth’s lowermost mantle

We performed high pressure–temperature ( P-T ) experiments on a model composition of hydrous subducted slabs in the MgO–Al 2 O 3 –Fe 2 O 3 –SiO 2 –H 2 O system using laser-heated diamond anvil cells. The phase assemblages were characterized combining in-situ synchrotron X-ray diffraction and ex-situ transmission electron microscope techniques. The hydrous δ-phase AlOOH–FeOOH–MgSiO 2 (OH) 2 –SiO 2 coexists with bridgmanite (Bdg), post-perovskite (pPv), or both in a broad P-T range of 104–126 GPa and 1750–2500 K. The hydrous pyrite-type FeOOHx phase was observed in-situ over a P-T range of 112–123 GPa and 1750–2300 K, coexisting with the pPv phase. Chemical analysis on recovered samples showed that considerable amount of Fe 2 O 2 (OH) 2 (8–13 mol%) and SiO 2 (9–13 mol%) in the δ-phase does not reduce its thermal stability compared to the Al-endmember, indicating that the Fe-bearing δ-phase can transport water to the lowermost mantle along the mantle geotherm. In this hydrous system, we observed that Al depletion in both the Bdg and pPv phases can significantly reduce the width of the Bdg to pPv transition in contrast to a wide two-phase coexistence region in a dry Al-rich system. Meanwhile, the Fe enrichment in the pPv phase relative to the coexisting Bdg phase lowers the transition pressure to the depth of the D" discontinuity. Accordingly, the depth and thickness of the Bdg to pPv transition in subducted basaltic crustal materials can explain the seismically detected D" discontinuity. Partial melting could be triggered by dehydration of the Fe-bearing hydrous phases due to a steep temperature gradient at the core mantle boundary (CMB), and therefore the ultralow-velocity zones (ULVZs) might be the regions where partial melting occurs at the lowermost mantle.

58 GEOSCIENCES↗

Sound velocities and thermal equation of state of fcc -iron-nickel alloys at high pressure and high temperature: Implications for the cores of Moon and several planets

Fcc-Fe-Ni alloy is believed to be the most dominant solid constitute of moderate-sized terrestrial planetary cores. Investigating the physical properties, especially the density and sound velocity of Fe-Ni alloys and comparing them with seismic observations is an indispensable approach to constructing compositional models for planetary interiors. In this study, we conducted sound velocity measurements on Fe-Ni alloys with 10 wt.% and 20 wt.% Ni up to ∼13.5 GPa and 1073 K, using the ultrasonic interferometry technique in a multi-anvil apparatus in conjunction with synchrotron radiation. By fitting the experimental data to finite strain equations, the bulk and shear moduli and their pressure and temperature derivatives are derived, yielding K S0 =145.8(14) GPa, G 0 = 73.2(7) GPa, K S0 ’ = 5.89(24), G 0 ’ = 2.89(8), (∂K S /∂T) P = -0.0181(12) GPa/K and (∂G/∂T) P = -0.0393(10) GPa/K for fcc-Fe 80 Ni 20 . An examination of the density-velocity relationship shows that compressional wave velocity is insensitive to temperature within the current pressure and temperature range, while shear wave velocity exhibits a large reduction with increasing temperature. Here, extrapolation of the sound velocities following the finite strain theories suggests that much slower Vs should be expected at pressure and temperature conditions corresponding to those of the lunar core. Possible core density and velocity profiles for other moderate planets and satellites, such as Mars, Mercury, and Ganymede are also calculated.

Equation of state↗

Application of pH-dependent environmental protection agency test method 1313 to low-activity nuclear waste glass

Here, in this study, we present results of a room temperature, pH-dependent leach test established by United States Environmental Protection Agency (EPA) Leaching Environmental Assessment Framework (LEAF) (EPA Method 1313 and further abbreviated EPA 1313) to study the alteration behavior of low-activity nuclear waste (LAW) glasses. Implementation of EPA 1313, which is closer to LAW disposal temperatures, could allow for more flexibility in the formulation of Hanford LAW glasses, specifically in terms of increased waste loading. A set of 48 LAW glasses were selected to test with EPA 1313. Briefly, the test method consists of adding glass to a number of solutions at different initial pH values and at a fixed mass-to-solution-volume ratio. The glass is then allowed to corrode for two days and the final solution pH and concentration of different glass elements in solution is measured. The test was modified from the published procedure to constrain particle size to provide more reproducible results, reduce the studied pH range to a more applicable range, and to reduce solution volume to reduce costs associated with material preparation. Compositional modeling was used to predict the measured elemental releases at various pH. The model predictions improved in tests where there was a relatively small difference between initial and final solution pH indicating that pH controls would result in better predictive models.

Accelerated aging↗