SEARCH · Engineering Papers
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.
Based on variances of log-amplitude and phase fluctuations of radio occultation data received from orbital and fly-by missions, structure constant for Venusian planetary atmosphere has been estimated with high-confidence factor. Analysis indicates that effects of inhomogeneity, finite size, and superrefractivity of atmospheric turbulence cannot be ignored.
An investigation of delamination crack growth in a unidirectional composite laminate under static and fatigue loading conditions is presented. An experimental program was conducted to study failure mechanisms and characterize the crack growth rate; analytical modeling based on microscopic observations was carried out using an advanced singular finite element method of the hybrid-stress model. Using a mixed-mode failure criterion, the delamination growth under a monotonically increasing load can be predicted; the rate of delamination crack propagation under fatigue loading is directly related to the amplitudes of cyclic opening and shearing mode stress intensity factors by a power law relationship.
A numerical model has been developed for investigating boundary layer transition control for a flat plate boundary layer. Active control of a periodically forced boundary layer in an incompressible fluid is studied using surface heating techniques. The spatially evolving boundary layer is simulated. The Navier-Stokes and energy equations are integrated using a fully implicit finite difference/spectral method. Temperature perturbations are introduced locally along finite heater strips to directly attenuate the instability waves in the flow. A feedback control loop is employed in which a downstream sensor is used to monitor wall shear stress fluctuations. Active control of small amplitude two-dimensional and three-dimensional disturbances is numerically simulated. With proper phase control, in-phase reinforcement and out-of-phase attenuation are demonstrated. A receptivity study of the localized temperature perturbations is made. It is shown that narrow heater strips are more receptive in that they maximize the amplitude level of the disturbances in the flow. Active control of the early stages of the fundamental breakdown process is also numerically simulated. Control is achieved with either two-dimensional or three-dimensional control inputs.
ABSTRACT There is a growing set of observational data demonstrating that cosmic rays exhibit small-scale anisotropies (5°–30°) with amplitude deviations lying between 0.01–0.1 per cent that of the average cosmic ray flux. A broad range of models have been proposed to explain these anisotropies ranging from finite-scale magnetic field structures to dark matter annihilation. The standard diffusion transport methods used in cosmic ray propagation do not capture the transport physics in a medium with finite-scale or coherent magnetic field structures. Here, we present a Monte Carlo transport method, applying it to a series of finite-scale magnetic field structures to determine the requirements of such fields in explaining the observed cosmic ray, small-scale anisotropies.
A finite difference method for the solution of the transonic flow about a harmonically oscillating wing is presented. The partial differential equation for the unsteady transonic flow was linearized by dividing the flow into separate steady and unsteady perturbation velocity potentials and by assuming small amplitudes of harmonic oscillation. The resulting linear differential equation is of mixed type, being elliptic or hyperbolic whereever the steady flow equation is elliptic or hyperbolic. Central differences were used for all derivatives except at supersonic points where backward differencing was used for the streamwise direction. Detailed formulas and procedures are described in sufficient detail for programming on high speed computers. To test the method, the problem of the oscillating flap on a NACA 64A006 airfoil was programmed. The numerical procedure was found to be stable and convergent even in regions of local supersonic flow with shocks.
Here we developed and implemented a method-independent, fully numerical, finite difference approach to calculating nuclear magnetic resonance shieldings, using gauge-including atomic orbitals. The resulting capability can be used to explore non-standard methods, given only the energy as a function of finite-applied magnetic fields and nuclear spins. For example, standard second-order Møller-Plesset theory (MP2) has well-known efficacy for 1H and 13C shieldings and known limitations for other nuclei such as 15N and 17O. It is, therefore, interesting to seek methods that offer good accuracy for 15N and 17O shieldings without greatly increased compute costs, as well as exploring whether such methods can further improve 1H and 13C shieldings. Using a small molecule test set of 28 species, we assessed two alternatives: κ regularized MP2 (κ-MP2), which provides energy-dependent damping of large amplitudes, and MP2.X, which includes a variable fraction, X, of third-order correlation (MP3). The aug-cc-pVTZ basis was used, and coupled cluster with singles and doubles and perturbative triples [CCSD(T)] results were taken as reference values. Our κ-MP2 results reveal significant improvements over MP2 for 13C and 15N, with the optimal κ value being element-specific. κ-MP2 with κ = 2 offers a 30% rms error reduction over MP2. For 15N, κ-MP2 with κ = 1.1 provides a 90% error reduction vs MP2 and a 60% error reduction vs CCSD. On the other hand, MP2.X with a scaling factor of 0.6 outperformed CCSD for all heavy nuclei. These results can be understood as providing renormalization of doubles amplitudes to partially account for neglected triple and higher substitutions and offer promising opportunities for future applications.
The roles of natural convection in the melt and the shape of the melt/solid interface on radial dopant segregation are analyzed for a prototype of vertical Bridgman crystal growth system by finite element methods that solve simultaneously for the velocity field in the melt, the shape of the solidification isotherm, and the temperature distribution in both phases. Results are presented for crystal and melt with thermophysical properties similar to those of gallium-doped germanium in Bridgman configurations with melt below (thermally destabilizing) and above (stabilizing) the crystal. Steady axisymmetric flow are classified according to Rayleigh number as either being nearly the growth velocity, having a weak cellular structure or having large amplitude cellular convention. The flows in the two Bridgman configurations are driven by different temperature gradients and are in opposite directions. Finite element calculations for the transport of a dilute dopant by these flow fields reveal radial segregation levels as large as sixty percent of the mean concentration. Segregation is found most severe at an intermediate value of Rayleigh number above which the dopant distribution along the interface levels as the intensity of the flow increases.
A theory of the measurement of the thermal diffusivity of a sample by the modified Angstrom method is developed for the case in which radiative heat loss from the end surface of the sample is not negligible, and applied to measurements performed on lunar samples. Formulas allowing sample thermal diffusivity to be determined from the amplitude decay and phase lag of a temperature wave traveling through the sample are derived for a flat disk sample for which only heat loss from the end surface is important, and a sample of finite diameter and length for which heat loss through the end and side surfaces must be considered. It is noted that in the case of a flat disk, measurements at a single angular frequency of the temperature wave are sufficient, while the sample of finite diameter and length requires measurements at two discrete angular frequencies. Comparison of the values of the thermal diffusivities of two lunar samples of dimensions approximately 1 x 1 x 2 cm derived by the present methods and by the Angstrom theory for a finite bar reveals them to differ by not more than 5%, and indicates that more refined data are required as the measurement theory becomes more complicated.
A program for predicting the linear stability of liquid propellant rocket engines is presented. The underlying model assumptions and analytical steps necessary for understanding the program and its input and output are also given. The rocket engine is modeled as a right circular cylinder with an injector with a concentrated combustion zone, a nozzle, finite mean flow, and an acoustic admittance, or the sensitive time lag theory. The resulting partial differential equations are combined into two governing integral equations by the use of the Green's function method. These equations are solved using a successive approximation technique for the small amplitude (linear) case. The computational method used as well as the various user options available are discussed. Finally, a flow diagram, sample input and output for a typical application and a complete program listing for program MODULE are presented.
The impingement of a row of finite-area vortices on an edge is presently used to efficiently simulate the shear layer/edge interaction, yielding support for the hypothesis that the pressure waves emitted from an impingement edge are generated by the vortices/edge interaction. A parametric application of this method shows that pressure wave amplitude is a function of the length of the succession of vortices and that frequency of their release; this amplitude decreases with decreasing vortex spacing while succession length remains constant, or when succession length decreases while the number of vortices remains constant.
The non-equilibrium critical-layer analysis of a system of frequency-detuned resonant-triads is presented. In this part of the analysis, the system of partial differential critical-layer equations derived in Part I is solved analytically to yield the amplitude equations which are analyzed using a combination of asymptotic and numerical methods. Numerical solutions of the inviscid non-equilibrium oblique-mode amplitude equations show that the frequency-detuned self-interaction enhances the growth of the lower-frequency oblique modes more than the higher-frequency ones. All amplitudes become singular at the same finite downstream position. The frequency detuning delays the occurrence of the singularity. The spanwise-periodic mean-flow distortion and low-frequency nonlinear modes are generated by the critical-layer interaction between frequency-detuned oblique modes. The nonlinear mean flow and higher harmonics as well as the primary instabilities become as large as the base mean flow in the inviscid wall layer in the downstream region where the distance from the singularity is of the order of the wavelength scale.
Quantum simulation of quantum field theory is a flagship application of quantum computers that promises to deliver capabilities beyond classical computing. The realization of quantum advantage will require methods that can accurately predict error scaling as a function of the resolution and parameters of the model and that can be implemented efficiently on quantum hardware. In this paper, we address the representation of lattice bosonic fields in a discretized field amplitude basis, develop methods to predict error scaling, and present efficient qubit implementation strategies. A low-energy subspace of the bosonic Hilbert space, defined by a boson occupation number cutoff, can be represented with exponentially good accuracy by a low-energy subspace of a finite-size Hilbert space. The finite representation construction and the associated errors are directly related to the accuracy of the Nyquist-Shannon sampling and the finite Fourier transforms of the boson number states in the field and the conjugate-field bases. We analyze the relation between the boson mass, the discretization parameters used for wave function sampling, and the finite representation size. Numerical simulations of small size Φ 4 problems demonstrate that the boson mass optimizing the sampling of the ground state wave function is a good approximation to the optimal boson mass yielding the minimum low-energy subspace size. However, we find that accurate sampling of general wave functions does not necessarily result in accurate representation. Finally, we develop methods for validating and adjusting the discretization parameters to achieve more accurate simulations.
The general one-dimensional non-linear equation for the acoustic velocity potential in a variable area duct carrying high subsonic Mach number flows is presented and solved numerically using an implicit finite difference scheme. For the linearized equation with no flow, the present scheme is compared with the exact solution and the fourth-order Runge-Kutta method with excellent agreement for dimensionless time periods. Non-linear solutions are more sensitive to Mach numbers and exciting amplitudes and less sensitive to exciting frequencies. In general, non-linear effects can be safely neglected for low Mach number flows but must be accounted for when high Mach number flows are encountered.
Results obtained from a numerical simulation of the final stages of transition to turbulence in plane channel flow are described. Three dimensional, incompressible Navier-Stokes equations are numerically integrated to obtain the time evolution of two and three dimensional finite amplitude disturbances. Computations are performed on CYBER-203 vector processor for a 32x51x32 grid. Results are presented for no-slip boundary conditions at the solid walls as well as for periodic suction blowing to simulate active control of transition by mass transfer. Solutions indicate that the method is capable of simulating the complex character of vorticity dynamics during the various stages of transition and final breakdown. In particular, evidence points to the formation of a lambda-shape vortex and the subsequent system of horseshoe vortices inclined to the main flow direction as the main elements of transition. Calculations involving periodic suction-blowing indicate that interference with a wave of suitable phase and amplitude reduces the disturbance growth rates.
Cyclic loading or other stresses can lead to development of cracks and crack growth in mechanical structures, leading to eventual failure. While ultrasound imaging can be used for non-destructive testing of such structures, conventional ultrasound techniques are often limited by crack size, density, and areal coverage. An effective characterization of real-world, large-area structures is required at an early damage stage to prevent catastrophic failure and predict remaining life. In this study, a new nonlinear ultrasonic testing (NUT) method is proposed for large-area monitoring of practical structures with arbitrary complexity by using multiple-mode guided-wave ultrasonic signals. The proposed guided-wave NUT technique requires single-element transducers, simple electronics, and a mixed time-frequency domain signal processing. As a proof-of-concept demonstration, numerical simulations and experiments are performed on an A36 carbon steel beam assembly with previously formed microstructural defects that cause nonlinearities in ultrasonic response. The quadratic dependence of the nonlinear wave excitation on the input ultrasonic signal amplitude is shown by numerical simulations, and such a nonlinear ultrasonic response is experimentally observed in the zone with a high density of microstructural defects.
Microscale, pick-and-place assembly is a non-lithographic assembly method poised to impact diverse fields including flexible electronics, microfluidics and robotics. However, a major technological challenge is the need to deterministically control adhesion between parts. Here, switchable adhesion involving 3D-printed, self-complementary surfaces is demonstrated. Mechanical properties of metasurfaces pressed against flat, rigid substrates are modeled using finite element methods. A series of flat slabs and metastructured slabs with 2D sinusoidal surfaces are printed using two-photon polymerization (2PP) of a shape-memory resin. The surface frequency of featured slabs was varied between $3.\bar3$ mm −1 and $26.\bar6$ mm −1 with similar amplitudes. Adhesion between printed metasurfaces and glass and between printed, self-complementary metasurfaces is studied above and below the cured resin’s glass transition temperature (∼45 °C). Simple heating of adhering surfaces to above 60 °C lowers adhesion, and compression of surfaces while above the glass transition temperature followed by cooling to room temperature elevates adhesion. The nominal adhesive strength between printed, self-complementary surfaces, as determined by the maximum observable pull-off stress, exceeds 3 MPa. Further tailoring complementary surfaces for adhesion control may facilitate microscale disassembly for recovery of components or precious metals.
Noise generated in aircraft engines is usually suppressed by acoustically treating the engine ducts. The optimization of this treatment requires an understanding of the transmission and attenuation of the acoustic waves. A critical review is presented of the state of the art regarding methods of determining the transmission and attenuation parameters and the effect on these parameters of (1) acoustic properties of liners, (2) the mean velocity, including uniform and shear profiles and nonparallel flow, (3) axial and transverse temperature gradients, (4) slowly and abruptly varying cross sections, and (5) finite-amplitude waves and nonlinear duct liners.
The decomposition of the four-point ABJM amplituhedron into negative geometries produces compact integrands of logarithmic of amplitudes such that the infrared divergence only comes from the last loop integration, from which we can compute the cusp anomalous dimension of the ABJM theory. In this note, we integrate L – 1 loop momenta of the L-loop negative geometries for all four-loop negative geometries and a special class of all-loop ladder-type negative geometries by a method based on Mellin transformation, and from these finite quantities we extract the corresponding contribution to the cusp anomalous dimension. We find that the infrared divergence of a box-type negative geometry at L = 4 is weaker than other negative geometries, then only tree-type negative geometries contribute to the cusp anomalous dimension at L = 4. For the all-loop ladder-type negative geometries, we prove and conjecture some recursive structures as integral equations in Mellin space and find that they cannot contribute zeta values like ζ 3 , ζ 5 to the cusp anomalous dimension.