Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “elastic fields”

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 973 records · Page 54

Progress of the CHILLAX Experiment (HEPCAT Year 1 Report)

A significant facet of modern experimental particle physics research is that of searching for hypothetical rare events. In this context, “rare” means “unusual to detect in the laboratory”, and can be a confirmed phenomenon (e.g. coherent elastic neutrino-nucleus scattering, which has a high cross section at low, difficult to detect energies), or an as-of-yet unobserved phenomenon that may not exist (e.g. dark matter - standard model matter non-gravitational interaction). A salient challenge for experimentalists is to devise detectors optimized for observing such rare events.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Imposition of local boundary conditions in peridynamics without a fictitious layer and unphysical stress concentrations

Here, this study introduces a general approach for the imposition of local boundary conditions in non- ordinary state-based peridynamics (NOSB PD) to eliminate the displacement kinks near the boundary without a fictitious layer under quasi-static loading conditions. It identifies the underlying reason for the unphysical displacement kinks. Under an imposed linear displacement field, the NOSB PD equilibrium equation is not satisfied near the boundary due to the unsymmetric horizon of material points. However, the equilibrium equation derived by using the PD differential operator is satisfied at such material points. Therefore, the material body is divided into three regions to satisfy the equilibrium equations and to impose displacement and tractions boundary conditions. This approach does not deviate from the original NOSB PD; however, it provides a simple solution to eliminate the displacement kink near the boundary, which leads to unphysical stress concentrations. Its efficacy is demonstrated by considering elastic rectangular and square plates subjected to various types of boundary conditions leading to homogeneous as well as nonhomogeneous deformations. The creep response of a rectangular plate further proves the robustness of the present approach. Also, a quasi-static crack propagation from a pre-existing crack in a square plate under mode-I, mode-II and mixed-mode loading conditions demonstrates its capability for failure prediction based on the critical stretch criteria. Finally, its applicability for 3D analysis is demonstrated by considering a rectangular prism under applied stretch and normal stress.

42 ENGINEERING↗

Active doping controls the mode of failure in dense colloidal gels

Mechanical properties of disordered materials are governed by their underlying free energy landscape. In contrast to external fields, embedding a small fraction of active particles within a disordered material generates nonequilibrium internal fields, which can help to circumvent kinetic barriers and modulate the free energy landscape. In this work, we investigate through computer simulations how the activity of active particles alters the mechanical response of deeply annealed polydisperse colloidal gels. We show that the “swim force” generated by the embedded active particles is responsible for determining the mode of mechanical failure, i.e., brittle vs. ductile. We find, and theoretically justify, that at a critical swim force the mechanical properties of the gel decrease abruptly, signaling a change in the mode of mechanical failure. The weakening of the elastic modulus above the critical swim force results from the change in gel porosity and distribution of attractive forces among gel particles, while below the critical swim force, the ductility enhancement is caused by an increase of gel structural disorder. Above the critical swim force, the gel develops a pronounced heterogeneous structure characterized by multiple pore spaces, and the mechanical response is controlled by dynamical heterogeneities. We contrast these results with those of a simulated monodisperse gel that exhibits a nonmonotonic trend of ductility modulation with increasing swim force, revealing a complex interplay between the gel energy landscape and embedded activity.

Zhou, Tingtao (ORCID:000000021766719X)↗

Theoretical Analysis of Fractional Viscoelastic Flow in Circular Pipes: General Solutions

