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 1,405 records · Page 78
All sourcesNASA NTRSDOE OSTI
NASA NTRS2021

Micromechanics-based Modeling of Laminated SiC/SiC Ceramic Matrix Composites

The behavior and response of ceramic matrix composites (CMCs), in particular silicon carbide fiber reinforced silicon carbide matrix (SiC/SiC), is affected by many factors such as variation of fiber volume fraction, residual stresses resulting from processing of the composites at high temperature, random microstructures, and the presence of matrix flaws (e.g., voids, pores, cracks etc.) as well as general material nonlinearity and heterogeneity that occurs randomly in a composite. Residual stresses arising from the phase change of constituents are evaluated in this paper and it is shown that they do influence composite strength and need to be properly accounted for. Additionally, the microstructures (location of fiber centers, coating thickness etc.) of advanced CMCs are usually disordered (or random) and fiber diameter and strength typically have a distribution. They rarely resemble the ordered fiber packing (square, rectangular, or hexagonal) that is generally assumed in micromechanics-based models with periodic boundary conditions for computational expediency. These issues raise the question of how should one model such systems effectively? Can an ordered hexagonal packed repeating unit cell (RUC) accurately represent the random microstructure behavior? How many fibers need to be included to enable accurate representation? Clearly, the number of fibers within an RUC must be limited to insure a balance between accuracy and efficiency. NASA’s in-house micromechanics-based code MAC/GMC provides a framework to analyze such RUCs for the overall composite behavior and the FEAMAC computer code provides linkage of MAC/GMC to the commercial FEA code, ABAQUS. The appropriate level of discretization of the RUC as well as the analysis method employed, i.e., Generalized Method of Cells (GMC) or High Fidelity Generalized Method of Cells (HFGMC), is investigated in this paper in the context of a unidirectional as well as a cross-ply laminated CMC. Results including effective composite properties, proportional limit stress (an important design parameter) and fatigue are shown utilizing both GMC as well as HFGMC. Finally, a few multiscale analyses are performed on smooth bar test coupons as well as test coupons with features such as open-hole and double notches using FEAMAC. Best practices and guidance are provided to take these phenomena into account and keep a proper balance between fidelity (accuracy) and efficiency. Following these guidelines can account for important physics of the problem and provide significant advantages when performing large multiscale composite structural analyses.

Ceramic Matrix Composites↗
NASA NTRS2022

Seasonal Tidewater Glacier Terminus Oscillations Bias Multi-Decadal Projections of Ice Mass Change

Numerical, process-based simulations of tidewater glacier evolution are necessary to project future sea-level change under various climate scenarios. Previous work has shown that nonlinearities in tidewater glacier and ice stream dynamics can lead to biases in simulated ice mass change in the presence of noisy forcings. Ice sheet modeling projections that are used in the IPCC Assessment Report 6 (AR6) utilize atmospheric and oceanic forcings at annual temporal resolution, omitting any higher frequency forcings. Here, we quantify the effect of seasonal (<1 year) tidewater glacier terminus oscillations on decadal-scale (30 years) mass change. We use an idealized geometry to mimic realistic tidewater glacier geometries and investigate the impact of the magnitude of seasonal oscillations, bed slope at the glacier terminus, and basal friction law. We find that omitting seasonal terminus motion results in biased mass change projections, with up to an 39% underestimate or 25% overestimate of mass loss when seasonality is neglected, depending on bed slope and magnitude of oscillations. The bias is most sensitive to the magnitude of the seasonal terminus oscillations and exhibits very little sensitivity to choice of friction law. Although omitting seasonal terminus motion may bias mass change projections of individual glaciers, further work is needed to investigate the impact on ice sheet mass change projections. In order to achieve this, seasonality in atmospheric and oceanic forcings must be adequately represented and observations of seasonal terminus positions and tidewater glacier thickness changes must be acquired to evaluate numerical models. Note: Zip folder contains all videos listed in the supplemental document and can be downloaded for viewing.

D. Felikson↗
NASA NTRS2004

