Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Variational approach”

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 847 records · Page 47

Machine Learning Inference of Random Medium Properties

Earth materials are heterogeneous across a range of spatial scales, but the resolvability of small structures is limited by sparse data coverage, noise, bandlimitedness, and other difficulties. In practice, heterogeneities below a certain size cannot be recovered from seismic data except through statistical medium descriptions, which even then can be difficult to uniquely determine. To improve the characterization of such heterogeneities, we develop a novel supervised machine learning (ML) model that provides insight about the recoverability of statistical medium properties from elastic waveform data and succeeds despite cycle-skipping and other challenges well known from elastic waveform inversion. We demonstrate the approach using random media generated by superimposing self-affine random variations on homogeneous and layered background structures. After training on sparsely-recorded, high-frequency waveforms from hundreds of different random medium realizations, we show the ability of our ML model to recover correlation lengths and other statistical properties of interest to near-surface and crustal seismology, among other fields. For frequency passbands and spatial offsets encountered in seismology, Gaussian correlation lengths and the amplitude of the random variations relative to the background model are recovered even in challenging scenarios involving unknown medium parameters, complex crustal structures, and low signal-to-noise ratio. In comparison, von Kármán correlation lengths, which are related to larger-wavelength variations of the medium than Gaussian correlation lengths, are not as well recovered. These results provide one of the first and most systematic investigations of the recoverability of statistical properties of heterogeneities below the resolution limit of deterministic seismic tomography, and suggest practical ML strategies for high-frequency waveform seismology.

58 GEOSCIENCES↗

Simulation Tools for Characterizing Stress Distribution in Laser Welded Dissimilar Joints

This project focuses on developing a thermo-metallurgical-mechanical modeling method to accurately predict the microstructural evolution and residual stress in laser welding between dissimilar metals, such as HSLA steel and high carbon equivalent (CE) gear steel. The method leverages a comprehensive material database to model the temperature and rate dependent phase transformations, along with their associated effects on material properties, such as thermal expansion and flow stress, throughout the welding process. A key innovation is the incorporation of phase transformation and phase-specific properties, which enhances the accuracy of residual stress predictions. The mixture material in the fusion zone due to the dissimilar metals will also be addressed in the numerical model. This is especially critical in scenarios involving phase transformations in the fusion zone and heat-affected zone (HAZ), where the phase changes can induce substantial residual stress variations. The material database has been generated using JMatPro. The modeling approach is implemented through a custom User Material (UMAT) subroutine, executed with the commercial finite element software Abaqus.

36 MATERIALS SCIENCE↗

Simulation Tools for Characterizing Stress Distribution in Laser Welded Dissimilar Joints

This project focuses on developing a thermo-metallurgical-mechanical modeling method to accurately predict the microstructural evolution and residual stress in laser welding between dissimilar metals, such as HSLA steel and high carbon equivalent (CE) gear steel. The method leverages a comprehensive material database to model the temperature and rate dependent phase transformations, along with their associated effects on material properties, such as thermal expansion and flow stress, throughout the welding process. A key innovation is the incorporation of phase transformation and phase-specific properties, which enhances the accuracy of residual stress predictions. The mixture material in the fusion zone due to the dissimilar metals will also be addressed in the numerical model. This is especially critical in scenarios involving phase transformations in the fusion zone and heat-affected zone (HAZ), where the phase changes can induce substantial residual stress variations. The material database has been generated using JMatPro. The modeling approach is implemented through a custom User Material (UMAT) subroutine, executed with the commercial finite element software Abaqus.

36 MATERIALS SCIENCE↗

Assessment of Self-interrogation Safeguards Signatures for Pebble Bed Reactor Fuel

New reactor designs, such as pebble-bed reactors, present challenges for the safeguarding of fissile material. Due to the multi-pass fuel circulation design and the non-uniform path a pebble may traverse through a pebble-bed reactor, there will be variations in the irradiation history of spent fuel pebbles. The standard approach of estimating actinide quantities based on fission gamma spectra or neutron emission paired with depletion modeling may yield uncertainties too large for safeguards and material accountancy. This project investigated, through modeling and experimentation, the potential of neutron self-interrogation of spent fuel pebbles as an innovative method to implement materials accountability. As an example, our feasibility studies indicate that the mass of U-235, U-238, Pu-239 and Pu-241 can be predicted to 4.1, 0.86, 13 and 13 % accuracy, respectively, when measuring 100 closely packed end-of-life spent fuel pebbles over approximately 12 days using a 4π counting geometry.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Exploring Nonmotorized Travel in New York State Using 2017 National Household Travel Survey Data

