Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical continuation”

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 1,441 records · Page 80
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

A time-parallel multiple-shooting method for large-scale quantum optimal control

Quantum optimal control plays a crucial role in quantum computing by providing the interface between compiler and hardware. Solving the optimal control problem is particularly challenging for multi-qubit gates, due to the exponential growth in computational complexity with the system's dimensionality and the deterioration of optimization convergence. To ameliorate the computational complexity of time-integration, this paper introduces a multiple-shooting approach in which the time domain is divided into multiple windows and the intermediate states at window boundaries are treated as additional optimization variables. Further, this enables parallel computation of state evolution across time-windows, significantly accelerating objective function and gradient evaluations. Since the initial state matrix in each window is only guaranteed to be unitary upon convergence of the optimization algorithm, the conventional gate trace infidelity is replaced by a generalized infidelity that is convex for non-unitary state matrices. Continuity of the state across window boundaries is enforced by equality constraints. A quadratic penalty optimization method is used to solve the constrained optimal control problem, and an efficient adjoint technique is employed to calculate the gradients in each iteration. We demonstrate the effectiveness of the proposed method through numerical experiments on quantum Fourier transform gates in systems with 2, 3, and 4 qubits, noting a speedup of 80x for evaluating the gradient in the 4-qubit case, highlighting the method's potential for optimizing control pulses in multi-qubit quantum systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Versatile tape-drive target for high-repetition-rate laser-driven proton acceleration

We present the development and characterization of a high-stability, multi-material, multi-thickness tape-drive target for laser-driven acceleration at repetition rates of up to 100 Hz. The tape surface position was measured to be stable on the sub-micrometre scale, compatible with the high-numerical aperture focusing geometries required to achieve relativistic intensity interactions with the pulse energy available in current multi-Hz and near-future higher repetition-rate lasers ( > kHz). Long-term drift was characterized at 100 Hz demonstrating suitability for operation over extended periods. The target was continuously operated at up to 5 Hz in a recent experiment for 70,000 shots without intervention by the experimental team, with the exception of tape replacement, producing the largest data-set of relativistically intense laser–solid foil measurements to date. This tape drive provides robust targetry for the generation and study of high-repetition-rate ion beams using next-generation high-power laser systems, also enabling wider applications of laser-driven proton sources.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

Analytical desmearing of Bonse–Hart ultra-small-angle neutron scattering data via truncated Abel inversion

A non-iterative analytical framework based on the truncated Abel inversion is developed for desmearing Bonse–Hart ultra-small-angle neutron scattering (USANS) data. The method directly inverts the slit-averaged intensity without empirical extrapolation or iterative regularization, establishing a closed-form relationship between the measured and intrinsic scattering profiles. Numerical benchmarks on representative models, including a rigid-line form factor, a Lorentzian function and a fractal structural model, demonstrate quantitative recovery of the ground-truth intensity across the full Q range. Application to a deuterated polystyrene/poly(2-vinylpyridine) blend further confirms that the approach yields smooth continuous profiles consistent with companion small-angle neutron scattering data. The truncated Abel inversion thus provides a stable, model-independent and physically transparent route for accurate desmearing of Bonse–Hart USANS measurements.

