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 523 records · Page 29

The excitation of electronic transverse energy levels in an intense magnetic field

Observations of the X-ray pulsar Hercules X-1 show a line emission feature at about 60 keV, which has been interpreted as the fundamental electron cyclotron line in a magnetic field of around six trillion gauss. In this interpretation, the line radiation results from transitions between transverse energy levels, which are quantized by the field. The expected line luminosity from the excitation of these levels by protons which are falling into the polar cap of a neutron star are calculated. They are assumed to attain kinetic energies up to around 200 MeV, the gravitational potential energy at the surface. The cross sections for high energy Coulomb encounters between small pitch angle protons and electrons in a strong field are measured and used to calculate the energy loss rate of the infalling protons. This rate, together with the rate of elastic nuclear proton collisions, is then used to calculate the number of line photons an infalling proton can be expected to produce, directly or indirectly. The results are applied to Hercules X-1.

Bussard, R. W.↗

Internal stresses in Phobos and other triaxial bodies

The unusual dynamical behavior of Phobos, its strange appearance, and its mysterious network of grooves all make it an intriguing object. Geophysical studies, though, have been hampered by the lack of suitable theories applicable to nonspherical bodies. In this paper the Martian satellites are modeled as homogeneous, elastic triaxial ellipsoids subject to tidal, rotational, and self-gravitational stresses. A novel semianalytical treatment then gives the stress and strain fields throughout their interiors. Yield phenomena and their possible surface expressions are also investigated. The results indicate that Phobos and Deimos have always been stable with respect to tidal fracture or disruption, but that Phobos will probably break up before colliding with Mars. Applications of the new formulation to other nonspherical bodies in the solar system are also discussed.

Dobrovolskis, A. R.↗

Magnetospheric interaction with Triton's ionosphere

The large electron densities measured by the Voyager radio occultation experiment are attributed to the precipitation of magnetospheric electrons with energy above 10 keV. Because the ionospheric electric Pedersen conductivity of Triton is about 10,000-20,000 mho and the Alfven conductance is about 3.5 mho, direct convective flow of plasma into the essentially infinitely conducting ionosphere is negligible. Magnetospheric electrons are transported to Triton's ionopause by curvature drift as a result of weak magnetic field line draping in a sub-Alfvenic plasma interaction with Triton. At the ionopause energetic electrons have a high probability of elastic and inelastic scattering and precipitate into the upper atmosphere. The average power dissipation is estimated to be about (2 - 3) x 10 to the 8th W.

Strobel, Darrell F.↗

Flux Compression Magnetic Nozzle

In pulsed fusion propulsion schemes in which the fusion energy creates a radially expanding plasma, a magnetic nozzle is required to redirect the radially diverging flow of the expanding fusion plasma into a rearward axial flow, thereby producing a forward axial impulse to the vehicle. In a highly electrically conducting plasma, the presence of a magnetic field B in the plasma creates a pressure B(exp 2)/2(mu) in the plasma, the magnetic pressure. A gradient in the magnetic pressure can be used to decelerate the plasma traveling in the direction of increasing magnetic field, or to accelerate a plasma from rest in the direction of decreasing magnetic pressure. In principle, ignoring dissipative processes, it is possible to design magnetic configurations to produce an 'elastic' deflection of a plasma beam. In particular, it is conceivable that, by an appropriate arrangement of a set of coils, a good approximation to a parabolic 'magnetic mirror' may be formed, such that a beam of charged particles emanating from the focal point of the parabolic mirror would be reflected by the mirror to travel axially away from the mirror. The degree to which this may be accomplished depends on the degree of control one has over the flux surface of the magnetic field, which changes as a result of its interaction with a moving plasma.

Thio, Y. C. Francis↗

Neutron elastic scattering on calcium isotopes from chiral nuclear optical potentials

We formulate microscopic neutron-nucleus optical potentials from many-body perturbation theory based on chiral two- and three-body forces. The neutron self-energy is first calculated in homogeneous matter to second order in perturbation theory, which gives the central real and imaginary terms of the optical potential. Here, the real spin-orbit term is calculated separately from the density matrix expansion using the same chiral interaction as in the self-energy. Finally, the full neutron-nucleus optical potential is derived within the improved local density approximation utilizing mean-field models consistent with the chiral nuclear force employed. We compare the results of the microscopic calculations to phenomenological models and experimental data up to projectile energies of E = 200 MeV. Experimental elastic differential scattering cross sections and vector analyzing powers are generally well reproduced by the chiral optical potential, but we find that total cross sections are overestimated at high energies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Efficient mapping between void shapes and stress fields using Deep Convolutional Neural Networks with sparse data

