Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “nonlinear phenomena”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

A Full-Induction Magnetohydrodynamics Solver for Liquid Metal Fusion Blankets in Vertex-CFD

Multiphysics modeling of liquid metal fusion blankets, which produce tritium and convert energy of neutrons created via fusion reactions into heat, is crucial for predicting performance, ensuring structural integrity, and optimizing energy production. While traditional blanket modeling of liquid metal flows during normal steady operating conditions commonly employs the inductionless approximation of the magnetohydrodynamics (MHD) equations, transient scenarios, when the plasma-confining magnetic field varies on millisecond time scales, require a full-induction MHD approach that dynamically evolves the magnetic field via the time-dependent induction equation. This paper presents the formulation, implementation, and initial verification of a full-induction MHD solver integrated within the open-source Vertex-CFD framework, which aims to achieve tight multiphysics coupling, a flexible software design enabling easy extension and addition of physics models, and performance portability across computing platforms. The solver utilizes finite element spatial discretization, implicit Runge–Kutta time integration, and an inexact Newton method to solve the resulting discrete nonlinear system, leveraging Trilinos packages for efficient computation. Verification against selected benchmark problems demonstrates accuracy and robustness of the solver. Furthermore, when the solver is applied to an idealized blanket model in 2.5D and full 3D, results obtained with Vertex-CFD are in good agreement with recently published quasi-2D simulations. These findings establish a computational foundation for future simulations of transient MHD phenomena in liquid metal blankets with Vertex-CFD, and open avenues for future extensions and performance optimizations.

Endeve, Eirik [ORNL] (ORCID:0000000312519507)↗
NASA NTRS1993

Feasibility study: Monodisperse polymer particles containing laser-excitable dyes

The optical properties associated with small particles, which include aerosols, hydrosols and solid microspheres have an impact on several areas of science and engineering. Since the advent of high-speed computers and lasers, the interaction of light with matter in the form of small particles with a discontinuous optical boundary relative to the surroundings has been much better understood. Various nonlinear optical effects have been observed involving interaction of a laser beam with both solid microspheres and liquid microdroplets. These include observation of second and third harmonic generation, four wave mixing, optical visibility, two photon absorption, observation of stimulated emission and lasing, and Stimulated Raman Scattering. Many of these effects are observed with laser intensities which are orders of magnitude less than that required by threshold condition for interactions in macroscopic bulk medium. The primary reason for this is twofold. The front surface of the microsphere acts as a thick lens to enhance the internal intensity of the input laser radiation, and the spherical shape of the droplet acts as an optical cavity to provide feedback at specific wavelengths corresponding to the whispering gallery modes or the morphology dependent resonances (MDR's). The most interesting and significant recent finding in this field is undoubtedly the existence of resonance peaks in linear and nonlinear optical spectra. Such resonance peaks are only dependent on the particle morphology, which means the size, shape and refractive index of the particle. Because of the simultaneous presence of these resonances, they have been referred to by many names, including structural resonances, whispering modes or whispering gallery modes, creeping waves, circumferential waves, surfaces modes, and virtual modes. All of these names refer to the same phenomena, i.e. morphology dependent resonances (MDR's) which has already been described and predicted precisely by electromagnetic theory and Loentz-Mie theory since 1908. MDR's can become important when the particle size (radius a) approaches and exceeds the wavelength of the electromagnetic wave (lambda) and the refractive index of the particle is greater than that of the surrounding medium. Such resonances are easiest to observe from a single particle with high symmetry, such as a sphere, spheroid, or cylinder. MDR's correspond to solutions of the characteristic equations of the electromagnetic fields in the presence of a sphere.

Venkateswarlu, Putcha↗
NASA NTRS1994

Sqeezing generated by a nonlinear master equation and by amplifying-dissipative Hamiltonians

In the first part of this contribution we show that the master equation derived from the generalized version of the nonlinear Doebner-Goldin equation leads to the squeezing of one of the quadratures. In the second part we consider two familiar Hamiltonians, the Bateman- Caldirola-Kanai and the optical parametric oscillator; going back to their classical Lagrangian form we introduce a stochastic force and a dissipative factor. From this new Lagrangian we obtain a modified Hamiltonian that treats adequately the simultaneous amplification and dissipation phenomena, presenting squeezing, too.

Dodonov, V. V.↗
DOE OSTI2024

A reduced-order model for nonlinear radiative transfer problems based on moment equations and POD-Petrov-Galerkin projection of the normalized Boltzmann transport equation