Unified Models of Turbulence and Nonlinear Wave Evolution in the Extended Solar Corona and Solar Wind

The PI (Cranmer) and Co-I (A. van Ballegooijen) made substantial progress toward the goal of producing a unified model of the basic physical processes responsible for solar wind acceleration. The approach outlined in the original proposal comprised two complementary pieces: (1) to further investigate individual physical processes under realistic coronal and solar wind conditions, and (2) to extract the dominant physical effects from simulations and apply them to a 1D model of plasma heating and acceleration. The accomplishments in Year 2 are divided into these two categories: 1a. Focused Study of Kinetic Magnetohydrodynamic (MHD) Turbulence. lb. Focused Study of Non - WKB Alfven Wave Rejection. and 2. The Unified Model Code. We have continued the development of the computational model of a time-study open flux tube in the extended corona. The proton-electron Monte Carlo model is being tested, and collisionless wave-particle interactions are being included. In order to better understand how to easily incorporate various kinds of wave-particle processes into the code, the PI performed a detailed study of the so-called "Ito Calculus", i.e., the mathematical theory of how to update the positions of particles in a probabilistic manner when their motions are governed by diffusion in velocity space.

Cranmer, Steven R.↗
NASA NTRS1975

Effect of Nozzle Nonlinearities upon Nonlinear Stability of Liquid Propellant Rocket Motors

A three dimensional, nonlinear nozzle admittance relation is developed by solving the wave equation describing finite amplitude oscillatory flow inside the subsonic portion of a choked, slowly convergent axisymmetric nozzle. This nonlinear nozzle admittance relation is then used as a boundary condition in the analysis of nonlinear combustion instability in a cylindrical liquid rocket combustor. In both nozzle and chamber analyses solutions are obtained using the Galerkin method with a series expansion consisting of the first tangential, second tangential, and first radial modes. Using Crocco's time lag model to describe the distributed unsteady combustion process, combustion instability calculations are presented for different values of the following parameters: (1) time lag, (2) interaction index, (3) steady-state Mach number at the nozzle entrance, and (4) chamber length-to-diameter ratio. In each case, limit cycle pressure amplitudes and waveforms are shown for both linear and nonlinear nozzle admittance conditions. These results show that when the amplitudes of the second tangential and first radial modes are considerably smaller than the amplitude of the first tangential mode the inclusion of nozzle nonlinearities has no significant effect on the limiting amplitude and pressure waveforms.

Padmanabhan, M. S.↗
DOE OSTI2021

Non-intrusive data-driven model reduction for differential–algebraic equations derived from lifting transformations

In this paper we present a non-intrusive data-driven approach for model reduction of nonlinear systems. The approach considers the particular case of nonlinear partial differential equations (PDEs) that form systems of partial differential–algebraic equations (PDAEs) when lifted to polynomial form. Such systems arise, for example, when the governing equations include Arrhenius reaction terms (e.g., in reacting flow models) and thermodynamic terms (e.g., the Helmholtz free energy terms in a phase-field solidification model). Using the known structured form of the lifted algebraic equations, the approach computes the reduced operators for the algebraic equations explicitly, using straightforward linear algebra operations on the basis matrices. The reduced operators for the differential equations are inferred from lifted snapshot data using operator inference, which solves a linear least squares regression problem. The approach is illustrated for the nonlinear model of solidification of a pure material. The lifting transformations reformulate the solidification PDEs as a system of PDAEs that have cubic structure. The operators of the lifted system for this solidification example have affine dependence on key process parameters, permitting us to learn a parametric reduced model with operator inference. Numerical experiments show the effectiveness of the resulting reduced models in capturing key aspects of the solidification dynamics.

42 ENGINEERING↗
NASA NTRS1992

A study of equation solvers for linear and non-linear finite element analysis on parallel processing computers

