Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite amplitude method”

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.

281 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Reactant conversion in homogeneous turbulence: Mathematical modeling, computational validations and practical applications

Closed form analytical expressions are obtained for predicting the limited rate of reactant conversion in a binary reaction of the type F + rO yields (1 + r) Product in unpremixed homogeneous turbulence. These relations are obtained by means of a single point Probability Density Function (PDF) method based on the Amplitude Mapping Closure. It is demonstrated that with this model, the maximum rate of the reactants' decay can be conveniently expressed in terms of definite integrals of the Parabolic Cylinder Functions. For the cases with complete initial segregation, it is shown that the results agree very closely with those predicted by employing a Beta density of the first kind for an appropriately defined Shvab-Zeldovich scalar variable. With this assumption, the final results can also be expressed in terms of closed form analytical expressions which are based on the Incomplete Beta Functions. With both models, the dependence of the results on the stoichiometric coefficient and the equivalence ratio can be expressed in an explicit manner. For a stoichiometric mixture, the analytical results simplify significantly. In the mapping closure, these results are expressed in terms of simple trigonometric functions. For the Beta density model, they are in the form of Gamma Functions. In all the cases considered, the results are shown to agree well with data generated by Direct Numerical Simulations (DNS). Due to the simplicity of these expressions and because of nice mathematical features of the Parabolic Cylinder and the Incomplete Beta Functions, these models are recommended for estimating the limiting rate of reactant conversion in homogeneous reacting flows. These results also provide useful insights in assessing the extent of validity of turbulence closures in the modeling of unpremixed reacting flows. Some discussions are provided on the extension of the model for treating more complicated reacting systems including realistic kinetics schemes and multi-scalar mixing with finite rate chemical reactions in more complex configurations.

Madnia, C. K.↗
DOE OSTI2020

Unitary Bethe-Salpeter Methods in Two- and Three-Body Systems

The study of resonances, unstable particles formed in particle scattering, is motivated by the questions about the strong interaction and the composition of hadrons. In order to advance this study in the intermediate energy region, a variety of techniques are used. In this thesis, we present work done on two different analyses that use the principle of unitarity and the Bethe-Salpeter equation to study the resonances Lambda(1405), Sigma(1385), and a1(1260). Firstly, we perform a simultaneous analysis of Sand P-waves of the Strangeness = -1 two-body meson-baryon scattering amplitude using all low-energy data. For the first time, differential cross section data are included for chiral unitary coupled-channel models. From this model, S- and P-wave amplitudes are extracted and we observe both well-known I(JP)=0(1/2^-) S-wave states as well as a new I(JP)=1(1/2^+) state absent in quark models and lattice QCD results. Multiple statistical and phenomenological tests suggest that, while the data clearly require an I =1 P-wave resonance, the new state just accounts for the absence of the decuplet Sigma(1385)3/2+ in the model which is subsequently included in the parameterization. Secondly, we formulate the final state interaction of the a1(1260) resonance decay in a manifestly three-body-unitary parameterization and fit it to the a1(1260) lineshape measured by the ALEPH experiment. Dalitz plots calculated from this fit are presented. The work demonstrates the feasibility to numerically solve a previously derived amplitude and its generalization to isobars with spin and coupled channels. The model is a good test case because it can also be applied to other meson decays including exotic states and modified for the finite-volume problem as it arises in lattice QCD due to its manifest unitarity. The understanding of the resonances being analyzed in both the two-body work and the three-body work gives a better understanding of intermediate energy QCD.

Sadasivan, Daniel↗
DOE OSTI2021

Real-time quantum calculations of phase shifts using wave packet time delays

We present a method to extract the phase shift of a scattering process using the real-time evolution in the early and intermediate stages of the collision in order to estimate the time delay of a wave packet. This procedure is convenient when using noisy quantum computers for which the asymptotic out-state behavior is unreachable. We demonstrate that the challenging Fourier transforms involved in the state preparation and measurements can be implemented in 1+1 dimensions with current trapped ion devices and IBM quantum computers. We compare quantum computations of the time delays obtained in the one-particle quantum mechanics limit and the scalable quantum field theory formulation with accurate numerical results. We discuss the finite volume effects in the Wigner formula connecting time delays to phase shifts. The results reported involve two- and four-qubit calculations, and we discuss the possibility of larger scale computations in the near future.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Nonlinear evolution, propagation, electron-trapping, and damping effects of ion-acoustic solitons using fully kinetic PIC simulations