This report presents systematic analysis of nonmotorized travel behavior within New York State (NYS) utilizing 2017 National Household Travel Survey data. As walking and biking assume increasingly prominent roles in advancing active transportation objectives, accessible mobility, and sustainable transportation systems, comprehensive understanding of these modal patterns becomes essential for evidence-based policy development and strategic planning initiatives. Recognizing the substantial geographic, demographic, and socioeconomic heterogeneity characterizing NYS—spanning from the concentrated urban fabric of New York City (NYC) to dispersed suburban and rural contexts—this analysis employs the structured "4Ws" analytical framework to examine participant demographics, spatial and temporal distributions, and motivational factors underlying active transportation choices relative to motorized alternatives. This methodological approach captures modal behavior patterns, user characteristics, trip purposes, and temporal variations across the state's diverse contexts. Additionally, the research examines behavioral differences among distinct user classifications, including walk-only and bike-only travelers, while conducting comparative analysis of nonmotorized travel patterns between ALICE (Asset Limited, Income Constrained, Employed) and non-ALICE household categories.

99 GENERAL AND MISCELLANEOUS↗

Real-Time Krylov Theory for Quantum Computing Algorithms

Quantum computers provide new avenues to access ground and excited state properties of systems otherwise difficult to simulate on classical hardware. New approaches using subspaces generated by real-time evolution have shown efficiency in extracting eigenstate information, but the full capabilities of such approaches are still not understood. In recent work, we developed the variational quantum phase estimation (VQPE) method, a compact and efficient real-time algorithm to extract eigenvalues on quantum hardware. Here we build on that work by theoretically and numerically exploring a generalized Krylov scheme where the Krylov subspace is constructed through a parametrized real-time evolution, which applies to the VQPE algorithm as well as others. We establish an error bound that justifies the fast convergence of our spectral approximation. We also derive how the overlap with high energy eigenstates becomes suppressed from real-time subspace diagonalization and we visualize the process that shows the signature phase cancellations at specific eigenenergies. We investigate various algorithm implementations and consider performance when stochasticity is added to the target Hamiltonian in the form of spectral statistics. To demonstrate the practicality of such real-time evolution, we discuss its application to fundamental problems in quantum computation such as electronic structure predictions for strongly correlated systems.

97 MATHEMATICS AND COMPUTING↗

Problem-tailored Simulation of Energy Transport on Noisy Quantum Computers

The transport of conserved quantities like spin and charge is fundamental to characterizing the behavior of quantum many-body systems. Numerically simulating such dynamics is generically challenging, which motivates the consideration of quantum computing strategies. However, the relatively high gate errors and limited coherence times of today's quantum computers pose their own challenge, highlighting the need to be frugal with quantum resources. In this work we report simulations on quantum hardware of infinite-temperature energy transport in the mixed-field Ising chain, a paradigmatic many-body system that can exhibit a range of transport behaviors at intermediate times. We consider a chain with L = 12 sites and find results broadly consistent with those from ideal circuit simulators over 90 Trotter steps, containing up to 990 entangling gates. To obtain these results, we use two key problem-tailored insights. First, we identify a convenient basis – the Pauli Y basis – in which to sample the infinite-temperature trace and provide theoretical and numerical justifications for its efficiency relative to, e.g., the computational basis. Second, in addition to a variety of problem-agnostic error mitigation strategies, we employ a renormalization strategy that compensates for global nonconservation of energy due to device noise. We discuss the applicability of the proposed sampling approach beyond the mixed-field Ising chain and formulate a variational method to search for a sampling basis with small sample-to-sample fluctuations for an arbitrary Hamiltonian. This opens the door to applying these techniques in more general models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Flight investigation of the roll requirements for transport airplanes in the landing approach

An in-flight evaluation of transport roll characteristics in the landing approach was made with a general purpose airborne simulator. The evaluation task consisted of an instrument approach with a visual correction for a (200-foot) lateral offset. Pilot evaluations and ratings were obtained for approaches made at 140 knots and 180 knots indicated airspeed with variations of wheel characteristics, maximum roll rate, and roll time constant.

Holleman, E. C.↗

Study of magnetic perturbations on SEC vidicon tubes

A laboratory measurements program was conducted to determine the tolerances that must be imposed to achieve optimum performance from SEC-vidicon data sensors in the LST mission. These measurements along with other data were used to formulate recommendations regarding the necessary telemetry and remote control for the television data sensors when in orbit. The study encompassed the following tasks: (1) Conducted laboratory measurements of the perturbations which an external magnetic field produces on a magnetically focused, SEC-vidicon. Evaluated shielding approaches. (2) Experimentally evaluated the effects produced on overall performance by variations of the tube electrode potentials, and the focus, deflection and alignment fields. (3) Recommended the extent of ground control of camera parameters and camera parameter telemetry required for optimizing the performance of the television system in orbit. The experimental data are summarized in a set of graphs.

Long, D. C.↗

Advanced microwave soil moisture studies