Establishing fast and accurate structure-to-property relationships is an important component in the design and discovery of advanced materials. Physics-based simulation models like the finite element method (FEM) are often used to predict deformation, stress, and strain fields as a function of material microstructure in material and structural systems. Such models may be computationally expensive and time intensive if the underlying physics of the system is complex. This limits their application to solve inverse design problems and identify structures that maximize performance. In such scenarios, surrogate models are employed to make the forward mapping computationally efficient to evaluate. However, the high dimensionality of the input microstructure and the output field of interest often renders such surrogate models inefficient, especially when dealing with sparse data. Deep convolutional neural network (CNN) based surrogate models have shown great promise in handling such high-dimensional problems. In this paper, a single ellipsoidal void structure under a uniaxial tensile load represented by a linear elastic, high-dimensional and expensive-to-query, FEM model. We consider two deep CNN architectures, a modified convolutional autoencoder framework with a fully connected bottleneck and a UNet CNN, and compare their accuracy in predicting the von Mises stress field for any given input void shape in the FEM model. Additionally, a sensitivity analysis study is performed using the two approaches, where the variation in the prediction accuracy on unseen test data is studied through numerical experiments by varying the number of training samples from 20 to 100.

surrogate modeling; convolutional neural networks;↗

Experimental Investigation of Low-Frequency Distributed Acoustic Strain-Rate Responses to Propagating Fractures

Summary Low-frequency distributed acoustic sensing (LF-DAS) responds to strain changes due to far-field fracture propagation. To better understand the LF-DAS response to fracture propagation, we performed laboratory-scale hydraulic fracture experiments with embedded optical strain sensors. The objectives of this research are to generate hydraulic fractures of known geometry, measure the strain response along the embedded fiber-optic cable comparable to LF-DAS measurements, and use the results to inform the interpretation of field-derived LF-DAS data. The experiments were conducted in unconfined transparent cubic blocks with a dimension of 8 in. on each side. The block was made of transparent epoxy to visualize the fracture propagation. Fiber-optic sensing cables were embedded in the block at different distances to the source of injection. We injected dyed water through an injection tubing to generate a transverse, radial fracture along an initial flaw. An optical interrogator recorded the response of offset fiber Bragg grating (FBG) strain sensors normal to the plane of the fracture. The strain data were visualized on a waterfall plot, akin to visualizations of field-derived LF-DAS data. Dimensional analysis was used to scale the laboratory results to field conditions. We compared the evolution of the strain response at the fiber-optic cable, injection pressure, and rate with known fracture geometry. The measured strains were compared with Sneddon’s (1946) linear elastic solution for a penny-shaped crack and found to follow this behavior. The generated radial fractures in transparent media can be modeled with Sneddon’s linear elastic radial fracture model and a Mode I critical stress intensity factor. The LF-DAS characteristic response of a narrowing region of extension surrounded by compression was exhibited as a fracture that approached and intersected the fiber-optic cable. The experimentally derived strain and strain-rate waterfall plots with known fracture geometry, injection rate, and pressure response provide insight in understanding LF-DAS responses in the field. Furthermore, we developed a method to estimate fracture geometry evolution from the fiber-optic strain data and validated the method against the experimental data.

Engineering↗

Stress Formulation in Three-Dimensional Elasticity

The theory of elasticity evolved over centuries through the contributions of eminent scientists like Cauchy, Navier, Hooke Saint Venant, and others. It was deemed complete when Saint Venant provided the strain formulation in 1860. However, unlike Cauchy, who addressed equilibrium in the field and on the boundary, the strain formulation was confined only to the field. Saint Venant overlooked the compatibility on the boundary. Because of this deficiency, a direct stress formulation could not be developed. Stress with traditional methods must be recovered by backcalculation: differentiating either the displacement or the stress function. We have addressed the compatibility on the boundary. Augmentation of these conditions has completed the stress formulation in elasticity, opening up a way for a direct determination of stress without the intermediate step of calculating the displacement or the stress function. This Completed Beltrami-Michell Formulation (CBMF) can be specialized to derive the traditional methods, but the reverse is not possible. Elasticity solutions must be verified for the compliance of the new equation because the boundary compatibility conditions expressed in terms of displacement are not trivially satisfied. This paper presents the variational derivation of the stress formulation, illustrates the method, examines attributes and benefits, and outlines the future course of research.

Patnaik, Surya N.↗

Unleashing the power of EFT in neutrino-nucleus scattering