We investigate ion acoustic solitary waves (solitons) of varying amplitudes in a one-dimensional plasma using fully kinetic particle-in-cell simulations. The initial soliton conditions are based on the Korteweg–de Vries (KdV) equation, treating ions as a cold species and electrons with finite temperature. Our findings reveal that KdV solitons evolve nonlinearly to a saturated state at higher amplitude, deviating from KdV predictions for ion density and electric potential, and from the Boltzmann relation for electron density. At this saturated state, the KdV model cannot accurately describe the soliton behavior. For small amplitudes, Sagdeev's model describes the saturated state, but not the soliton width; for larger amplitudes, it models the width accurately, but not the amplitude. These discrepancies arise from assuming a Boltzmann relation for electron density, while electron trapping creates non-Boltzmann densities—a deviation that increases with soliton amplitude. Additionally, we observe that the soliton amplitude oscillates roughly at the electron bounce frequency. The soliton is better described by Schamel's electron density formulation and a modified KdV equation incorporating electron trapping. The soliton velocity matches best with predictions from Sagdeev's and Schamel's models. Moreover, the soliton speed–amplitude relationship differs from existing theoretical predictions. Finally, we find minimal ion and electron Landau damping effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2018

Global High Resolution Crustal Magnetic Field at the Surface of the Moon from Low-Altitude Lunar Prospector Magnetic Gradient Data

We derived new vector gradients based models of crustal magnetic field at the lunar surface with data from the Lunar Prospector (LP) satellite using two model parameterization approaches: a global set of 35820 1° spaced (~30 km) equal area monopoles at 20 km below the surface (O’Brien and Parker, 1994; Olsen et al., 2017) and combined results of subsets of 100000 0.66° spaced monopoles at the same depth. We use the scheme of iteratively reweighted least-squares inversion to compute the initial model. Then the amplitudes of these monopoles are determined by minimizing the misfit to the components together with the global average of |Br| at the ellipsoid surface (i.e. applying a L1 model regularization of Br). In previous approaches using vector fields for modeling, we found that external field contamination leads to spurious anomalies in the downward continued field models even with stringent data selection criteria and ad-hoc noise removal techniques (e.g., satellite’s position in the Moon’s wake w.r.t. the solar wind and in the Earth’s magnetotail, internal/external dipoles fields removal, low-order polynomial removal, joint equivalent source cross-validation technique and visually removing remaining anomalous segments). On the other hand, with the use of gradients-only data (along-track first finite differences), we were able to completely bypass the ad-hoc techniques. Similar processing of Kaguya magnetic data, which have only higher altitude coverage over the most of the Moon except in the region of the South Pole–Aitken (SPA) basin, completely misses some of the anomalies seen in the Lunar Prospector data. The combined Lunar Prospector and Kaguya gradient-based models also severely degrade the derivation of the anomaly fields in many regions. Euler analysis of isolated anomaly features from the Reiner Gamma swirl suggests top depths of about 0.3 to 1.5 km and center depths of 10-14 km; in the region between Stein and Vallier craters north of the SPA basin our analysis suggests top depths of around 1.5 km and center depths of 13-15 km. With the spectral depth determination techniques, the SPA basin region yields depths to the base of magnetization ranging between 15 and 40 km. Three-dimensional modeling and the bulk magnetization determinations of the sources constrained by the Euler and spectral methods is underway.

Gradient↗
NASA NTRS2004

Reduced-Order Blade Mistuning Analysis Techniques Developed for the Robust Design of Engine Rotors

The primary objective of this research program is to develop vibration analysis tools, design tools, and design strategies to significantly improve the safety and robustness of turbine engine rotors. Bladed disks in turbine engines always feature small, random blade-to-blade differences, or mistuning. Mistuning can lead to a dramatic increase in blade forced-response amplitudes and stresses. Ultimately, this results in high-cycle fatigue, which is a major safety and cost concern. In this research program, the necessary steps will be taken to transform a state-of-the-art vibration analysis tool, the Turbo-Reduce forced-response prediction code, into an effective design tool by enhancing and extending the underlying modeling and analysis methods. Furthermore, novel techniques will be developed to assess the safety of a given design. In particular, a procedure will be established for using eigenfrequency curve veerings to identify "danger zones" in the operating conditions--ranges of rotational speeds and engine orders in which there is a great risk that the rotor blades will suffer high stresses. This work also will aid statistical studies of the forced response by reducing the necessary number of simulations. Finally, new strategies for improving the design of rotors will be pursued. Several methods will be investigated, including the use of intentional mistuning patterns to mitigate the harmful effects of random mistuning, and the modification of disk stiffness to avoid reaching critical values of interblade coupling in the desired operating range. Recent research progress is summarized in the following paragraphs. First, significant progress was made in the development of the component mode mistuning (CMM) and static mode compensation (SMC) methods for reduced-order modeling of mistuned bladed disks (see the following figure). The CMM method has been formalized and extended to allow a general treatment of mistuning. In addition, CMM allows individual mode mistuning, which accounts for the realistic effects of local variations in blade properties that lead to different mistuning values for different mode types (e.g., mistuning of the first torsion mode versus the second flexural mode). The accuracy and efficiency of the CMM method and the corresponding Turbo-Reduce code were validated for an example finite element model of a bladed disk.

