Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nonlinear processes”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Characterization of nonlinear ultrasonic waves behavior while interacting with poor interlayer bonds in large-scale additive manufactured materials

Over the past decades, researchers have developed several nonlinear ultrasonic techniques for quality control of materials commonly used in different applications. Owing to the superior sensitivity of nonlinear ultrasound waves to small defects such as micro-cracks, their applicability in different nondestructive testing (NDT) problems has been investigated in numerous studies. These studies utilize frequency domain analysis to detect the generation of higher harmonics because of the formation of defects in the inspected medium. Frequency domain analysis based on the Fourier transform is a significant approach used in linear systems; however, it may not perform adequately on nonlinear systems. Hence, studies on nonlinear dynamics and physics consider analyzing systems' behavior in the phase-space domain. In contrast to the frequency domain analysis, which can result in information loss, analysis in the phase-space domain retains all the information regarding a system's states. Here, we investigate the nonlinearities induced by poor interlayer bonding in polymer-based additive manufactured parts in the phase-space domain. It is convenient to characterize the nonlinearity in the phase-space domain because it provides a geometrical representation of a system's states. Two types of low quality interlayer bond are considered. The first type is simulated artificially while the second type is manufactured by reducing the bond quality during the printing process. The analysis verified that the received ultrasonic signals exhibit classical nonlinear behavior in the phase-space domain while interacting with simulated poor interlayer bonds. In addition, the results showed that the behavior of ultrasonic waves is amplitude-dependent and evolves into models that have not been previously reported. Furthermore, Largest Lyapunov Exponent (LLE) is used to quantify the behavior of nonlinear ultrasonic waves while interacting with poor interlayer bonds. Using LLE, it was observed that the divergence rate of the phase-space trajectories depends on the amplitude of the excitation. This observation quantitatively proves that nonlinear behavior of ultrasound while interacting with poor interlayer bonds can be amplitude-dependent. The results of both simulated and inherent poor interlayer bond cases showed that LLE can be used as a reliable quantitative damage-sensitive feature to detect and potentially characterize weak bonds, which are difficult to detect using conventional approaches. Additionally, the reported results in the phase-space domain provide a basis for proposing a new mathematical model for ultrasonic waves interacting with poor interlayer bonds.

36 MATERIALS SCIENCE↗
NASA NTRS1997

Validity Using Pump-Probe Pulses to Determine the Optical Response of Niobate Crystals

A variety of niobate crystals have found their places in nonlinear optical applications as well as in laser devices. In recent years much attention has been paid to study the ultrafast optical response in a variety of photorefractive crystals such as KTa(1-x)Nb(x)O3 and KNbO3 crystals, glasses, semiconductors and polymers for applications in optical switching, information processing, optical computing, and all-optical device systems. Third-order optical nonlinearity is the most important property for realization of all-optical switching. Therefore experiments have been performed on the third order susceptibility using a variety of techniques such as the third-order harmonic generation, EFISH and degenerate four-wave mixing(DFWM). The latter has been conducted with a variety of pump wavelengths and with nanosecond, picosecond and femtosecond pulses. Niobate crystals, such as potassium niobate KNbO3, potassium tantalate niobate KTN family (KTa(1-x)Nb(x)O3), strontium barium niobate SBN (Sr(x)Ba(1-x)Nb2O6) and potassium-sodium niobate SBN (KNSBN) are attractive due to their photorefractive properties for application in optical storage and processing. The pulsed probe experiments performed on theses materials have suggested two types of time responses. These responses have been associated with an coherent response due to Chi(sup 3), and a long lived component due to excited state population. Recent study of DFWM on KNbO3 and KTN family reveals that the long lived component of those crystals depends on the crystal orientation. A slowly decaying signal is observable when the grating vector K(sub g) is not perpendicular to the C-axis of those photorefractive crystals', otherwise the optical response signal would be only a narrow coherent peak with FWHM equal to the cross-correlation width of the write beam pulses. Based on this understanding, we study the photodynamical process of a variety of niobate crystals using DFWM in a Kg perpindicular to C geometry with a ps-YAG: Nd laser operating at 532nm. However, the discrepancies in numerical estimations of Chi(sup 3) in these materials and other nonlinear optical media have resulted in a number of discussions. In order to better understand the photodynamical process of niobate crystals after an short pulse excitation we analyze the factors governing the coherent signal and present the DFWM spectra of niobate crystals.

Liu, Huimin↗
DOE OSTI2025