Neutrino physics is advancing into a precision era with the construction of new experiments, particularly in the few GeV energy range. Within this energy range, neutrinos exhibit diverse interactions with nucleons and nuclei. This study delves in particular into neutrino-nucleus quasi-elastic cross sections, taking into account both standard and, for the first time, non-standard interactions, all within the framework of effective field theory (EFT). The main uncertainties in these cross sections stem from uncertainties in the nucleon-level form factors, and from the approximations necessary to solve the nuclear many-body problem. We explore how these uncertainties influence the potential of neutrino experiments to probe new physics introduced by left-handed, right-handed, scalar, pseudoscalar, and tensor interactions. For some of these interactions the cross section is enhanced, making long-baseline experiments an excellent place to search for them. Our results, including tabulated cross sections for all interaction types and all neutrino flavors, can serve as the foundation for such searches.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Graph neural networks for mechanical property prediction of 2D fiber composites

This work investigates the ability of graph neural networks (GNNs) to homogenize 2D fiber composite microstructures. We use different inhomogeneity and anisotropy indices to motivate and show that the Volume Elements (VEs) used in ML methods should ideally be far from their Representative Volume Element (RVE) size limit and, consequently, are notably anisotropic. Hence, training only the isotropic limit properties may not be acceptable. Another aspect is the need to normalize elastic stiffness values for ML, especially when high elastic contrast ratios are encountered between composite phases or in the material set. We introduce a normalization technique based on the mean-field method (MFM) to handle such high contrast ratios and train for the entire stiffness tensor. We show that the proposed GNN approaches exhibit high accuracy and efficiency compared to traditional methods and convolutional neural networks, utilizing unstructured graphs constructed from microstructure topology. Our model successfully predicts the stiffness tensor, peak strength under bulk damage, and brittle fracture initiation strength across diverse microstructure configurations while maintaining high accuracy even for extreme material contrasts and volume fractions. We also present a method to improve prediction accuracy for small dataset sizes using Voronoi partitioning.

Brittle strength↗

Monoenergetic neutron scattering on 238 U at Triangle Universities Nuclear Laboratory

In this report, we summarize the progress made in FY23 for experiments fielded at Triangle Universities Nuclear Laboratory (TUNL) of Duke University on the neutron scattering measurements of 238 U. The preliminary results on the elastic and inelastic reaction channels will be presented in addition to details of experiment and data analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Random walk theory applied to electron avalanche formation

Use of microscopic detail in random walk theory describing the initial formations of a large number of avalanches shows that concomitant electron transport coefficients quickly relax to equilibrium values. This enables the use of random walks having step sizes and probabilities based only on local electric field strengths and densities. A self-consistent avalanche solution which accounts for collective long range Coulomb interactions as well as short range elastic and inelastic collisions between electrons and background atoms is demonstrated for helium. Avalanche growth retardation followed by an abrupt growth augmentation as time proceeds is shown to be associated with the formation of regions of charge density extrema near the avalanche axis and within the axial distance covered by the electron swarm.

Englert, G. W.↗

SCF treatment of charge polarization effects in intermediate-energy electron scattering calculations with applications to N2

Converged rotational close coupling calculations of differential, integral, and momentum-transfer cross sections are obtained for seven model potentials for electron-N2 scattering at an impact energy of 30 eV. The model potentials involve a static potential calculated by the INDO/1s or INDOXI/1s method, and exchange potential calculated by the semiclassical exchange approximation from the INDO/1s or INDOXI/1s unperturbed electron density and a polarization potential. The polarization potentials used include the Buckley-Burke semiempirical one and various modifications of the INDOXI and INDO SCF (self-consistent-field) adiabatic polarization potentials. Without adjustment of parameters, good agreement is obtained with the angle dependence of the experimentally measured sum of the elastic and rotational excitation differential cross sections; however, the absolute value of the calculated cross sections is about 20-30% larger than the measured values in the best case, perhaps indicating that the model potentials are too strong or should have a nonzero imaginary part.

Onda, K.↗

Advances and trends in plate buckling research

Advances and current trends in plate bucklin research are summarized. The field is divided into three parts: (1) classical buckling studies, including plates of rectangular, circular and other shapes; (2) classical complicating effects, including elastic foundation, anisotropic material, variable thickness, shear deformation and nonhomogeneous material; and (3) nonclassical considerations, including postbuckling, imperfections, parametric excitation, follower forces, magnetoelastic buckling and inelastic buckling.

Leissa, A.↗

Wave propagation in intense flux tubes