Fractional calculus is a relatively old yet emerging field of mathematics with the widest range of engineering and biomedical applications. Despite being an incredibly powerful tool, it, however, requires promotion in the engineering community. Rheology is undoubtedly one of the fields where fractional calculus has become an integral part of cutting-edge research. There exists extensive literature on the theoretical, experimental, and numerical treatment of various fractional viscoelastic flows in constraint geometries. However, the general theoretical approach that unites several most commonly used models is missing. Here we present exact analytical solutions for fractional viscoelastic flow in a circular pipe. We find velocity profiles and shear stresses for fractional Maxwell, Kelvin–Voigt, Zener, Poynting–Thomson, and Burgers models. The dynamics of these quantities are studied with respect to normalized pipe radius, fractional orders, and elastic moduli ratio. Three different types of behavior are identified: monotonic increase, resonant, and aperiodic oscillations. The models developed are applicable in the widest material range and allow for the alteration of the balance between viscous and elastic properties of the materials.

Gritsenko, Dmitry (ORCID:0000000248181512)↗

How well can we Measure the Vertical Profile of Tropospheric Aerosol Extinction?

The recent Department of Energy Atmospheric Radiation Measurement (ARM) Aerosol Intensive Operations Period (MOP, May 2003) yielded one of the best measurement sets obtained to-date to assess our ability to measure the vertical profile of ambient aerosol extinction sigma(sub ep)(lambda) in the lower troposphere. During one month, a heavily instrumented aircraft with well characterized aerosol sampling ability carrying well proven and new aerosol instrumentation, devoted most of the 60 available flight hours to flying vertical profiles over the heavily instrumented ARM Southern Great Plains (SGP) Climate Research Facility (CRF). This allowed us to compare vertical extinction profiles obtained from 6 different instuments: airborne Sun photometer (AATS-14), airborne nephelometer/absorption photometer, airborne cavity ring-down system, ground-based Raman lidar and 2 ground-based elastic backscatter lidars. We find the in-situ measured sigma(sub ep)(lambda) to be lower than the AATS-14 derived values. Bias differences are 0.002 - 0.004 K/m equivalent to 12-17% in the visible, or 45% in the near-infrared. On the other hand, we find that with respect to AATS-14, the lidar sigma(sub ep)(lambda) are higher. An unnoticed loss of sensitivity of the Raman lidar had occurred leading up to AIOP and we expect better agreement from the recently restored system looking at the collective results from 6 field campaigns conducted since 1996, airborne in situ measurements of sigma(sub ep)(lambda) tend to be biased slightly low (17% at visible wavelengths) when compared to airborne Sun photometer sigma(sub ep)(lambda). On the other hand, sigma(sub ep)(lambda) values derived from lidars tend to have no or positive biases. From the bias differences we conclude that the typical systematic error associated with measuring the tropospheric vertical profile of the ambient aerosol extinction with current state of-the art instrumentation is 15-20% at visible wavelengths and potentially larger in the UV and near-infrared.

Schmid, Beat↗

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.↗

Elastic-Plastic Nonlinear Response of a Space Shuttle External Tank Stringer: Stringer-Feet Imperfections and Assembly - Part 1

Elastic-plastic, large-deflection nonlinear stress analyses are performed for the external hat-shaped stringers (or stiffeners) on the intertank portion of the Space Shuttle s external tank. These stringers are subjected to assembly strains when the stringers are initially installed on an intertank panel. Four different stringer-feet configurations including the baseline flat-feet, the heels-up, the diving-board, and the toes-up configurations are considered. The assembly procedure is analytically simulated for each of these stringer configurations. The location, size, and amplitude of the strain field associated with the stringer assembly are sensitive to the assumed geometry and assembly procedure. The von Mises stress distributions from these simulations indicate that localized plasticity will develop around the first eight fasteners for each stringer-feet configuration examined. However, only the toes-up configuration resulted in high assembly hoop strains.

Knight, Norman F., Jr.↗

Level-set topology optimization with PDE generated conformal meshes