Topology-imprinting in nonlinear metasurfaces

Flat optical components, or metasurfaces, have transformed optical imaging, data storage, information processing, and biomedical applications by providing unprecedented control over light-matter interactions. These nano-engineered structures enable compact, multidimensional manipulation of light’s amplitude, phase, polarization, and wavefront, producing scalar and vector beams with unique properties such as orbital angular momentum and knotted topologies. This flexibility has potential applications in optical communication and imaging, particularly in complex environments such as atmospheric turbulence and undersea scattering. However, designing metasurfaces for shorter wavelengths, such as visible and ultraviolet light, remains challenging due to fabrication limitations and material absorption. Here, we introduce an innovative concept called topology imprinting using judiciously designed all-dielectric nonlinear optical metasurfaces to replicate desired waveforms at fundamental and harmonic frequencies, opening promising avenues for advanced photonic applications.

42 ENGINEERING↗
DOE OSTI2026

A Generalized Grain-Scale Model for the Non-Plasma and Plasma-Assisted Hydrogen Direct Reduction of Iron Ore

Direct Reduction of Iron ore using hydrogen (H-DRI) is a promising pathway towards efficient steelmaking and accurate predictive models are a necessity for scale-up and optimization of this technology. However, accurate models of this process remain limited because existing models oversimplify grain-scale phenomena, such as nonlinearity inside grain, self-sufficient porosity, surface reactions, and the role of plasma species. These phenomena are important for flash steelmaking and plasma-assisted H-DRI processes. To address this need, we present a phenomenological model for simulating H-DRI at the scale of a single micron-sized grain of the iron ore. We call this the Transient Reactive Grain Model (TRGM). TRGM incorporates key physical process: gas species transport, a chemical kinetics of material conversion, nanopore structural evolution and, adsorption-desorption surface kinetics at the reactive nanopore surface. The important contribution of this work is that the model provides a dependence on different reductant species, specifically hydrogen atoms versus molecules, so that role of hydrogen plasma reduction can be clarified compared to the use of pure hydrogen gas reduction. TRGM predictions agree well with experimental data for both molecular H2 reduction of Fe2O3 and plasma hydrogen reduction of Fe3O4. Results reveal species concentration gradients with a diffuse reaction zone, and enhanced hydrogen diffusion at the grain outer surface due to evolving porosity. These findings challenge common assumptions in existing models, including sharp reaction fronts, quasi-steady diffusion and kinetics, and the neglect of surface chemistry. As a generalized grain-scale model for H-DRI processes, TRGM has practical applications in flash steelmaking and in-flight reduction using both molecular and plasma hydrogen.

08 HYDROGEN↗
DOE OSTI2025

Understanding plasma turbulence through exact coherent structures

Plasma turbulence is a key challenge in understanding transport phenomena in magnetically confined plasmas. This work presents a generalized framework to analyze plasma turbulence that utilizes periodic orbit theory. In periodic orbit theory, doubly periodic solutions (coherent structures) of the governing equation(s) serve as building blocks of the considered turbulent dynamics. To illustrate the concept and method, the particularly simple Kuramoto–Sivashinsky (referred to here as LMRT for the original authors: LaQuey, Mahajan, Rutherford, and Tang) trapped-ion mode toy model is used. By applying numerical optimization techniques to the LMRT equation, we extract coherent spacetime patterns that represent the library of allowable fundamental structures of the equation. These structures provide a framework to systematically describe turbulence as a composition of recurrent solutions, revealing an underlying order within chaotic plasma motion. Although illustrated here using the simplified LMRT model for clarity, this framework provides a general strategy that can be extended to more complex and realistic models of plasma turbulence, including gyrokinetic systems. This offers a new method for predicting and potentially controlling transport processes in fusion plasmas by providing a bridge between nonlinear dynamical systems theory and plasma physics in the form of a generalized framework with which to analyze and understand spatially extended nonlinear partial differential equations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Quantum simulations for strong-field QED

Quantum field theory in the presence of strong background fields contains interesting problems where quantum computers may someday provide a valuable computational resource. In the noisy intermediate-scale quantum era it is useful to consider simpler benchmark problems in order to develop feasible approaches, identify critical limitations of current hardware, and build new simulation tools. Here we perform quantum simulations of strong-field QED (SFQED) in 3 + 1 dimensions, using real-time nonlinear Breit-Wheeler pair production as a prototypical process. The strong-field QED Hamiltonian is derived and truncated in the Furry-Volkov mode expansion, and the interactions relevant for Breit-Wheeler are transformed into a quantum circuit. Quantum simulations of a “null double slit” experiment are found to agree well with classical simulations following the application of various error mitigation strategies, including an asymmetric depolarization algorithm which we develop and adapt to the case of Trotterization with a time-dependent Hamiltonian. We also discuss longer-term goals for the quantum simulation of SFQED. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2020