Min, James B.↗
DOE OSTI2025

High-Speed and High-Quality Field Welding Repair Based on Advanced Non-Destructive Evaluation and Numerical Modeling

Creep strength-enhanced ferritic (CSEF) steels such as Grade 91 (9Cr-1Mo-V) and Grade 92 (Fe-9Cr-2W-0.5Mo) steels are widely used in the fossil-fuel-fired and nuclear power plants. The weld integrity of these steels is crucial for power plants' safe and reliable operations. Due to harsh service conditions, the steel weld can become susceptible to environmental degradation. Field welding repair is used to restore the degraded weld’s performance where a controlled temper-bead welding technique is commonly used to temper the freshly formed martensite during welding. However, knowledge of weld repairability is limited and experimental trial and error optimization to achieve desired microstructure and joint properties is expensive and time-consuming. Many existing computational models, e.g., finite element models, are limited to solving heat conduction equation and ignoring convective heat transfer due to molten metal flow. These models can result in over-prediction of peak temperatures of weld pool and heat-affected zone (HAZ), which in turn can affect the accuracy of tempering prediction. Moreover, these finite element models require an input of the deposit profiles in advance and thus limits the usability of these models. Here, a molten pool-based, multi-pass multi-layer model has been developed based on computational fluid dynamics (CFD) approach with the Volume of Fluid (VOF) method. The model calculates the bead formation, thereby eliminating the need for pre-determined bead profiles required by finite element models. For computational efficiency, a coordinate system attached to the moving heat source is utilized. A subroutine is developed to convert the temperature profiles in the reference frame stationary to the heat source to that stationary to the workpiece. The converted thermal cycles are then imported into a microstructure model to compute the tempering kinetics and resultant hardness using a Johnson-Mehl-Avrami-Kolmogorov (JMAK), and modified Grange-Baughman parameter. The modeling approach is first developed and validated on single- and multi-pass deposition of stainless steel filler metal onto a SA-533 high strength steel substrate. The models are then applied to a multi-pass V-groove repair weld of Grade 91 steel plate as well as directed energy deposition of Grade 92 steel. Non-destructive characterization of microstructures was performed on Grade 91 and 92 steel welds. Two welding processes, cold metal transfer (CMT) and flux-cored arc welding (FCAW), were investigated for the Grade 91 steel weld samples. For the Grade 92 weld samples, three different heat inputs (low, medium, and high) of gas tungsten arc welding (GTAW) were utilized to replicate traditional field welding processes. The non-destructive evaluation (NDE) method used for this research was immersion ultrasonic testing (UT) using a micro-resolution ultrasonic imaging methodology specifically designed to operate in the through-transmission configuration operating at 20 MHz of frequency. The system used a focused ultrasonic beam spot size diameter between 250-300 μm, and a 6 μm laser vibrometer spot size for detection, to produce highly defined images with longitudinal and mode-converted shear waves. From the micro-resolution ultrasonic C-scan images, three microstructural regions, i.e., weld metal (WM), HAZ, and base metal (BM), were clearly identifiable. Various levels of ultrasonic amplitudes distributed over the three regions were correlated with electron beam backscattered diffraction (EBSD) images using grain size, grain boundaries, and dislocation densities. The results showed that areas with relatively higher ultrasonic amplitude levels were associated with smaller grains and higher dislocation densities, while areas with lower amplitude levels were associated with larger grains and lower dislocation densities. In addition, ultrasonic velocity data obtained across the three different weld microstructural regions of Grade 91 test samples were correlated with optical metallographic images and hardness measurements. The results showed distinctive decreases in ultrasonic velocity and hardness over the HAZ region, where weld failures often occur during service.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Modulation of intensity emerging from zero effort (MIEZE) with extended Fourier time at large scattering angle