This paper presents a level-set topology optimization approach that uses conformal meshes for the analysis of the displacement field. The structure’s boundary is represented by the iso-contour of a level-set field discretized on a fixed background design mesh. The conformal mesh is updated for each design iteration via a PDE based mesh morphing process that identifies the set of facets in the background mesh that are homeomorphic to the boundary and relaxes the homeomorphic mesh to conform to the structure’s boundary and ensure high element quality. The conformal mesh allows for a more accurate computation of the response versus density and some level-set based methods which interpolate material properties using the volume fraction. Numerical examples illustrate the proposed approach by optimizing linear-elastic two- and three-dimensional structures, wherein insight into the performance of the mesh morphing process is provided. The examples also highlight the scalability of the approach.

42 ENGINEERING↗

The photon content of the proton in the CT18 global analysis

Recently, two photon PDF sets based on implementations of the LUX ansatz into the CT18 global analysis were released. In CT18lux, the photon PDF is calculated directly using the LUX master formula for all scales, \mu μ . In an alternative realization, CT18qed, the photon PDF is initialized at the starting scale, \mu_0 μ 0 , using the LUX formulation and evolved to higher scales \mu(>\mu_0) μ ( > μ 0 ) with a combined QED+QCD kernel at \mathcal{O}(\alpha),~\mathcal{O}(\alpha\alpha_s) 𝒪 ( α ) , 𝒪 ( α α s ) and \mathcal{O}(\alpha^2) 𝒪 ( α 2 ) . In the small- x x region, the photon PDF uncertainty is mainly induced by the quark and gluon PDFs, through the perturbative DIS structure functions. In comparison, the large- x x photon uncertainty comes from various low-energy, nonperturbative contributions, including variations of the inelastic structure functions in the resonance and continuum regions, higher-twist and target-mass corrections, and elastic electromagnetic form factors of the proton. We take the production of doubly-charged Higgs pairs, (H^{++}H^{--}) ( H + + H − − ) , as an example of scenarios beyond the Standard Model to illustrate the phenomenological implications of these photon PDFs at the LHC.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

An equivalent domain integral for analysis of two-dimensional mixed mode problems

An equivalent domain integral (EDI) method for calculating J-integrals for two-dimensional cracked elastic bodies subjected to mixed mode loading is presented. The total and product integrals consist of the sum of an area or domain integral and line integrals on the crack faces. The EDI method gave accurate values of the J-integrals for two mode I and two mixed mode problems. Numerical studies showed that domains consisting of one layer of elements are sufficient to obtain accurate J-integral values. Two procedures for separating the individual modes from the domain integrals are presented. The procedure that uses the symmetric and antisymmetric components of the stress and displacement fields to calculate the individual modes gave accurate values of the integrals for all the problems analyzed.

Raju, I. S.↗

A numerical model for predicting crack path and modes of damage in unidirectional metal matrix composites

A finite element-based numerical technique has been developed to simulate damage growth in unidirectional composites. This technique incorporates elastic-plastic analysis, micromechanics analysis, failure criteria, and a node splitting and node force relaxation algorithm to create crack surfaces. Any combination of fiber and matrix properties can be used. One of the salient features of this technique is that damage growth can be simulated without pre-specifying a crack path. In addition, multiple damage mechanisms in the forms of matrix cracking, fiber breakage, fiber-matrix debonding and plastic deformation are capable of occurring simultaneously. The prevailing failure mechanism and the damage (crack) growth direction are dictated by the instantaneous near-tip stress and strain fields. Once the failure mechanism and crack direction are determined, the crack is advanced via the node splitting and node force relaxation algorithm. Simulations of the damage growth process in center-slit boron/aluminum and silicon carbide/titanium unidirectional specimens were performed. The simulation results agreed quite well with the experimental observations.

Bakuckas, J. G.↗

Temporal and Spatial Evolution of Non-Elastic Strain Accumulation in Stanstead Granite During Brittle Creep