Optimal Network Topology for Node-Breaker Representations With AC Power Flow Constraints

It has been demonstrated that network topology optimization (NTO) may change the topology of power system networks, and consequently, provide additional flexibility to reduce network congestion and violations. Most NTO problems are formulated based on the bus-branch model in which it is challenging to represent a realistic picture of all substation configurations. In this paper, we explore advantages of substation reconfiguration modeling based on node-breaker representations for NTO problem with full nonlinear alternating current power flow. It also proposes a tailored solution algorithm to solve this nonconvex mixed-integer nonlinear programming through the outer approximation method. The proposed solution approach iterates between a mixed-integer linear programming and a nonlinear subproblem. Additional enhancements to further accelerate the iteration process are illustrated. Numerical case studies demonstrate the relative economic and operational impact of optimal network topology with node-breaker representations.

42 ENGINEERING↗
DOE OSTI2023

Efficient GPU Implementation of Automatic Differentiation for Computational Fluid Dynamics

Many scientific and engineering applications require repeated calculations of derivatives of output functions with respect to input parameters. Automatic Differentiation (AD) is a method that automates derivative calculations and can significantly speed up code development. In Computational Fluid Dynamics (CFD), derivatives of flux functions with respect to state variables (Jacobian) are needed for efficient solutions of the nonlinear governing equations. AD of flux functions on graphics processing units (GPUs) is challenging as flux computations involve many intermediate variables that create high register pressure and require significant memory traffic because of the need to store the derivatives. This paper presents a forward-mode AD method based on multivariate dual numbers that addresses these challenges and simultaneously reduces the floating-point operation count. The dimension of the multivariate dual numbers is optimized for performance. The flux computations are restructured to minimize the number of temporary variables and reduce register pressure. For effective utilization of memory bandwidth, shared memory is used to store the local flux Jacobian. This AD implementation is compared with several other Jacobian implementations on an NVIDIA V100 GPU (V100). For three-dimensional perfect-gas compressible-flow equations implemented in a practical CFD code, the AD implementation of a flux Jacobian based on multivariate dual numbers of dimension 5 outperforms all other GPU AD implementations on V100. Its performance is comparable with the optimized hand-differentiated version. Finally, the implementation achieves 75% of the peak floating-point throughput and 61 % of the peak global device memory bandwidth usage.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2026

High speed laser-induced refractive index change in ophthalmic HEMA hydrogel using a line focus

Laser-induced refractive index change (LIRIC) techniques typically focus femtosecond laser pulses to a tight spot to achieve the high intensities required to alter ophthalmic materials. A line-focus implementation using 800 nm pulses at 1 kHz repetition rate has been shown to greatly increase the effective writing speed of the process, but the high pulse energies used led to nonlinear propagation effects, presenting challenges for the fabrication of ophthalmic correctors. By taking advantage of the lower effective order of multiphoton absorption for 514 nm pulses, power levels may be reduced, allowing overlapped pulse writing. We demonstrate that this allows high-quality phase bars and devices to be made at high speeds.