Concurrent computing environments provide the means to achieve very high performance for finite element analysis of systems, provided the algorithms take advantage of multiple processors. The authors have examined several algorithms for both linear and nonlinear finite element analysis. The performance of these algorithms on an Alliant FX/80 parallel supercomputer has been studied. For single load case linear analysis, the optimal solution algorithm is strongly problem dependent. For multiple load cases or nonlinear analysis through a modified Newton-Raphson method, decomposition algorithms are shown to have a decided advantage over element-by-element preconditioned conjugate gradient algorithms.

Watson, Brian C.↗
DOE OSTI2022

Size-Dependent Photon Avalanching in Tm 3+ Doped LiYF 4 Nano, Micro, and Bulk Crystals

Photon avalanche (PA) is a highly nonlinear mode of upconversion that is characterized by 100–1000-fold increase in luminescence intensity upon minute increments of pumping power. The practical realization of numerous possible nano-bio-technology applications utilizing the PA phenomenon will require information on its susceptibility to the material volume and surface. In this report these parameters are investigated via experimental and theoretical PA. The two-color, highly nonlinear PA emission at 475 and 800 nm is clearly observed in bulk single crystal, individual microcrystals, and ensembles of colloidal core and core–shell nanoparticles of LiYF 4 host doped with either 3 or 8% of thulium ions. The properties of PA emission, such as PA nonlinearity, PA gain, PA intensity, and luminescence kinetics in these materials show dependence on crystal volume and surface quenching. Theoretical simulations provide understanding of key physical processes that influence PA performance. Moreover, photon avalanche single beam super-resolution imaging is realized for the first time in 3% Tm 3+ doped LiYF 4 core–shell nanoparticles. The obtained insights and predictions form a solid background for further development and applications of new optimized PA materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Detecting isolated resonance curves using fixed frequency voltage control tests

Isolated resonance curves, or isolas, are resonance branches of the harmonically forced system that exist separately from the main nonlinear forced response curve, leading to excessive vibrations. Traditional stepped or swept sine simulations and tests rely on continuation along the frequency parameter, typically resulting in a jump phenomenon along the primary resonance branch, prior to the disconnected isola. The main objective of this research is to propose an approach to identify isolated resonance curves by performing continuation along the input amplitude that initializes the response from a low-amplitude solution in the linear regime. Furthermore, this is achieved with the open-loop fixed frequency voltage control method that continues along the shaker voltage parameter and measures the so-called S-curves, which are theoretically a continuous solution branch that connect to the isola. The methodology is demonstrated on a fixture-wing-pylon assembly with a vibro-impact nonlinearity localized in a pylon subcomponent attachment. Multi-harmonic balance simulations are deployed to compute both the nonlinear forced response curves and S-curves to demonstrate the isola detection strategy on a reduced-order finite element model of the nonlinear system. Swept sine and fixed frequency voltage control tests are then conducted on the physical structure to demonstrate the isola detection experimentally, revealing the existence of the large amplitude vibrations that are undetected in the forces levels and frequencies measured with traditional frequency sweeping.

Characterization and Analytical Technique↗
DOE OSTI2022

Coupled effect of Cr and Al on interactions between a prismatic interstitial dislocation loop and an edge dislocation line in Fe-Cr-Al alloy

Segregation of alloying elements to a prismatic dislocation loop under irradiation is an important phenomenon in understanding the role of loops in radiation effects in alloys. In this work, various segregation sequences of Cr and Al atoms to a $\langle$100$\rangle$ and a 1/2 $\langle$111$\rangle$ prismatic dislocation loop in a Fe-Cr-Al alloy are explored using molecular dynamics method and ab initio energy calculation. Here, the results reveal that Cr can readily segregate to the loops while Al atoms cannot by themselves, but the segregated Cr atoms are able to promote Al segregation at the loops, defined as a coupled segregation of Cr and Al. After interacting with Cr and Al, a previously mobile 1/2 $\langle$111$\rangle$ loop becomes sessile. Al segregation through the coupled segregation process can further affect the pinning behavior of a dislocation loop to an edge dislocation motion. By varying the Al fraction in the segregated solutes, a nonlinear dependence of the pinning effect on dislocation motion is determined. At lower fractions, the pinning strength increases with the increasing Al fraction until Al fraction reaches a critical value of about 19%, above which the pinning strength decreases with further increase of Al fraction. All these results suggest that, in addition to the effect of Al on the surface resistance to oxidation, Al effect on the dislocation loop behavior is important to optimize the composition of Fe-Cr-Al alloys for their applications in accident tolerant fuel concepts.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Deep Learning Segmentation of Complex Features in Atomic-Resolution Phase-Contrast Transmission Electron Microscopy Images