Understanding the long-term behavior of brittle rocks requires fundamental consideration of time-dependent strain evolution and brittle creep processes. Previous studies have evaluated sub-critical crack growth during time-dependent deformation and damage evolution in brittle rocks; however, there is an incomplete knowledge of how damage evolves spatially and temporally within the body of intact rocks, where distributed regions of damage interact and coalesce during creep. This paper presents laboratory research focusing on evaluating brittle creep damage processes in Stanstead granite (SG) using 2-dimensional digital image correlation (2D-DIC). In the laboratory, the prismatic SG specimens were loaded beyond an estimated Crack Damage stress threshold (CD) level and then maintained a constant stress to initiate the creep process. DIC was used to characterize full-field spatiotemporal strain evolution, which was then interpreted in the context of local regions of “damage”, determined according to a strain-based criterion. Here, a method was proposed for identifying “existing” and “new” damage regions over specified intervals during the test, followed by spatial clustering of these regions to assess their spatiotemporal evolution. The clustering analysis results demonstrated the extension of existing damage regions was the main damage process during brittle creep, which is consistent with existing models of sub-critical crack growth. In addition, temporal analysis of tensile and shear strains on a point-by-point basis revealed both new damage formation and the strain concentration within existing damaged regions significantly contribute to overall specimen strain during primary creep. In contrast, during secondary creep, increases in specimen deformation are influenced by the accumulation of strains within already damaged regions.

58 GEOSCIENCES↗

Proximate ferromagnetic state in the Kitaev model material α-RuCl 3

α-RuCl 3 is a major candidate for the realization of the Kitaev quantum spin liquid, but its zigzag antiferromagnetic order at low temperatures indicates deviations from the Kitaev model. We have quantified the spin Hamiltonian of α-RuCl 3 by a resonant inelastic x-ray scattering study at the Ru L 3 absorption edge. In the paramagnetic state, the quasi-elastic intensity of magnetic excitations has a broad maximum around the zone center without any local maxima at the zigzag magnetic Bragg wavevectors. This finding implies that the zigzag order is fragile and readily destabilized by competing ferromagnetic correlations. The classical ground state of the experimentally determined Hamiltonian is actually ferromagnetic. The zigzag state is stabilized by quantum fluctuations, leaving ferromagnetism – along with the Kitaev spin liquid – as energetically proximate metastable states. The three closely competing states and their collective excitations hold the key to the theoretical understanding of the unusual properties of α-RuCl 3 in magnetic fields.

36 MATERIALS SCIENCE↗

Nucleon axial-vector form factor and radius from future neutrino experiments

Precision measurements of antineutrino elastic scattering on hydrogen from future neutrino experiments offer a unique opportunity to access the low-energy structure of protons and neutrons. We discuss the determination of the nucleon axial-vector form factor and radius from antineutrino interactions on hydrogen that can be collected at the future Long-Baseline Neutrino Facility and study the sources of theoretical and experimental uncertainties. The projected accuracy would improve existing measurements by 1 order of magnitude and be competitive with contemporary lattice-QCD determinations, potentially helping to resolve the corresponding tension with measurements from (anti)neutrino elastic scattering on deuterium. We find that the current knowledge of the nucleon vector form factors could be one of the dominant sources of uncertainty. We also evaluate the constraints that can be simultaneously obtained on the absolute ν ¯ μ flux normalization. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Dynamic response of 17-4 stainless steel as a function of manufacturing method and heat treatment

We present a series of plate-impact experiments on 17-4 stainless steel to study the effect of manufacturing method and heat treatment on the Hugoniot elastic limit (HEL), Hugoniot, phase transformation stress, and spallation strength. Two traditional manufacturing methods were considered, wrought processing and casting, as well as two additive manufacturing methods, laser powder-bed fusion (LPBF) and wire-fed electron beam (EBAM). For both LPBF and EBAM 17-4 stainless steel variants, two billets were printed, enabling the application of two unique heat treatments. The HEL stress depended heavily on the thermal history, with the HEL increasing after the formation of Cu-rich precipitates via heat treatment. The Hugoniot response both below and above the phase transition was unaffected by the manufacturing method or heat treatment. The phase transition stress depended heavily on the thermal history, with its variation being attributed to the presence of various microstructural features. This is supported by a marked increase in the phase transition stress after precipitation hardening. These results suggest that the notion of the phase transition stress being dictated by bulk composition is an oversimplification and the stress fields generated by the meso-scale structure are a dominant force. The spallation strength was lower in the cast material compared to all other 17-4 stainless steel variants due to the presence of brittle δ-ferrite inclusions. Additionally, a drop in the tensile strain-rate was observed in the spallation response above the phase transition stress, which was hypothesized to stem from the kinetics of the reversion to the low-pressure phase during spall.