The nature of nonadiabatic wave propagation in a slender magnetic flux tube is explored. The results of the theory are compared with the observations of Giovanelli et al. (1978), and found to be in general agreement. Those observations, of tubes in the photosphere-chromosphere, show outwardly propagating waves, with periods of 300 s, which take some 19 s to propagate from one level of line formation to another level higher in the atmosphere. In sharp contrast to this, is the time of 7 s for a similar disturbance outside the tube to propagate between the same two levels of line formation, estimated to be some 600 km apart in the field-free atmosphere. It is argued that the sharply contrasting propagation times for the tube and its environment is principally due to the elasticity of the tube and its subsequent propensity for propagation. A nonadiabatic disturbance may be essentially propagating within the tube but essentially nonpropagating outside, with considerably slower phase speeds thus arising inside the tube. The theory suggests that the observed disturbances are nonadiabatic, acoustic-gravity waves channelled along a magnetic flux tube and modulated by external pressure variations.

Roberts, B.↗

Modeling of periodic great earthquakes on the San Andreas fault: Effects of nonlinear crustal rheology

We analyze the cycle of great earthquakes along the San Andreas fault with a finite element numerical model of deformation in a crust with a nonlinear viscoelastic rheology. The viscous component of deformation has an effective viscosity that depends exponentially on the inverse absolute temperature and nonlinearity on the shear stress; the elastic deformation is linear. Crustal thickness and temperature are constrained by seismic and heat flow data for California. The models are for anti plane strain in a 25-km-thick crustal layer having a very long, vertical strike-slip fault; the crustal block extends 250 km to either side of the fault. During the earthquake cycle that lasts 160 years, a constant plate velocity v(sub p)/2 = 17.5 mm yr is applied to the base of the crust and to the vertical end of the crustal block 250 km away from the fault. The upper half of the fault is locked during the interseismic period, while its lower half slips at the constant plate velocity. The locked part of the fault is moved abruptly 2.8 m every 160 years to simulate great earthquakes. The results are sensitive to crustal rheology. Models with quartzite-like rheology display profound transient stages in the velocity, displacement, and stress fields. The predicted transient zone extends about 3-4 times the crustal thickness on each side of the fault, significantly wider than the zone of deformation in elastic models. Models with diabase-like rheology behave similarly to elastic models and exhibit no transient stages. The model predictions are compared with geodetic observations of fault-parallel velocities in northern and central California and local rates of shear strain along the San Andreas fault. The observations are best fit by models which are 10-100 times less viscous than a quartzite-like rheology. Since the lower crust in California is composed of intermediate to mafic rocks, the present result suggests that the in situ viscosity of the crustal rock is orders of magnitude less the rock viscosity determined in the laboratory.

Reches, Ze'ev↗

Free-field Ground Motion Induced by Underground Explosions at Aqueduct Mesa with Predictions for Physical Experiment One (PE1)

We fit standard power-law models for ground acceleration and particle velocity using scaled ground motion observations from underground nuclear explosions at Aqueduct Mesa on the Nevada National Security Site. The models are then validated using small-scale high-explosive (HE) test data from nearby Rainier Mesa. The comparison between model predictions and observations is good and can be made more favorable if the assumed yield for the HE tests are doubled. Further validation is made with a numerical experiment that additionally shows a transition in ground motion attenuation to elastic propagation (1 = r). An extended model is provided to incorporate the long-range transition in attenuation, which results in increased ground motions at farther range. We then make predictions for the field experiment, Physical Experiment One (PE1), to take place in P-tunnel on Aqueduct Mesa, and suggest that the models could be used in planning various operations for PE1.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

An evaluation of a coupled microstructural approach for the analysis of functionally graded composites via the finite-element method

A comparison is presented between the predictions of the finite-element analysis and a recently developed higher-order theory for functionally graded materials subjected to a thorough-thickness temperature gradient. In contrast to existing micromechanical theories that utilize classical (i.e., uncoupled) homogenization schemes to calculate micro-level and macro-level stress and displacement fields in materials with uniform or nonuniform fiber spacing (i.e., functionally graded materials), the new theory explicitly couples the microstructural details with the macrostructure of the composite. Previous thermo-elastic analysis has demonstrated that such coupling is necessary when: the temperature gradient is large with respect to the dimension of the reinforcement; the characteristic dimension of the reinforcement is large relative to the global dimensions of the composite and the number of reinforcing fibers or inclusions is small. In these circumstances, the standard micromechanical analyses based on the concept of the representative volume element used to determine average composite properties produce questionable results. The comparison between the predictions of the finite-element method and the higher-order theory presented herein establish the theory's accuracy in predicting thermal and stress fields within composites with a finite number of fibers in the thickness direction subjected to a thorough-thickness thermal gradient.

Pindera, Marek-Jerzy↗