Phase-contrast transmission electron microscopy (TEM) is a powerful tool for imaging the local atomic structure of materials. TEM has been used heavily in studies of defect structures of two-dimensional materials such as monolayer graphene due to its high dose efficiency. However, phase-contrast imaging can produce complex nonlinear contrast, even for weakly scattering samples. It is, therefore, difficult to develop fully automated analysis routines for phase-contrast TEM studies using conventional image processing tools. For automated analysis of large sample regions of graphene, one of the key problems is segmentation between the structure of interest and unwanted structures such as surface contaminant layers. In this study, we compare the performance of a conventional Bragg filtering method with a deep learning routine based on the U-Net architecture. We show that the deep learning method is more general, simpler to apply in practice, and produces more accurate and robust results than the conventional algorithm. We provide easily adaptable source code for all results in this paper and discuss potential applications for deep learning in fully automated TEM image analysis.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Data-Efficient Methods for Determining Flory–Huggins χ Parameters in Multicomponent Polymer Formulations

Polymer formulations are essential in diverse applications including personal care products, coatings, paints, adhesives, and plastic materials. Designing these formulations requires navigating large, complex design spaces, where phase and self-assembly behavior critically impact performance. The Flory–Huggins χ parameter, which quantifies segmental miscibility, is widely used to parametrize the excess free energy of mixing in formulation models. In this work, we introduce two data-efficient, top-down methods for estimating χ parameters using the Random Phase Approximation (RPA): (i) Boundary Nonlinear Regression (Boundary-NLR), which fits theoretical spinodal boundaries to experimental phase boundaries, and (ii) Surrogate Model Inverse Parameter Estimation (SMIPE), which uses a Gaussian Process Classifier to fit sparse phase maps via a surrogate model. Both methods allow rapid parametrization of polymer field-theoretic models without the need for additional experiments. We evaluate these approaches on data sets involving polymer–solvent–nonsolvent ternary mixtures and block copolymer–solvent systems, demonstrating their robustness to experimental noise and their relevance for real-world formulation design.

copolymers↗
DOE OSTI2020

A Two-Stage Decomposition Approach for AC Optimal Power Flow

The alternating current optimal power flow (AC-OPF) problem is critical to power system operations and planning, but it is generally hard to solve due to its nonconvex and large-scale nature. Furthermore, this paper proposes a scalable decomposition approach in which the power network is decomposed into a master network and a number of subnetworks, where each network has its own AC-OPF subproblem. This formulates a two-stage optimization problem and requires only a small amount of communication between the master and subnetworks. The key contribution is a smoothing technique that renders the response of a subnetwork differentiable with respect to the input from the master problem, utilizing properties of the barrier problem formulation that naturally arises when subproblems are solved by a primal-dual interior-point algorithm. Consequently, existing efficient nonlinear programming solvers can be used for both the master problem and the subproblems. The advantage of this framework is that speedup can be obtained by processing the subnetworks in parallel, and it has convergence guarantees under reasonable assumptions. The formulation is readily extended to instances with stochastic subnetwork loads. Numerical results show favorable performance and illustrate the scalability of the algorithm which is able to solve instances with more than 11 million buses.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Advances in SRF cavity architectures for quantum computing