A data-driven projection-based reduced-order model (ROM) for nonlinear thermal radiative transfer (TRT) problems is presented. The TRT ROM is formulated by (i) a hierarchy of low-order quasidiffusion (aka variable Eddington factor) equations for moments of the radiation intensity and (ii) the normalized Boltzmann transport equation (BTE). The multilevel system of moment equations is derived by projection of the BTE onto a sequence of subspaces which represent elements of the phase space of the problem. Exact closure for the moment equations is provided by the Eddington tensor. A Petrov-Galerkin (PG) projection of the normalized BTE is formulated using a proper orthogonal decomposition (POD) basis representing the normalized radiation intensity over the whole phase space and time. The Eddington tensor linearly depends on the solution of the normalized BTE. By linear superposition of the POD basis functions, a low-rank expansion of the Eddington tensor is constructed with coefficients defined by the PG projected normalized BTE. The material energy balance (MEB) equation is coupled with the effective gray low-order equations which exist on the same dimensional scale as the MEB equation. The resulting TRT ROM is structure and asymptotic preserving. A detailed analysis of the ROM is performed on the classical Fleck-Cummings (F-C) TRT multigroup test problem in 2D geometry. Numerical results are presented to demonstrate the ROM's effectiveness in the simulation of radiation wave phenomena. Importantly, the ROM is shown to produce solutions with sufficiently high accuracy while using low-rank approximation of the normalized BTE solution. Essential physical characteristics of supersonic radiation wave are preserved in the ROM solutions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Fascinating physics at the edge of magnetic fusion devices

Abstract The physics of the processes at the edge of magnetic fusion devices is multifaceted and exhibits complex, nonlinear synergistic effects. Even though this region occupies only a small portion of the whole device, it plays a crucial role in overall plasma confinement, heat exhaust, and plasma–wall interactions. The latter affects not only the performance but also the lifetime of plasma-facing components and therefore remains an outstanding challenge for future fusion reactors. At the edge of fusion devices, researchers are dealing with phenomena including classical and anomalous plasma transport, atomic physics effects, and physics of plasma-facing material under strong irradiation by particle and energy fluxes. The diversity of edge physics makes it particularly attractive for young scientists. Working in this field, they can find endless possibilities to demonstrate their talents and creativity. This short review describes just some of the basic scrape-off layer and divertor plasma phenomena including divertor plasma detachment, intermittent bursts of anomalous cross-field plasma transport, plasma–material interactions, and dust in fusion plasmas, which are of particular interest in relation to fusion reactors.

Physics↗
DOE OSTI2026

Large-scale white noise and cosmology

The generation of white noise on large scales is a generic property of the dynamics of physical systems described by local nonlinear partial differential equations. Nonlinearities prevent the small-scale dynamics from being erased by smoothing. Unresolved small-scale dynamics act as an uncorrelated (white or Poissonian) noise (seemingly stochastic but actually deterministic) contribution to large-scale dynamics. This white noise exists even when the dynamics is very nearly linear. In cases where the power spectrum is sub-Poissonian on large scales, this noise will dominate on the largest scale power no matter the amplitude of the inhomogeneities. Such is the case in the standard model of cosmology, where the primordial density power spectrum is expected to have an almost Harrison-Zel’dovich, P [ k ] ∼ k , spectrum on a much broader range of scales than can be observed. Even though linear gravitational evolution dominates nonlinear corrections by a factor of ∼ 10 5 , the nonobservation of white noise on the Hubble scale precludes the extrapolation of this power law below the comoving 1 pc scale. More generally, observation or nonobservation of large-scale white noise provides a powerful probe of the Universe on very small scales in the early Universe. Gravitational radiation, phase transitions, vorticity, and running of the spectral index are all phenomena that can be probed with large-scale white noise. Large-scale white noise is a nonoptional feature of all cosmological models but one which has not heretofore been appreciated.