Manning, Zachary A. [University of Rochester, NY (↗
NASA NTRS1972

A comparison of the effects of small nonlinearities on several estimation schemes.

A variety of techniques are available for estimating the states of nonlinear dynamic systems from noisy data. These procedures are generally equivalent when applied to linear systems. This paper investigates the difference between several of these procedures in the presence of small dynamic and observational nonlinearities. In particular, it examines one least square batch processing algorithm, and three recursive algorithms similar to the Kalman filter. To first order, all the estimators have the same covariance. Expressions for the means, however, show that each estimator has a different bias. The examples presented show that the biases can be a strong function of such parameters as initial covariances and number of data points being considered.

Conrad, B.↗
NASA NTRS1975

Analysis of nonlinear structures via mode synthesis

An effective procedure for NASTRAN was developed that permits any number of substructures of any size to be synthesized for the purpose of developing normal modes of vibration of the complete structural system. The technique is extended to permit modal transient analysis of the subdivided system. This latter procedure permits the use of NASTRAN's ability to include nonlinear forces in the problem. The five-phase process is accomplished using standard NASTRAN rigid formats with problem-independent alter packages and DMAP sequences.

Gieseke, R. K.↗
NASA NTRS1975

A variable-metric algorithm employing linear and quadratic penalties

A variable-metric algorithm is described that uses both linear and quadratic penalty terms for handling nonlinear constraints. Quadratic penalty coefficients are adjusted in a process which maintains a positive-definite matrix of second partial derivatives of the function without generating the large positive eigenvalues which cause zigzagging and slow convergence. The schemes suggested use inferred second-order properties not only in terms of the variable metric of the Davidson-Fletcher-Powell algorithm (or its relatives) but by estimating of second directional derivatives by fitting cubics to various functions along search directions.

Kelley, H. J.↗
NASA NTRS1982

Application of a transonic potential flow code to the static aeroelastic analysis of three-dimensional wings

Since the aerodynamic theory is nonlinear, the method requires the coupling of two iterative processes - an aerodynamic analysis and a structural analysis. A full potential analysis code, FLO22, is combined with a linear structural analysis to yield aerodynamic load distributions on and deflections of elastic wings. This method was used to analyze an aeroelastically-scaled wind tunnel model of a proposed executive-jet transport wing and an aeroelastic research wing. The results are compared with the corresponding rigid-wing analyses, and some effects of elasticity on the aerodynamic loading are noted.

Whitlow, W., Jr.↗
NASA NTRS1982

Application of a transonic potential flow code to the static aeroelastic analysis of three-dimensional wings

A method for including elastic effects in steady, transonic wing analysis is presented. Since the aerodynamic theory is nonlinear, the method requires the coupling of two iterative processes - an aerodynamic analysis and a structural analysis. A full potential analysis code, FL022, is combined with a linear structural analysis to yield aerodynamic load distributions on and deflections of elastic wings. This method was used to analyze an aeroelastically-scaled wind tunnel model of a proposed executive-jet transport wing and an aeroelastic research wing. The results are compared with the corresponding rigid-wing analyses, and some effects of elasticity on the aerodynamic loading are noted.

Whitlow, W., Jr.↗
NASA NTRS1981

Multigrid solvers on parallel computers

Massively parallel computers, as considered in this investigation, are not yet available. However, a large-scale parallel computer cannot usefully be designed before the hypothetical algorithms which will employ it are studied. Most of the studies of parallel partial differential equations (PDE) solvers are based on solution techniques much slower (on sequential machines) than multigrid methods. Multigrid methods are highly parallelizable. Each of their processes can simultaneously be performed at all grid points. The present investigation is concerned with a preliminary exploration of the potential of multigrid, or, more generally, Multi-Level Adaptive Techniques (MLAT) on computers with many processors. Basic processes are considered, taking into account coarse-grid approximation, relaxation, coarse-grid corrections, full multigrid algorithms, nonlinear problems and eigenvalue problems, fine-to-coarse correction, and chains of problems. Details of parallel multigrid processing are also examined.

Brandt, A.↗
NASA NTRS1985

Induced emission of radiation near 2(omega sub e) by a synchrotron-maser instability

In the literature, the emission of radiation at 2(omega sub e), where omega sub e denotes the electron plasma frequency, is usually explained as having been produced by the nonlinear interaction of two Langmuir waves via a backscattering process. Since the emission is frequently observed in solar radio bursts, the mechanism has attracted considerable theoretical interest. In the present paper a model is proposed based on a synchrotron-maser instability excited by a hollow beam of moderately relativistic electrons in a plasma, in which the plasma frequency is much higher than the gyrofrequency. An important conclusion is that, as a result of this instability, unpolarized electromagnetic waves with frequencies near 2(omega sub e) may be amplified.

Wu, C. S.↗
NASA NTRS1988

Technology advances for Space Shuttle processing

One of the major initial tasks of the Space Systems Integration and Operations Research Applications (SIORA) Program was the application of automation and robotics technology to all aspects of the Shuttle tile processing and inspection system. The SIORA Program selected a nonlinear systems engineering methodology which emphasizes a team approach for defining, developing, and evaluating new concepts and technologies for the operational system. This is achieved by utilizing rapid prototyping testbeds whereby the concepts and technologies can be iteratively tested and evaluated by the team. The present methodology has clear advantages for the design of large complex systems as well as for the upgrading and evolution of existing systems.

Wiskerchen, M. J.↗