Modulation of Intensity Emerging from Zero Effort (MIEZE) is a neutron resonant spin echo technique that allows one to measure time correlation scattering functions in materials by implementing radio-frequency (RF) intensity modulation at the sample and the detector. The technique avoids neutron spin manipulation between the sample and the detector and, thus, could find applications in cases where the sample depolarizes the neutron beam. However, the finite sample size creates a variance in the path length between the locations where scattering and detection happen, which limits the contrast in intensity modulation that one can detect, in particular, toward long correlation times or large scattering angles. We propose a modification to the MIEZE setup that will enable one to extend those detection limits to longer times and larger angles. Furthermore, we use Monte Carlo simulations of a neutron scattering beamline to show that by tilting the RF flippers in the primary spectrometer with respect to the beam direction, one can shape the wave front of the intensity modulation at the sample to compensate for the path variance from the sample and the detector. The simulation results indicate that this change enables one to operate a MIEZE instrument at much increased RF frequencies, thus improving the effective energy resolution of the technique. For the MIEZE instrument simulated, it shows that for an incident beam with the maximum divergence of 0.33°, the maximum Fourier time can be increased by a factor of 3.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

A combined numerical and experimental study to elucidate primary breakup dynamics in liquid metal droplet-on-demand printing

Droplet-on-demand liquid metal jetting is emerging as a powerful technology for the additive manufacturing of metallic parts. The success of this method hinges on overcoming several technological challenges. Additionally, the principal one among these challenges is the controlled repeatable ejection of single uniform droplets. Due to the high density and surface tension of liquid metals, the droplet ejection process occurs near the minimal extremes of the printability phase diagram, defined by acceptable ranges for the Weber (We) and Ohnesorge (Oh) numbers. In this work, we experimentally demonstrate the satellite-free ejection of pneumatically actuated molten tin droplets in this extreme corner of printability and use a combination of high-speed video analysis and volume-of-fluid modeling to elucidate the droplet dynamics. While the simulations at low Oh and We can correctly describe several aspects of the breakup process, such as an increasing tail and pinch-point near the nozzle, no single parameter set can completely capture the droplet shape at breakup. Instead, the experimental droplet dynamics appear to include features from both high and low Oh breakup. This disagreement is ascribed to the incomplete description of the droplet ejection process including wetting and exit effects near the nozzle opening and surface effects such as transient cooling and oxide formation.

36 MATERIALS SCIENCE↗
NASA NTRS2013

Engineering Evaluation of International Low Impact Docking System Latch Hooks

The international Low Impact Docking System (iLIDS) provides a structural arrangement that allows for visiting vehicles to dock with the International Space Station (ISS) (Fig 1). The iLIDS docking units are mechanically joined together by a series of active and passive latch hooks. In order to preserve docking capability at the existing Russian docking interfaces, the iLIDS latch hooks are required to conform to the existing Russian design. The latch hooks are classified as being fail-safe. Since the latch hooks are fail-safe, the hooks are not fracture critical and a fatigue based service life assessment will satisfy the structural integrity requirements. Constant amplitude fatigue testing to failure on four sets of active/passive iLIDS latch hooks was performed at load magnitudes of 10, 11, and 12 kips. Failure analysis of the hook fatigue failures identified multi-site fatigue initiation that was effectively centered about the hook mid-plane (consistent with the 3D model results). The fatigue crack initiation distribution implies that the fatigue damage accumulation effectively results in a very low aspect ratio surface crack (which can be simulated as thru-thickness crack). Fatigue damage progression resulted in numerous close proximity fatigue crack initiation sites. It was not possible to determine if fatigue crack coalescence occurs during cyclic loading or as result of the fast fracture response. The presence of multiple fatigue crack initiation sites on different planes will result in the formation of ratchet marks as the cracks coalesce. Once the stable fatigue crack becomes unstable and the fast fracture advances across the remaining ligament and the plane stress condition at a free-surface will result in failure along a 45 deg. shear plane (slant fracture) and the resulting inclined edge is called a shear lip. The hook thickness on the plane of fatigue crack initiation is 0.787". The distance between the shear lips on this plane was on the order of 0.48" and it was effectively centered about the mid-plane of the section. The numerous ratchet marks between the shear lips on the fracture initiation plane are indicative of multiple fatigue initiation sites within this region. The distribution of the fatigue damage about the centerline of the hook is consistent with the analytical results that demonstrate peak stress/strain response at the mid-plane that decreases in the direction of the hook outer surfaces. Scanning electron microscope images of the failed sections detected fatigue crack striations in close proximity to the free surface of the hook radius. These findings were documented at three locations on the fracture surface : 1) adjacent to the left shear lip, 2) adjacent to the right shear lip, and 3) near the centerline of the section. The features of the titanium fracture surface did not allow for a determination of a critical crack size via identification of the region where the fatigue crack propagation became unstable. The fracture based service life projections where benchmarked with strain-life analyses. The strainrange response in the hook radius was defined via the correlated finite element models and the modified method of universal slopes was incorporated to define the strain-life equation for the titanium alloy. The strain-life assessment confirmed that the fracture based projections were reasonable for the loading range of interest. Based upon the analysis and component level fatigue test data a preliminary service life capability for the iLIDS active and passive hooks of 2 lifetimes is projected (includes a scatter factor of 4).