Barenboim, Gabriela [Valencia U., IFIC; Valencia U↗
DOE OSTI2020

MASTODON: An Open-Source Software for Seismic Analysis and Risk Assessment of Critical Infrastructure

Seismic analysis and risk assessment of safety-critical infrastructure like hospitals, nuclear power plants, dams, and facilities handling radioactive materials involve computationally intensive numerical models and coupled multiphysics scenarios. They are also performed in a strict regulatory environment that requires high software quality assurance standards, and in the case of safety-related nuclear facilities, a conformance to the American Society of Mechanical Engineers Nuclear Quality Assurance (NQA-1) standard. This paper introduces the open-source finite-element software, MASTODON (Multi-hazard Analysis of Stochastic Time-Domain Phenomena), which implements state-of-the-art seismic analysis and risk assessment tools in a quality-controlled environment. MASTODON is built on MOOSE (Multi-physics Object-Oriented Simulation Environment), which is a highly parallelizable, NQA-1 conforming, coupled multiphysics, finite-element framework developed at Idaho National Laboratory. MASTODON is capable of fault rupture and source-to-site wave propagation using the domain reduction method, nonlinear site response, and soil-structure interaction analysis, implicit and explicit time integration, automated stochastic simulations, and seismic probabilistic risk assessment. When coupled with other MOOSE applications, MASTODON can also solve strongly and weakly coupled multiphysics problems. This paper presents a summary of the capabilities of MASTODON and some demonstrative examples.

42 ENGINEERING↗
NASA NTRS1994

Nonlinear Resonances in the Solar System

Orbital resonances are ubiquitous in the Solar system. They play a decisive role in the long term dynamics, and in some cases the physical evolution, of the planets and of their natural satellites, as well as the evolution of small bodies (including dust) in the planetary system. The few-body gravitational problem of hierarchical planetary-type systems allows for a complex range of dynamical timescales, from the fast orbital periods to the very slow orbit precession rates. The interaction of fast and slow degrees of freedom produces a rich diversity of resonance phenomena. Weak dissipative effects - such as tides or radiation drag forces - also produce unexpectedly rich dynamical behaviors. This paper provides a mostly qualitative discussion of simple dynamical models for the commonly encountered orbital resonance phenomena in the Solar system.

Malhotra, Renu↗
NASA NTRS1983

Electrostatic Kelvin-Helmholtz instability in a radially injected plasma cloud

An elecrostatic finite-sized particle simulation model is used to study the early time-scale phenomena asociated with an impulsively injected plasma, expanding radially, normal to a strong ambient magnetic field. Results of the simulation show the early formation of a radial polarization electric field due to ion-electron charge separation which causes electrons to drift in the azimuthal direction. Velocity shear within this motion gives rise to a Kelvin-Helmholtz (dioxotron) instability, creating a radial fluted pattern as the plasma expands. In the nonlinear stage of instability, azimuthal variations in the electric field are observed which cause electrons to drift across the magnetic field, thereby reducing the space charge created by energetic ion expansion. The results of a linear stability analysis reveal that for a decreasing amount of charge separation in the plasma, the number of unstable azimuthal modes at maximum growth rate increases. Also, as the thickness of the electron ring increases, a fewer number of unstable modes develops for a given amount of charge separation. It is concluded that electrons become unstable at the very early time scale, within an ion gyroperiod, and that unstable electrons lead to modification of the ion dynamics.

Sydora, R. D.↗
NASA NTRS1988

The role of unsteady aerodynamics in aeroacoustics

The role of acoustics and unsteady aerodynamics research in understanding the fundamental physics of time-dependent fluid phenomena is reviewed. The key issues are illustrated by considering the sound radiation of turbulent jets and the aeroacoustics of rotating bodies such as helicopter rotors. The importance of computational methods as a link between aerodynamics and acoustics is also discussed. It is noted that where acoustic analogy techniques are sufficiently accurate, unsteady aerodynamics can be used for acoustic prediction. In supersonic problems where acoustics and aerodynamics are coupled, an integrated nonlinear analysis can provide an accurate problem solution.

Pao, S. Paul↗
NASA NTRS1988

Linear response theory for annealing of radiation damage in semiconductor devices

A theoretical study of the radiation/annealing response of MOS ICs is described. Although many experiments have been performed in this field, no comprehensive theory dealing with radiation/annealing response has been proposed. Many attempts have been made to apply linear response theory, but no theoretical foundation has been presented. The linear response theory outlined here is capable of describing a broad area of radiation/annealing response phenomena in MOS ICs, in particular, both simultaneous irradiation and annealing, as well as short- and long-term annealing, including the case when annealing is nearing completion. For the first time, a simple procedure is devised to determine the response function from experimental radiation/annealing data. In addition, this procedure enables us to study the effect of variable temperature and dose rate, effects which are of interest in spaceflight. In the past, the shift in threshold potential due to radiation/annealing has usually been assumed to depend on one variable: the time lapse between an impulse dose and the time of observation. While such a suggestion of uniformity in time is certainly true for a broad range of radiation annealing phenomena, it may not hold for some ranges of the variables of interest (temperature, dose rate, etc.). A response function is projected which is dependent on two variables: the time of observation and the time of the impulse dose. This dependence on two variables allows us to extend the theory to the treatment of a variable dose rate. Finally, the linear theory is generalized to the case in which the response is nonlinear with impulse dose, but is proportional to some impulse function of dose. A method to determine both the impulse and response functions is presented.

Litovchenko, Vitaly↗
NASA NTRS1992

Investigation of forced unsteady separated flows using velocity-vorticity form of Navier-Stokes equations

The phenomenon of forced unsteady separation and eruption of boundary-layer vorticity is a highly-complex, high-Reynolds number flow phenomenon, which abruptly leads to the formation of a dynamic stall vortex as demonstrated earlier by the authors for a NACA 0015 airfoil undergoing constant rate pitch-up motion. This, as well as the results of other researchers, have convincingly demonstrated a complex vortical structure within the state of unsteady separation prior to the evolution of dynamic stall. This phenomenon of vortex eruption, although observed in studying dynamic stall phenomena, is also associated with transition from laminar to turbulence flow and its generic nature has been stressed by many researchers including the present investigators. An unsteady Navier-Stokes (NS) analysis is developed for arbitrarily maneuvering bodies using velocity-vorticity variables; this formulation is nearly form-invariant under a generalized non-inertial coordinate transformation. A fully-implicit uniformly second-order accurate method is used, with the nonlinear convective terms approximated using a biased third-order upwind differencing scheme to be able to simulate higher-Re flows. No explicit artificial dissipation is added. The numerical method is fully vectorized and currently achieves a computational index of 7 micro-seconds per time step per mesh point, using a single processor on a CRAY Y-MP. The simulation results show that the energetic free shear from the leading edge is responsible for the wall viscous layer to abruptly erupt near the center of the counterclockwise rotating eddy in the unsteady boundary layer. Primary, secondary, tertiary and quaternary vortices have been observed before the dynamic stall vortex evolves and gathers its maximum strength. This study will discuss the simulation results of Reynolds number up to Re = 45,000 and will also discuss the efforts of initial acceleration in a specific maneuver, on the evolution of the stall vortex.

Ghia, K. N.↗
DOE OSTI2022

Scientific challenges to characterizing the wind resource in the marine atmospheric boundary layer

Abstract. With the increasing level of offshore wind energy investment, it is correspondingly important to be able to accurately characterize the wind resource in terms of energy potential as well as operating conditions affecting wind plant performance, maintenance, and lifespan. Accurate resource assessment at a particular site supports investment decisions. Following construction, accurate wind forecasts are needed to support efficient power markets and integration of wind power with the electrical grid. To optimize the design of wind turbines, it is necessary to accurately describe the environmental characteristics, such as precipitation and waves, that erode turbine surfaces and generate structural loads as a complicated response to the combined impact of shear, atmospheric turbulence, and wave stresses. Despite recent considerable progress both in improvements to numerical weather prediction models and in coupling these models to turbulent flows within wind plants, major challenges remain, especially in the offshore environment. Accurately simulating the interactions among winds, waves, wakes, and their structural interactions with offshore wind turbines requires accounting for spatial (and associated temporal) scales from O(1 m) to O(100 km). Computing capabilities for the foreseeable future will not be able to resolve all of these scales simultaneously, necessitating continuing improvement in subgrid-scale parameterizations within highly nonlinear models. In addition, observations to constrain and validate these models, especially in the rotor-swept area of turbines over the ocean, remains largely absent. Thus, gaining sufficient understanding of the physics of atmospheric flow within and around wind plants remains one of the grand challenges of wind energy, particularly in the offshore environment. This paper provides a review of prominent scientific challenges to characterizing the offshore wind resource using as examples phenomena that occur in the rapidly developing wind energy areas off the United States. Such phenomena include horizontal temperature gradients that lead to strong vertical stratification; consequent features such as low-level jets and internal boundary layers; highly nonstationary conditions, which occur with both extratropical storms (e.g., nor'easters) and tropical storms; air–sea interaction, including deformation of conventional wind profiles by the wave boundary layer; and precipitation with its contributions to leading-edge erosion of wind turbine blades. The paper also describes the current state of modeling and observations in the marine atmospheric boundary layer and provides specific recommendations for filling key current knowledge gaps.

17 WIND ENERGY↗
NASA NTRS1981

Simulation of space plasma phenomena

The importance of computer modeling in investigating physical phenomena is discussed, noting that modern computers can follow the time evolution of systems containing many millions of degrees of freedom, all of which are simultaneously interacting with each other. The two types of models used in simulating plasma, that is, particle and fluid, are described. Attention is given to particle models and some recent results obtained concerning the diffuse aurora. The simulations are shown to give clear evidence that a type of nonlinear cyclotron resonance process for the cold particles is responsible for the saturation observed in these calculations. It is also responsible for the cold electron heating and it is found that at saturation the rate of heating of the cold electrons essentially balances the rate at which energy is fed into unstable waves.

Dawson, J. M.↗
NASA NTRS1993

Projections on structures and material strength in the computational context

The ideas associated with the subsequent viewgraphs are summarized. The primary motivation behind this presentation is to observe that certain macroscopic, microscopic, and submicroscopic phenomena are being understood that have basic influence on understanding the durability and high temperature sensitivity of polymers and polymer-based composites. This understanding covers important issues of long term stability with respect to residual stresses and deformations which can have very deleterious effects on structures used for long periods of time as a result of the heat-involving manufacturing process. Beyond this, important progress is being made in understanding the nonlinear material response of polymers in the fracture context, because the nonlinear mechanics of the material at the tip of a crack, either propagating or ready to do so, is being understood with increasing precision.

Knauss, Wolfgang G.↗
DOE OSTI2022

Exact Coherent Structures and Phase Space Geometry of Preturbulent 2D Active Nematic Channel Flow

Confined active nematics exhibit rich dynamical behavior, including spontaneous flows, periodic defect dynamics, and chaotic “active turbulence.” Here, we study these phenomena using the framework of exact coherent structures, which has been successful in characterizing the routes to high Reynolds number turbulence of passive fluids. Exact coherent structures are stationary, periodic, quasiperiodic, or traveling wave solutions of the hydrodynamic equations that, together with their invariant manifolds, serve as an organizing template of the dynamics. We compute the dominant exact coherent structures and connecting orbits in a preturbulent active nematic channel flow, which enables a fully nonlinear but highly reduced-order description in terms of a directed graph. Using this reduced representation, we compute instantaneous perturbations that switch the system between disparate spatiotemporal states occupying distant regions of the infinite-dimensional phase space. Overall, our results lay the groundwork for a systematic means of understanding and controlling active nematic flows in the moderate- to high-activity regime.

36 MATERIALS SCIENCE↗
NASA NTRS2011

Discretely Conservative Finite-Difference Formulations for Nonlinear Conservation Laws in Split Form: Theory and Boundary Conditions

Simulations of nonlinear conservation laws that admit discontinuous solutions are typically restricted to discretizations of equations that are explicitly written in divergence form. This restriction is, however, unnecessary. Herein, linear combinations of divergence and product rule forms that have been discretized using diagonal-norm skew-symmetric summation-by-parts (SBP) operators, are shown to satisfy the sufficient conditions of the Lax-Wendroff theorem and thus are appropriate for simulations of discontinuous physical phenomena. Furthermore, special treatments are not required at the points that are near physical boundaries (i.e., discrete conservation is achieved throughout the entire computational domain, including the boundaries). Examples are presented of a fourth-order, SBP finite-difference operator with second-order boundary closures. Sixth- and eighth-order constructions are derived, and included in E. Narrow-stencil difference operators for linear viscous terms are also derived; these guarantee the conservative form of the combined operator.

Fisher, Travis C.↗
DOE OSTI2024

Nonlinear wavepacket dynamics in proximity to a stationary inflection point

A stationary inflection point (SIP) in the Bloch dispersion relation of a periodic waveguide is an exceptional point degeneracy where three Bloch eigenmodes coalesce, forming the so-called frozen mode with a divergent amplitude and vanishing group velocity of its propagating component. For this work, we have developed a theoretical framework to study the time evolution of wavepackets centered at an SIP. Analysis of the evolution of the statistical moment distribution of linear pulses shows a strong deviation from the conventional ballistic wavepacket dynamics in dispersive media. The presence of nonlinear interactions dramatically changes the situation, resulting in a mostly ballistic propagation of nonlinear wavepackets with the speed and even the direction of propagation essentially dependent on the wavepacket amplitude. Such a behavior is unique to nonlinear wavepackets centered at an SIP and can be used for the realization of a novel family of beam power routers for classical waves.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