Huang, Guan-Rong [National Tsing Hua University, T↗
DOE OSTI2020

FY20 SAM Code Developments and Validations for Transient Safety Analysis of Advanced non-LWRs

The System Analysis Module (SAM) is under development at Argonne National Laboratory as a modern system-level modeling and simulation tool for advanced non-light water reactor safety analyses. It utilizes the object-oriented application framework MOOSE to leverage the modern software environment and advanced numerical methods. The capabilities of SAM are being extended to enable the transient modeling, analysis, and design of various advanced nuclear reactor systems. This report summarizes major progress in SAM code development, capability enhancements, demonstration, and validation to support transient safety analysis of advanced non-LWRs. Rapid developments continued in fiscal year 2020 (FY20) to support various needs of the advanced reactor community, especially the NRC and industry on the licensing safety analysis of advanced reactor designs. Significant code changes were made to provide various capability enhancements, bug fixes, and user friendliness improvements. Major code updates are summarized in Section 1, while four important enhancements are detailed in Sections 2-5, including a multi-dimension flow model; reactivity feedback and decay heat models; control and trip system modeling, and additional fluid and solid thermophysical property models. Code validation activities in FY20 include using test data from the Fast Flux Test Facility (FFTF), the High Temperature Test Facility (HTTF), and several separate effects test facilities for pebble-bed modeling.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS1972

Regional variations in degradation and density of Martian craters.

Martian craters visible on Mariner 6 and 7 imagery show a spectrum of topographic types from very fresh to highly degraded. A method of numerical scoring of rim, wall, and floor is proposed to yield a degradation number to classify each crater. Plots of degradation class vs density of large craters are similar for all four regions studied, whereas similar plots for small craters show marked differences between regions. The data suggest general continuity of crater formation and degradation, along with locally sporadic formation and/or degradation of the smallest craters classified. Deucalionis Regio, with an excess of fresh, small craters, experienced an episode of small crater formation (or nondegradation) most recently; Margaritifer Sinus was similarly disturbed at some more remote time. Meridiani Sinus and Hellespontus-Noachis show little or no sign of excess fresh, small craters.

Mcgill, G. E.↗
NASA NTRS1974

Research study on stabilization and control modern sampled-data control theory

The methods of continuous and discrete describing function analysis were applied to predicting the existence of self-sustained oscillations in the single-axis model of the large space telescope system with nonlinear control moment gyroscope friction characteristics. It is shown that the stability equations may be solved by a numerical-iterative technique using the describing function analysis, instead of the usual graphical methods. The numerical method is found to be effective in leading to a convergent solution rapidly, with an appropriate guess of the initial condition.

Kuo, B. C.↗
NASA NTRS1950

Linear theory of boundary effects in open wind tunnels with finite jet lengths

In the first part, the boundary conditions for an open wind tunnel (incompressible flow) are examined with special reference to the effects of the closed entrance and exit sections. Basic conditions are that the velocity must be continuous at the entrance lip and that the velocities in the upstream and downstream closed portions must be equal. In the second part, solutions are derived for four types of two-dimensional open tunnels, including one in which the pressures on the two free surfaces are not equal. Numerical results are given for every case. In general, if the lifting element is more than half the tunnel height from the inlet, the boundary effect at the lifting element is the same as for an infinitely long open tunnel. In the third part, a general method is given for calculating the boundary effect in an open circular wind tunnel of finite jet length. Numerical results are given for a lifting element concentrate at a point on the axis.

Katzoff, S↗
NASA NTRS1994

Direct numerical simulation of incompressible axisymmetric flows

In the present work, we propose to conduct direct numerical simulations (DNS) of incompressible turbulent axisymmetric jets and wakes. The objectives of the study are to understand the fundamental behavior of axisymmetric jets and wakes, which are perhaps the most technologically relevant free shear flows (e.g. combuster injectors, propulsion jet). Among the data to be generated are various statistical quantities of importance in turbulence modeling, like the mean velocity, turbulent stresses, and all the terms in the Reynolds-stress balance equations. In addition, we will be interested in the evolution of large-scale structures that are common in free shear flow. The axisymmetric jet or wake is also a good problem in which to try the newly developed b-spline numerical method. Using b-splines as interpolating functions in the non-periodic direction offers many advantages. B-splines have local support, which leads to sparse matrices that can be efficiently stored and solved. Also, they offer spectral-like accuracy that are C(exp O-1) continuous, where O is the order of the spline used; this means that derivatives of the velocity such as the vorticity are smoothly and accurately represented. For purposes of validation against existing results, the present code will also be able to simulate internal flows (ones that require a no-slip boundary condition). Implementation of no-slip boundary condition is trivial in the context of the b-splines.

Loulou, Patrick↗
NASA NTRS2001

Rheology of Foam Near the Order-Disorder Transition

The first part of our research results are summarized in the recent journal publication: J. Gregory McDaniel and R. Glynn Holt, 'Measurement of aqueous foam rheology by acoustic levitation', Phys. Rev. E 61, 2204 (2000). This aspect of the work was a combination of experiment and analysis. We built a levitation system capable of acoustically levitating small samples of aqueous foam of arbitrary gas and liquid volume fractions. We then modulated the acoustic field to induce normal mode oscillations of the foam samples. The observables from the experiment were frequency and mode number. For dry (roughly > 70% gas by volume) foams and small deformations, we developed an effective medium, normal-modes analysis which took the frequency and mode number from experiment, and gave us the shear elastic modulus of the foam as a function of Poisson's ratio. The second part of our results may be found in a soon-to-be submitted manuscript 'Dynamics of aqueous foam drops', I.Sh. Akhatov, J.G. McDaniel and R.G. Holt, describing our modeling in the wet foam limit by considering the acoustic problem. This aspect of the research is purely theoretical. Beginning from a mass-conserving mixture law, the fully nonlinear equations of motion for a wet (roughly < 10% gas by volume) foam drop of initially spherical shape were derived. The frequencies for normal mode oscillations were derived in the linear inviscid limit. The nonlinear equations were numerically solved to elicit the motion of a foam drop under acoustic excitation. The role of the time-varying void fraction in breathing-mode oscillations is of particular interest. As of the end of the current (NAG#3-2121) grant, this work was not yet concluded. We continue to work on this aspect in order to extend the analysis to cover the transition regime of gas volume fractions, as well as to compare to experiments in the wet regime.

Holt, R. Glynn↗
NASA NTRS2003

Detailed Vibration Analysis of Pinion Gear with Time-Frequency Methods

In this paper, the authors show a detailed analysis of the vibration signal from the destructive testing of a spiral bevel gear and pinion pair containing seeded faults. The vibration signal is analyzed in the time domain, frequency domain and with four time-frequency transforms: the Short Time Frequency Transform (STFT), the Wigner-Ville Distribution with the Choi-Williams kernel (WV-CW), the Continuous Wavelet' Transform (CWT) and the Discrete Wavelet Transform (DWT). Vibration data of bevel gear tooth fatigue cracks, under a variety of operating load levels and damage conditions, are analyzed using these methods. A new metric for automatic anomaly detection is developed and can be produced from any systematic numerical representation of the vibration signals. This new metric reveals indications of gear damage with all of the time-frequency transforms, as well as time and frequency representations, on this data set. Analysis with the CWT detects changes in the signal at low torque levels not found with the other transforms. The WV-CW and CWT use considerably more resources than the STFT and the DWT. More testing of the new metric is needed to determine its value for automatic anomaly detection and to develop fault detection methods for the metric.

Mosher, Marianne↗
NASA NTRS2016

Overview of the CHarring Ablator Response (CHAR) Code

An overview of the capabilities of the CHarring Ablator Response (CHAR) code is presented. CHAR is a one-, two-, and three-dimensional unstructured continuous Galerkin finite-element heat conduction and ablation solver with both direct and inverse modes. Additionally, CHAR includes a coupled linear thermoelastic solver for determination of internal stresses induced from the temperature field and surface loading. Background on the development process, governing equations, material models, discretization techniques, and numerical methods is provided. Special focus is put on the available boundary conditions including thermochemical ablation, surface-to-surface radiation exchange, and flowfield coupling. Finally, a discussion of ongoing development efforts is presented.

Amar, Adam J.↗
NASA NTRS2016

Overview of the CHarring Ablator Response (CHAR) Code

An overview of the capabilities of the CHarring Ablator Response (CHAR) code is presented. CHAR is a one-, two-, and three-dimensional unstructured continuous Galerkin finite-element heat conduction and ablation solver with both direct and inverse modes. Additionally, CHAR includes a coupled linear thermoelastic solver for determination of internal stresses induced from the temperature field and surface loading. Background on the development process, governing equations, material models, discretization techniques, and numerical methods is provided. Special focus is put on the available boundary conditions including thermochemical ablation and contact interfaces, and example simulations are included. Finally, a discussion of ongoing development efforts is presented.

Amar, Adam J.↗
NASA NTRSDate not supplied

The Challenger tragedy was caused by an Apollo mistake, terminating risk analysis

NASA’s view of risk changed between early Apollo and the Space Shuttle. Risk was a known serious problem at the beginning of Apollo and the risk estimates were disturbingly high. To avoid public concern, risk analysis was discontinued. Risk analysis was avoided in Shuttle, leading to an unnecessarily risky design. The immediate cause of the Challenger tragedy was the mistaken decision to launch in cold weather. The fundamental cause was the high risk of the Shuttle design. Before Challenger, management thought and testified that the probability of an accident was 1 in 100,000. After Challenger, Probabilistic Risk Analysis (PRA) found a roughly 1 in 100 chance of a Shuttle failure. The recent Orion design uses the safer Apollo approach, with a hardened capsule, launch abort escape, and the crew placed above the rocket tanks and engines. During Apollo it was estimated that, “assuming all elements from propulsion to rendezvous and life support were done as well or better than ever before, that 30 astronauts would be lost before 3 were returned safely to the Earth.” The chance of astronaut survival was only 10%. After the Apollo 1 tragedy, the awareness of risk led to an intense focus on achieving safety. “The only possible explanation for the astonishing success – no losses in space and on time – was that every participant at every level in every area far exceeded the norm of human capabilities.” During Apollo, a NASA PRA found that the chance of success was “less than 5 percent.” The NASA Administrator felt that “the numbers could do irreparable harm,” and discontinued numerical risk assessment. This led to decreasing understanding of risk. The head of Apollo reliability and safety decided, “Statistics don’t count for anything,” and that risk is reduced by “attention taken in design.” The great and initially unexpected success of Apollo appeared to validate the neglect of PRA. Continuing to neglect the mathematical estimation of risk led Shuttle into a high risk design that produced tragic results. The initial design of the Shuttle emphasized increasing capability and reducing cost without analysis or even mention of risk. A retired NASA official stated, “some NASA people began to confuse desire with reality. … One result was to assess risk in terms of what was thought acceptable without regard for verifying the assessment. … Note that under such circumstances real risk management is shut out.” Not computing risk led to removing launch abort, removing crew escape, selecting less reliable Solid Rocket Boosters, placing the crew compartment next to the rocket boosters, and accepting more stressed shielding tile designs. Accepting these specific risks directly caused the shuttle disasters. The Challenger tragedy is frequently taught as a case of management failure. The focus is on the Challenger launch decision hours before, which is a dramatic example of bad management. However, the true cause of the Challenger disaster occurred decades earlier in the Apollo era. When the easily predictable failures occurred, failure investigations focused on how they might have been avoided. The Shuttle was cancelled after the space station was completed because of its high risk. The ultimate cause of the Shuttle tragedies was the choice by the Apollo-era NASA administrator to avoid a negative public reaction to realistic risk analysis.

Harry W. Jones↗
DOE OSTI2021

Closed-Form Solutions for the Equations of Motion of the Heavy Symmetrical Top with One Point Fixed

The equations of motion (EOM) for the heavy symmetrical top with one point fixed are highly nonlinear. The literature describes the numerical methods that are used to resolve this classical system, including modern tools, such as the Runge–Kutta fourth–order method. Finding the derivate of closed-form solutions for the EOM is more difficult and, as mentioned in the literature, discovering the solution is not always possible for all the EOM. Fortunately, a few examples are available that serve as a guide to move further in this topic. The purpose of this paper is to find a methodology that will produce the solutions for a given subset of EOMs that fulfill certain requisites. This paper summarizes the literature available on this topic and then follows with the derivation of the EOM using the Euler–Lagrange method. The Routhian method will be used to reduce the size of the expression, and it continues with the formulation of the classical cubic function, ƒ(u), through a novel process. The roots of ƒ(u) are of the utmost importance in finding the EOM closed-form solution, and once the final roots are selected, the general method that will produce the closed-form solutions is presented. Two sets of examples are included to show the validity of the process, and comparisons of the results from the closed-form solutions vs. the numerical results for these examples are shown.

Laos, Hector↗
NASA NTRS1994

Testing the gravitational instability hypothesis?

We challenge a widely accepted assumption of observational cosmology: that successful reconstruction of observed galaxy density fields from measured galaxy velocity fields (or vice versa), using the methods of gravitational instability theory, implies that the observed large-scale structures and large-scale flows were produced by the action of gravity. This assumption is false, in that there exist nongravitational theories that pass the reconstruction tests and gravitational theories with certain forms of biased galaxy formation that fail them. Gravitational instability theory predicts specific correlations between large-scale velocity and mass density fields, but the same correlations arise in any model where (a) structures in the galaxy distribution grow from homogeneous initial conditions in a way that satisfies the continuity equation, and (b) the present-day velocity field is irrotational and proportional to the time-averaged velocity field. We demonstrate these assertions using analytical arguments and N-body simulations. If large-scale structure is formed by gravitational instability, then the ratio of the galaxy density contrast to the divergence of the velocity field yields an estimate of the density parameter Omega (or, more generally, an estimate of beta identically equal to Omega(exp 0.6)/b, where b is an assumed constant of proportionality between galaxy and mass density fluctuations. In nongravitational scenarios, the values of Omega or beta estimated in this way may fail to represent the true cosmological values. However, even if nongravitational forces initiate and shape the growth of structure, gravitationally induced accelerations can dominate the velocity field at late times, long after the action of any nongravitational impulses. The estimated beta approaches the true value in such cases, and in our numerical simulations the estimated beta values are reasonably accurate for both gravitational and nongravitational models. Reconstruction tests that show correlations between galaxy density and velocity fields can rule out some physically interesting models of large-scale structure. In particular, successful reconstructions constrain the nature of any bias between the galaxy and mass distributions, since processes that modulate the efficiency of galaxy formation on large scales in a way that violates the continuity equation also produce a mismatch between the observed galaxy density and the density inferred from the peculiar velocity field. We obtain successful reconstructions for a gravitational model with peaks biasing, but we also show examples of gravitational and nongravitational models that fail reconstruction tests because of more complicated modulations of galaxy formation.

Babul, Arif↗
NASA NTRS2011

Geostationary Microwave Sounders: Science, Applications and the Geostar Instrument Concept

Microwave atmospheric sounders have long provided some of the most imporant data for use in numerical weather prediction (NWP) and have played an important role in atmospheric weather and climate research. With 7 US satellites now carrying such sensors, we are in a 'golden age' of microwave remote sensing of the atmosphere. However, as this fleet ages and is replaced by a smaller number of new sensors in the coming yars, the main shortcoming of sensors in low Earth orbit -i.e. poor spacial and temporal converage and sampling - will become more apparent. Placing such sensors on geostationary satellites, enabling time-continuous views of large portions of the Earth disc, would solve this problem. but the GEO orbit is approximately 40 times higher than a typical LEO orbit, which requires antenna apertures also about 40 times larger than for LEO systems to maintain spatial resolution, and it has not been feasible to develop such systems. Recently, a solution to this problem has appeared in the form of aperture synthesis.

numerical weather prediction (NWP)↗
NASA NTRS2008

Numerically Modeling the Erosion of Lunar Soil by Rocket Exhaust Plumes

In preparation for the Apollo program, Leonard Roberts of the NASA Langley Research Center developed a remarkable analytical theory that predicts the blowing of lunar soil and dust beneath a rocket exhaust plume. Roberts assumed that the erosion rate was determined by the excess shear stress in the gas (the amount of shear stress greater than what causes grains to roll). The acceleration of particles to their final velocity in the gas consumes a portion of the shear stress. The erosion rate continues to increase until the excess shear stress is exactly consumed, thus determining the erosion rate. Roberts calculated the largest and smallest particles that could be eroded based on forces at the particle scale, but the erosion rate equation assumed that only one particle size existed in the soil. He assumed that particle ejection angles were determined entirely by the shape of the terrain, which acts like a ballistic ramp, with the particle aerodynamics being negligible. The predicted erosion rate and the upper limit of particle size appeared to be within an order of magnitude of small-scale terrestrial experiments but could not be tested more quantitatively at the time. The lower limit of particle size and the predictions of ejection angle were not tested. We observed in the Apollo landing videos that the ejection angles of particles streaming out from individual craters were time-varying and correlated to the Lunar Module thrust, thus implying that particle aerodynamics dominate. We modified Roberts theory in two ways. First, we used ad hoc the ejection angles measured in the Apollo landing videos, in lieu of developing a more sophisticated method. Second, we integrated Roberts equations over the lunar-particle size distribution and obtained a compact expression that could be implemented in a numerical code. We also added a material damage model that predicts the number and size of divots which the impinging particles will cause in hardware surrounding the landing rocket. Then, we performed a long-range ballistics analysis for the ejected particulates.

Source record↗