Martinez, J.↗
NASA NTRS2018

Implementation of J-A Methodology Elastic-Plastic Crack Instability Analysis Capability into the WARP-3D Code

Characterization of the near crack-tip stress/strain fields is the foundation of fracture mechanics. The description of the near tip stress field and the prediction of when fracture occurs is well established for brittle materials that exhibit linear elastic behavior. However, in ductile materials or conditions that violate linear elastic assumptions (Aluminum alloys, Al 2024-T3, Al 2024- T351 etc.), the elastic-plastic crack-tip stress fields are characterized by the Hutchison-Rice-Rosengren (HRR) field. The J-integral is commonly used to characterize amplitude of the HRR field under elastic-plastic conditions. The J-integral has been demonstrated for crack-tip fields that are under high constraint conditions (i.e., small-scale plasticity where the J-dominance is maintained). However, as the external load increases, yielding changes from small- to largescale plasticity and usually a loss of constraint (i.e., reduction in the triaxial stress field along the crack front). The loss of constraint leads to the deviation of the crack-tip stress fields from that given by the HRR field. Hence, the J-dominance will be gradually lost and additional parameter(s) are required to quantify the crack-tip stress fields and predict fracture behavior. The assessment objectives were to: 1) implement a two-parameter (i.e., J-A) fracture criterion into an elastic-plastic three-dimensional (3D) finite element analysis (FEA), 2) validate the implementation by comparison with the A parameter from literature data, 3) conduct material characterization tests to quantify the material behavior and provide fracture data for validation of the J-A fracture criteria, and (4) perform evaluations to establish if the J-A criteria can be used to predict fracture in a ductile metallic material (e.g., aluminum alloys). The A parameter in these criteria is the second parameter in a three-term elastic-plastic asymptotic expansion of the neartip stress behavior. A series of extensive FEAs were performed using WARP3D software package to obtain solutions for the A parameter for different specimen configurations. The methodology needed for the estimation of the A parameter in the asymptotic expansion was developed and implemented using Matlab®. A user material (UMAT) routine was used to model the material stress-strain response using a Ramberg-Osgood power law with a hardening exponent (n) and a material coefficient (alpha). This UMAT routine was successfully implemented in WARP3D software and validated through comparison with the experimental data. Three configurations were extracted from published results: 1) center cracked plate (CCP), 2) single edge-cracked plate (SECP), and 3) double edge-cracked plate (DECP). These configurations and four other configurations (three-hole tension (THT)), three-point bend (3PTB), three-hole compact tension (3PCT), and compact tension (CT)) were analyzed to verify the methodology that was developed and implemented into WARP3D. Solutions of the A parameter were obtained for remote tension loading conditions that started with small-scale yielding and continued into the large-scale plasticity regime. The results indicate that the methodology developed can be used to calculate the elastic-plastic J-A parameters for test specimens with a range of crack geometries, material strain hardening behaviors, and loading conditions. The J-A parameters were implemented as fracture criteria and used to predict the test results. For comparison, other fracture criteria were used to predict the same test results. Major findings include: The A constraint parameter A varies with specimen type and applied load thus accurate determination is crucial in predicting the failure load, and the A parameter is asymptotic as the failure load is approached, making an accurate determination difficult (i.e., small differences in the A parameter can cause large variations in failure load) for materials exhibiting elastic-plastic behavior. The failure predictions from J-A methodology were more accurate than the traditionally used KC and J methods, and have comparable scatter to that observed when using the crack-tip opening angle (CTOA) method. However, the J-A methodology requires considerable effort (expertise level and labor) to implement and to evaluate the A parameter for different specimen types and materials, or to apply this methodology to part-through crack (e.g., 3D problems) structural applications.

Hamm, Kenneth R., Jr.↗