Superconducting radio frequency (SRF) cavities offer an excellent platform for storing and processing quantum information due to their exceptionally long lifetimes and large accessible Hilbert spaces. A common strategy to manipulate the quantum states is to use a nonlinear element like a transmon. Nevertheless, constructing a 3D SRF architecture with sustained cavity lifetimes presents several challenges. In this talk, we present our successful integration of transmons with single-cell Nb SRF cavities, enabling the preparation of several non-classical states. Preliminary results reveal the potential of this platform to build a qudit-based quantum computer, departing from traditional qubit-based architectures. We discuss strategies to enhance coherence times, develop high-fidelity gate schemes, and expand the system for constructing a multi-qudit quantum processor. Furthermore, we explore how this architecture can bolster dark matter detection efficiency through multi-photon quantum state preparation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1978

The effect of ground hydrology on climate sensitivity to solar constant variations

A zonally averaged climate model is used to investigate the effects of different evaporation parameterizations on climate sensitivity to solar constant variations. The model is based on the annual mean zonally averaged hemispheric climate model of Ohring and Adler (1978), with some modifications on the heating parameterizations. The dynamic model utilizes the two-level quasi-geostrophic potential vorticity equations. The heating processes are related to solar radiation, long-wave radiation, convection, evaporation, latent heat release, and ocean transport. Evaporation parameterization has been improved by using a nonlinear formulation in which the Bowen ratio is determined by the vertical gradients of water vapor and temperature near the earth's surface.

Chow, S.-H.↗
NASA NTRS1978

Nonlinear and linear bottom interaction effects in shallow water

The paper examines wave-energy dissipation rates in shallow water calculated from measured wave spectra at different distances from the shore. Different linear and nonlinear transfer and dissipation mechanisms are discussed. The various data sets are interpreted in terms of prevailing mechanisms at the respective sites. The incorporation of different processes in a predictive shallow-water model is outlined. The analysis suggests that bottom motion is primarily responsible for wave-energy dissipation in the Delta Region of the Gulf of Mexico, that friction is mainly responsible for wave-energy dissipation in Marineland, Panama City and Melkbosstrand, and that percolation is probably the dominant mechanism in the JONSWAP area of the North Sea.

Shemdin, O.↗
NASA NTRS1980

The Seasat scanning multichannel microwave radiometer /SMMR/ - Radiometric calibration algorithm development and performance

Radiometric calibration algorithms for processing Seasat scanning multichannel microwave radiometer (SMMR) flight data are presented. An equation is derived which accounts for losses and reradiation in the microwave components and nonlinearities in the radiometer. The application of pre-launch calibration data to this derivation is described. A performance evaluation of the post-launch operation of the algorithm is presented, indicating temperature resolution of 0.63 K at 6.6 GHz and 1.39 K at 37 GHz.

Swanson, P. N.↗
NASA NTRS1980

Estimating the chlorophyll content in the waters of Guanabara Bay from the LANDSAT multispectral scanning digital data

Detection of water quality in Guanabara Bay using multispectral scanning digital data taken from LANDSAT satellites was examined. To test these processes, an empirical (statistical) approach was choosen to observe the degree of relationship between LANDSAT data and the in situ data taken simultaneously. The linear and nonlinear regression analyses were taken from among those developed by INPE in 1978. Results indicate that the major regression was in the number six MSS band, atmospheric effects, which indicated a correction coefficient of 0.99 and an average error of 6.59 micrograms liter. This error was similar to that obtained in the laboratory. The chlorophyll content was between 0 and 100 micrograms/liter, as taken from the MSS of LANDSAT.

Dejesusparada, N.↗
NASA NTRS1982

An experience in mesh generation for three-dimensional calculation of potential flow around a rotating propeller

An attempt is made to develop a three-dimensional, finite volume computational code for highly swept, twisted, small aspect ratio propeller blades with supersonic tip speeds, in a way that accounts for cascade effects, hub-induced flow, and nonlinear transonic effects. Attention is presently given to the generation of a computational mesh for such a complex propeller configuration, with the aim of sharing developmental process experience. The problem treated is unique, in that blade chord, blade length, hub length and blade-to-blade distance represent several characteristic length scales among which there is considerable disparity. An ad hoc mesh-generation scheme is accordingly developed.

Jou, W.-H.↗