Compressive stress↗

Electron‐Assisted Generation and Straight Movement of Skyrmion Bubble in Kagome TbMn 6 Sn 6

Abstract Topological magnetic textures are promising candidates as binary data units for the next‐generation memory device. The precise generation and convenient control of nontrivial spin topology at zero field near room temperature endows the critical advantages in skyrmionic devices but is not simultaneously integrated into one material. Here, in the Kagome plane of quantum TbMn 6 Sn 6 , the expedient generation of the skyrmion bubbles in versatile forms of lattice, chain, and isolated one by converging the electron beam, where the electron intensity gradient contributes to the dynamic generation from local anisotropy variation near spin reorientation transition (SRT) is reported. Encouragingly, by utilizing the dynamic shift of the SRT domain interface, the straight movement is actualized with the skyrmion bubble slave to the SRT domain interface forming an elastic composite object, avoiding the usual deflection from the skyrmion Hall effect. The critical contribution of the SRT domain interface via conveniently electron‐assisted heating is further theoretically validated in micromagnetic simulation, highlighting the compatible application possibility in advanced devices.

Li, Zhuolin↗

Crystal structure of calcium-ferrite type NaAlSiO4 up to 45 GPa

Abstract Alkali-rich aluminous high-pressure phases including calcium-ferrite (CF) type NaAlSiO4 are thought to constitute ~20% by volume of subducted mid-ocean ridge basalt (MORB) under lower mantle conditions. As a potentially significant host for incompatible elements in the deep mantle, knowledge of the crystal structure and physical properties of CF-type phases is therefore important to understanding the crystal chemistry of alkali storage and recycling in the Earth’s mantle. We determined the evolution of the crystal structure of pure CF-NaAlSiO4 and Fe-bearing CF-NaAlSiO4 at pressures up to ~45 GPa using synchrotron-based, single-crystal X-ray diffraction. Using the high-pressure lattice parameters, we also determined a third-order Birch-Murnaghan equation of state, with V0 = 241.6(1) Å3, KT0 = 220(4) GPa, and KT0′ = 2.6(3) for Fe-free CF, and V0 = 244.2(2) Å3, KT0 = 211(6) GPa, and KT0′ = 2.6(3) for Fe-bearing CF. The addition of Fe into CF-NaAlSiO4 resulted in a 10 ± 5% decrease in the stiffest direction of linear compressibility along the c-axis, leading to stronger elastic anisotropy compared with the Fe-free CF phase. The NaO8 polyhedra volume is 2.6 times larger and about 60% more compressible than the octahedral (Al,Si)O6 sites, with K0NaO8 = 127 GPa and K0(Al,Si)O6 ~304 GPa. Raman spectra of the pure CF-type NaAlSiO4 sample shows that the pressure coefficient of the mean vibrational mode, 1.60(7) cm–1/GPa, is slightly higher than 1.36(6) cm−1/GPa obtained for the Fe-bearing CF-NaAlSiO4 sample. The ability of CF-type phases to contain incompatible elements such as Na beyond the stability field of jadeite requires larger and less-compressible NaO8 polyhedra. Detailed high-pressure crystallographic information for the CF phases provides knowledge on how large alkali metals are hosted in alumina framework structures with stability well into the lowermost mantle.

Geochemistry & Geophysics↗