Comparisons of low level L-band brightness temperature (TB) and thermal infrared (TIR) data as well as the following data sets: soil map and land cover data; direct soil moisture measurement; and a computer generated contour map were statistically evaluated using regression analysis and linear discriminant analysis. Regression analysis of footprint data shows that statistical groupings of ground variables (soil features and land cover) hold promise for qualitative assessment of soil moisture and for reducing variance within the sampling space. Dry conditions appear to be more conductive to producing meaningful statistics than wet conditions. Regression analysis using field averaged TB and TIR data did not approach the higher sq R values obtained using within-field variations. The linear discriminant analysis indicates some capacity to distinguish categories with the results being somewhat better on a field basis than a footprint basis.

Dalsted, K. J.↗

Comparison of Euler and full potential marching techniques for flows over complex configurations

Two recently developed aerodynamic prediction techniques based on the steady full potential equation and the unsteady Euler equations have been applied to a variety of three-dimensional supersonic flow problems exhibiting embedded subsonic regions. Both techniques utilize planar Gauss-Seidel relaxation in the marching direction and approximate factorization in the cross-flow plane. A conservative switching scheme and flux bias technique are employed in the full potential method to transition from the supersonic marching procedure to a subsonic relaxation algorithm and vice versa. A new unified approach with finite volume, high accuracy (up to third order) Total Variation Diminishing formulation (based on Roe's scheme) is used in the Euler solver. In the supersonic regions of the flow an 'infinitely large' time step is employed, and a finite time step is applied in the subsonic regions of the flow to reach the steady-state as a time-asymptote. Numerical solutions are obtained for a number of complex configurations, including: (1) an elliptic waverider, (2) a realistic fighter configuration, (3) the Space Shuttle, and (4) a Shuttle-like configuration. Both the Full Potential and Euler numerical results are in good agreement with available experimental data.

Szema, K. Y.↗

The acousto-ultrasonic approach

The nature and underlying rationale of the acousto-ultrasonic approach is reviewed, needed advanced signal analysis and evaluation methods suggested, and application potentials discussed. Acousto-ultrasonics is an NDE technique combining aspects of acoustic emission methodology with ultrasonic simulation of stress waves. This approach uses analysis of simulated stress waves for detecting and mapping variations of mechanical properties. Unlike most NDE, acousto-ultrasonics is less concerned with flaw detection than with the assessment of the collective effects of various flaws and material anomalies. Acousto-ultrasonics has been applied chiefly to laminated and filament-wound fiber reinforced composites. It has been used to assess the significant strength and toughness reducing effects that can be wrought by combinations of essentially minor flaws and diffuse flaw populations. Acousto-ultrasonics assesses integrated defect states and the resultant variations in properties such as tensile, shear, and flexural strengths and fracture resistance. Matrix cure state, porosity, fiber orientation, fiber volume fraction, fiber-matrix bonding, and interlaminar bond quality are underlying factors.

Vary, Alex↗

The acousto-ultrasonic approach

The nature and underlying rationale of the acousto-ultrasonic approach is reviewed, needed advanced signal analysis and evaluation methods suggested, and application potentials discussed. Acousto-ultrasonics is an NDE technique combining aspects of acoustic emission methodology with ultrasonic simulation of stress waves. This approach uses analysis of simulated stress waves for detecting and mapping variations of mechanical properties. Unlike most NDE, acousto-ultrasonics is less concerned with flaw detection than with the assessment of the collective effects of various flaws and material anomalies. Acousto-ultrasonics has been applied chiefly to laminated and filament-wound fiber reinforced composites. It has been used to assess the significant strength and toughness reducing effects that can be wrought by combinations of essentially minor flaws and diffuse flaw populations. Acousto-ultrasonics assesses integrated defect states and the resultant variations in properties such as tensile, shear, and flexural strengths and fracture resistance. Matrix cure state, porosity, fiber orientation, fiber volume fraction, fiber-matrix bonding, and interlaminar bond quality are underlying factors.

Vary, Alex↗

The looped adhesive strip: An example of coplanar delamination interaction

The phenomenon of peeling and debonding of thin layers is a subject of interest to those concerned with adhesives, thin films, and layered materials. In recent years much attention has been focused on such problems as a result of increased interest and application of advanced composites and thin film coatings. A related problem which is of interest for its own sake but also represents a simple example of a tangled adhesive strip and of coplanar delamination interaction, is the problem of a looped adhesive strip. This is the subject of the present study. Researchers consider here the problem of an elastic strip which possesses an adherend on (at least) one of its surfaces. If the strip is deformed so that two portions of such a surface are brought into contact, a position of the strip becomes bonded and a loop is formed. Researchers are interested in determining the equilibrium configuration of such a strip and investigating the behavior of the strip when its edges are pulled apart. The problem is approached as a moving interior boundary problem in the calculus of variations with the strip modeled as an inextensible elastica and the bond strength characterized by its surface energy. A Griffith type energy criterion is employed for debonding, and solutions corresponding to the problem of interest obtained. The solution obtained will be seen to predict the interesting phenomenon of bond point propagation, as well as the more standard peeling type behavior. Numerical results demonstrating the phenomena of interest are presented as well and will be seen to reveal both stable and unstable propagation of the boundaries of the bonded portion of the strip, depending upon the loading conditions.

Bottega, W. J.↗

Numerical simulations of drainage flows on Mars

Data collected by Viking Landers have shown that the meteorology of the near surface Martian environment is analogous to desertlike terrestrial conditions. Geological evidence such as dunes and frost streaks indicate that the surface wind is a potentially important factor in scouring of the martian landscape. In particular, the north polar basin shows erosional features that suggest katabatic wind convergence into broad valleys near the margin of the polar cap. The pattern of katabatic wind drainage off the north polar cap is similar to that observed on Earth over Antarctica or Greenland. The sensitivity is explored of Martian drainage flows to variations in terrain slope and diurnal heating using a numerical modeling approach. The model used is a 2-D sigma coordinate primitive equation system that has been used for simulations of Antarctic drainage flows. Prognostic equations include the flux forms of the horizontal scalar momentum equations, temperature, and continuity. Parameterization of both longwave (terrestrial) and shortwave (solar) radiation is included. Turbulent transfer of heat and momentum in the Martian atmosphere remains uncertain since relevant measurements are essentially nonexistent.

Parish, Thomas R.↗

Anomalous meridional thermospheric neutral winds in the AE-E NATE data: Effects of the equatorial nighttime pressure bulge

The work described here makes it possible to identify anomalous wind behavior such as the nighttime meridional wind abatements that occur at F-region heights. A new analysis technique uses a simple empirical wind model to simulate measurements of 'normal' winds (as measured by the Neutral Atmosphere and Temperature Experiment (NATE) that flew on the Atmosphere Explorer-E (AE-E)) to highlight anomalous wind measurements made by the satellite while in circular orbits at 270-290 km altitude. Our approach is based on the recognition that the 'in orbit' wind variation must show the combined effects of the diurnal wind variation as seen from the ground with the latitude variation of the satellite orbit. For the data period 77250-78035 examined thus far, the wind abatement always occurred with a corresponding pressure or temperature maximum, and was detected on 12 out of the 36 nights with data. This study has revealed that the wind abatement occur only during or shortly after increases in solar extreme ultraviolet (EUV) flux, as indicated by daily radio flux measurements. In the past, nighttime wind reversals at mid-latitudes have been associated with increased geomagnetic activity. This study indicates that intensified solar EUV heating may be responsible for anomalous thermospheric nighttime winds at mid-latitudes.

Goembel, L.↗

A 2-D Interface Element for Coupled Analysis of Independently Modeled 3-D Finite Element Subdomains

Over the past few years, the development of the interface technology has provided an analysis framework for embedding detailed finite element models within finite element models which are less refined. This development has enabled the use of cascading substructure domains without the constraint of coincident nodes along substructure boundaries. The approach used for the interface element is based on an alternate variational principle often used in deriving hybrid finite elements. The resulting system of equations exhibits a high degree of sparsity but gives rise to a non-positive definite system which causes difficulties with many of the equation solvers in general-purpose finite element codes. Hence the global system of equations is generally solved using, a decomposition procedure with pivoting. The research reported to-date for the interface element includes the one-dimensional line interface element and two-dimensional surface interface element. Several large-scale simulations, including geometrically nonlinear problems, have been reported using the one-dimensional interface element technology; however, only limited applications are available for the surface interface element. In the applications reported to-date, the geometry of the interfaced domains exactly match each other even though the spatial discretization within each domain may be different. As such, the spatial modeling of each domain, the interface elements and the assembled system is still laborious. The present research is focused on developing a rapid modeling procedure based on a parametric interface representation of independently defined subdomains which are also independently discretized.

Kandil, Osama A.↗

Determination of Creep Behavior of Thermal Barrier Coatings Under Laser Imposed High Thermal and Stress Gradient Conditions

A laser sintering/creep technique has been established to determine the creep behavior of thermal barrier coatings under steady-state high heat flux conditions. For a plasma sprayed zirconia-8 wt. % yttria coating, a significant primary creep strain and a low apparent creep activation energy were observed. Possible creep mechanisms involved include stress induced mechanical sliding and temperature and stress enhanced cation diffusion through the splat and grain boundaries. The elastic modulus evolution, stress response, and total accumulated creep strain variation across the ceramic coating are simulated using a finite difference approach. The modeled creep response is consistent with experimental observations.

Zhu, Dongming↗