Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “classic nonlinearity”

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

Universality of a truncated sigma-model

Bosonic quantum field theories, even when regularized using a finite lattice, possess an infinite dimensional Hilbert space and, therefore, cannot be simulated in quantum computers with a finite number of qubits. A truncation of the Hilbert space is then needed and the physical results are obtained after a double limit: one to remove the truncation and another to remove the regulator (the continuum limit). A simpler alternative is to find a model with a finite dimensional Hilbert space belonging to the same universality class as the continuum model (a “qubitization”), so only the space continuum limit is required. A qubitization of the 1 + 1 dimensional asymptotically free O(3) nonlinear σ-model based on ideas of non-commutative geometry was previously proposed [1] and, in this paper, we provide evidence that it reproduces the physics of the σ-model both in the infrared and the ultraviolet regimes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Trivalent Rare Earth Adsorption at Phosphonic Acid Monolayers

The increasing need for rare earth separations requires a detailed understanding of trivalent ion behavior at charged aqueous interfaces. Here, neodymium (Nd) adsorption on Langmuir monolayers of octadecylphosphonic acid (ODPA), a single‐chain phosphonic acid capable of double deprotonation, at the air/water interface, is investigated. Combining sum frequency generation (SFG) spectroscopy with X‐ray fluorescence near total reflection (XFNTR), both the interfacial water ordering and ion density are examined. Under ambient conditions, Nd ions induce enhanced deprotonation of ODPA headgroups, leading to interfacial ion densities as high as 1 Nd per 30 Å 2 . This adsorption behavior arises from a complex interplay between direct electrostatic interactions, ion pairing, and hydration effects, which cannot be fully captured by classical Gouy–Chapman–Stern models. These insights into trivalent ion adsorption mechanisms provide a pathway toward more effective separation processes for rare earth metals.

adsorptions↗
DOE OSTI2024

Wall modes and the transition to bulk convection in rotating Rayleigh-Bénard convection

We investigate states of rapidly rotating Rayleigh-Bénard convection in a cylindrical cell over a range of Rayleigh numbers 3 × 10 5 ≤ Ra ≤ 5 × 10 9 and Ekman numbers 10 − 6 ≤ Ek ≤ 10 − 4 for Prandtl number Pr = 0.8 and aspect ratios 1 / 5 ≤ Γ ≤ 5 using direct numerical simulations. We characterize, for perfectly insulating sidewall boundary conditions, the first transition to convection via wall mode instability and the nonlinear growth and instability of the resulting wall mode states, including a secondary transition to time dependence. We show how the radial structure of the vertical velocity u z and the temperature T is captured well by the linear eigenfunctions of the wall mode instability where the radial width of u z is δ u z ∼ Ek 1 / 3 r / H whereas δ T ∼ e − k r ( k is the wave number of a laterally infinite wall mode state). The disparity in spatial scales for Ek = 10 − 6 means that the heat transport is dominated by the radial structure of u z since T varies slowly over the radial scale δ u z . We further describe how the transition to a state of bulk convection is influenced by the presence of the wall mode states. We use temporal and spatial scales as measures of the local state of convection and the Nusselt number Nu as representative of global transport. Our results elucidate the evolution of the wall state of rotating convection and confirm that wall modes are strongly linked with the boundary zonal flow being the robust remnant of nonlinear wall mode states. We also show how the heat transport ( Nu ) contributions of wall modes and bulk modes are related and discuss approaches to disentangling their relative contributions. Published by the American Physical Society 2024

58 GEOSCIENCES↗
DOE OSTI2025

Inference of phase field fracture models

The phase field approach to modeling fracture uses a diffuse damage field to represent cracks. This representation mollifies singularities that arise in computations with sharp interface models and some of the resultant difficulties in the mathematical and numerical treatment of fracture. Phase field fracture models have proven effective at representing crack propagation, branching, and merging. Specific formulations, beginning with brittle fracture, have also been shown to converge to classical solutions. Extensions to cover the range of material failure, including ductile and cohesive fracture, lead to an array of possible models. There exists a large body of literature focusing on this class of models and on the impact of model form on the predicted crack evolution. However, there have not been systematic studies into how optimal models may be chosen. Here, we take a first step in this direction by developing formal methods for identification of the best parsimonious model of phase field fracture given full-field data on the damage and deformation fields. We consider some of the main models that have been used for the degradation of elastic response due to damage and its propagation. Our approach builds upon Variational System Identification (VSI), a weak form variant of the Sparse Identification of Nonlinear Dynamics (SINDy). Furthermore, in this first communication we focus on synthetically generated data but we also consider central issues associated with the use of experimental full-field data, such as data sparsity and noise.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1981

Slosh dynamics of a spin-stabilized spacecraft comprising off-axis tanks filled partially with liquid propellant

The fundamental nonlinear equations of motion were derived and the specialized to a steady-state rotation of the vehicle about a given axis of rotation. A thrust about the spin axis was introduced. A perturbation solution was derived which linearizes the problem. The effect of the centrifugal and coriolis accelerations together with vorticity are implicitly taken into consideration in the formulation. A variational formulation of the associated boundary conditions is presented. For practical cases it is shown that the simple classical pendulum representation for slosh is not very appealing for a spinning spacecraft unless severe restrictions are allowed.

Fontenot, L. L.↗
NASA NTRS1980

Sound field near the edge of a liner in uniform flow

The classical concept of energy flux in sound propagation in a uniform medium without flow was extended to sound propagation with mean flow. A simple case of constant velocity mean flow is discussed for a two-dimensional duct with constant cross-section and a liner of finite length. Results do not show acoustic energy conservation and it is concluded that the liner converts flow energy into acoustic energy. It is important that results do not take into consideration nonlinear terms which are usually considered responsible for sound generation. Numerical computations estimated orders of magnitude and confirm the conclusions.

Moehring, W.↗
DOE OSTI2023

Multifidelity validation of digital surrogates using variable-density turbulent mixing models

In this study, ensembles of experimental data are presented and utilized to compare and validate two models used in the simulation of variable-density (Atwood = 0.22), compressible turbulent mixing. Though models of this kind (Reynolds averaged NavierStokes and large-eddy simulations) have been validated extensively with more canonical flows in previous studies, here the present approach offers novelty in the complexity of the geometry, the ensemble-based validation, and the uniformity of the computational framework on which the models are tested. Moreover, all experimental and computational tasks were completed by the authors which has led to a tightly coupled experimental configuration with its “digital twin.” The experimental divergent-shock-tube facility and its data acquisition methods are described and replicated in simulation space. A 2D Euler model which neglects the turbulent mixing at the interface is optimized to experimental data using a Gaussian process. This model then serves as the basis for both the 2D RANS and 3D LES studies that make comparisons to the mixing-layer data from the experiment. A relatively simple RANS model is shown to produce good agreement with experimental data only at late flow development times. The LES ensembles generally show good agreement with experimental data but display sensitivity to the characterization of initial conditions. Resolution-dependent behavior is also observed for certain higher-order statistics of interest. Overall, the LES model successfully captures the effects of divergent geometry, compressibility, and combined nonlinear instabilities inherent to the problem. The successful prediction of mixing width and its growth rate highlight the existence of three distinct regimes in the development of the instability, each with similarities to previously studied instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Dynamics and transport in the boundary-driven dissipative Klein-Gordon chain

Motivated by experiments on chains of superconducting qubits, we consider the dynamics of a classical Klein-Gordon chain coupled to coherent driving and subject to dissipation solely at its boundaries. As the strength of the boundary driving is increased, this minimal classical model recovers the main features of the "dissipative phase transition" seen experimentally. Between the transmitting and non-transmitting regimes on either side of this transition (which support ballistic and diffusive energy transport respectively), we observe additional dynamical regimes of interest. These include a regime of superdiffusive energy transport at weaker driving strengths, together with a "dissipative breather crystal" regime at stronger driving strengths, which is characterized by emergent translation symmetry, ballistic energy transport, and coherent oscillations of a nonlinear normal mode. Here, we propose a non-local Lyapunov exponent as an experimentally measurable diagnostic of many-body chaos in this system, and more generally in open systems that are only coupled to an environment at their boundaries.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2003

Analysis Method for Inelastic, Adhesively Bonded Joints with Anisotropic Adherends

A one-dimensional analysis method for evaluating adhesively bonded joints composed of anisotropic adherends and adhesives with nonlinear material behavior is presented in the proposed paper. The strain and resulting stress field in a general, bonded joint overlap are determined by using a variable-step, finite-difference solution algorithm to iteratively solve a system of first-order differential equations. Applied loading is given by a system of combined extensional, bending, and shear forces that are applied to the edge of the joint overlap. Adherends are assumed to behave as linear, cylindrically bent plates using classical laminated plate theory that includes the effects of first-order transverse shear deformation. Using the deformation theory of plasticity and a modified von-Mises yield criterion, inelastic material behavior is modeled in the adhesive layer. Results for the proposed method are verified against previous results from the literature and shown to be in excellent agreement. An additional case that highlights the effects of transverse shear deformation between similar adherends is also presented.

Smeltzer, Stanley S., III↗
DOE OSTI2021

Enhanced nonlinear interaction of polaritons via excitonic Rydberg states in monolayer WSe 2 .

Strong optical nonlinearities play a central role in realizing quantum photonic technologies. Exciton-polaritons, which result from the hybridization of material excitations and cavity photons, are an attractive candidate to realize such nonlinearities. While the interaction between ground state excitons generates a notable optical nonlinearity, the strength of such interactions is generally not sufficient to reach the regime of quantum nonlinear optics. Excited states, however, feature enhanced interactions and therefore hold promise for accessing the quantum domain of single-photon nonlinearities. Here we demonstrate the formation of exciton-polaritons using excited excitonic states in monolayer tungsten diselenide (WSe 2 ) embedded in a microcavity. The realized excited-state polaritons exhibit an enhanced nonlinear response ~[Formula: see text] which is ~4.6 times that for the ground-state exciton. The demonstration of enhanced nonlinear response from excited exciton-polaritons presents the potential of generating strong exciton-polariton interactions, a necessary building block for solid-state quantum photonic technologies.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Eigenstate Entanglement: Crossover from the Ground State to Volume Laws

For the typical quantum many-body systems that obey the eigenstate thermalization hypothesis (ETH), we argue that the entanglement entropy of (almost) all energy eigenstates is described by a single crossover function. The ETH implies that the crossover functions can be deduced from subsystem entropies of thermal ensembles and have universal properties. These functions capture the full crossover from the ground-state entanglement regime at low energies and small subsystem size (area or log-area law) to the extensive volume-law regime at high energies or large subsystem size. For critical one-dimensional systems, a universal scaling function follows from conformal field theory and can be adapted for nonlinear dispersions. We use it to also deduce the crossover scaling function for Fermi liquids in d > 1 dimensions. Here, the analytical results are complemented by numerics for large noninteracting systems of fermions in d ≤ 3 dimensions and have also been confirmed for bosonic systems and nonintegrable spin chains.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1983

Ion-beam-driven electrostatic ion cyclotron instabilities

Results are presented of a particle simulation study of the electrostatic ion-cyclotron (EIC) instability driven by a parallel ion beam. The results of this simulation study demonstrate the nonlinear consequences of nonresonant EIC waves destabilized by an ion beam parallel to the magnetic field for the case of a large beam velocity. As a consequence of the instability, it is shown that the beam ions are heated strongly in the perpendicular direction and suffer a strong anomalous friction via EIC waves which leads to the beam slowing down. Simulation results indicate that the anomalous slowing down of beam ions by EIC waves is much larger than that from the classical electron-drag, and the perpendicular collision frequency measured from perpendicular beam heating is as large as that from Bohm diffusion. It is concluded that the ion beam driven EIC wave is a viable mechanism for the transfer of ion parallel beam energy to the ion perpendicular energy.

Okuda, H.↗
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

Geometric learning for computational mechanics Part II: Graph embedding for interpretable multiscale plasticity

The history-dependent behaviors of classical plasticity models are often driven by internal variables evolved according to phenomenological laws. The difficulty to interpret how these internal variables represent a history of deformation, the lack of direct measurement of these internal variables for calibration and validation, and the weak physical underpinning of those phenomenological laws have long been criticized as barriers to creating realistic models. In this work, geometric machine learning on graph data (e.g. finite element solutions) is used as a means to establish a connection between nonlinear dimensional reduction techniques and plasticity models. Geometric learning-based encoding on graphs allows the embedding of rich time-history data onto a low-dimensional Euclidean space such that the evolution of plastic deformation can be predicted in the embedded feature space. Finally, a corresponding decoder can then convert these low-dimensional internal variables back into a weighted graph such that the dominating topological features of plastic deformation can be observed and analyzed.

42 ENGINEERING↗
NASA NTRS2013

Flutter Analysis of the Thermal Protection Layer on the NASA HIAD

A combination of classical plate theory and a supersonic aerodynamic model is used to study the aeroelastic flutter behavior of a proposed thermal protection system (TPS) for the NASA HIAD. The analysis pertains to the rectangular configurations currently being tested in a NASA wind-tunnel facility, and may explain why oscillations of the articles could be observed. An analysis using a linear flat plate model indicated that flutter was possible well within the supersonic flow regime of the wind tunnel tests. A more complex nonlinear analysis of the TPS, taking into account any material curvature present due to the restraint system or substructure, indicated that significantly greater aerodynamic forcing is required for the onset of flutter. Chaotic and periodic limit cycle oscillations (LCOs) of the TPS are possible depending on how the curvature is imposed. When the pressure from the base substructure on the bottom of the TPS is used as the source of curvature, the flutter boundary increases rapidly and chaotic behavior is eliminated.

Goldman, Benjamin D.↗
NASA NTRS2015

A Numerical and Experimental Study of Coflow Laminar Diffusion Flames: Effects of Gravity and Inlet Velocity

In this work, the influence of gravity, fuel dilution, and inlet velocity on the structure, stabilization, and sooting behavior of laminar coflow methane-air diffusion flames was investigated both computationally and experimentally. A series of flames measured in the Structure and Liftoff in Combustion Experiment (SLICE) was assessed numerically under microgravity and normal gravity conditions with the fuel stream CH4 mole fraction ranging from 0.4 to 1.0. Computationally, the MC-Smooth vorticity-velocity formulation of the governing equations was employed to describe the reactive gaseous mixture; the soot evolution process was considered as a classical aerosol dynamics problem and was represented by the sectional aerosol equations. Since each flame is axisymmetric, a two-dimensional computational domain was employed, where the grid on the axisymmetric domain was a nonuniform tensor product mesh. The governing equations and boundary conditions were discretized on the mesh by a nine-point finite difference stencil, with the convective terms approximated by a monotonic upwind scheme and all other derivatives approximated by centered differences. The resulting set of fully coupled, strongly nonlinear equations was solved simultaneously using a damped, modified Newton's method and a nested Bi-CGSTAB linear algebra solver. Experimentally, the flame shape, size, lift-off height, and soot temperature were determined by flame emission images recorded by a digital camera, and the soot volume fraction was quantified through an absolute light calibration using a thermocouple. For a broad spectrum of flames in microgravity and normal gravity, the computed and measured flame quantities (e.g., temperature profile, flame shape, lift-off height, and soot volume fraction) were first compared to assess the accuracy of the numerical model. After its validity was established, the influence of gravity, fuel dilution, and inlet velocity on the structure, stabilization, and sooting tendency of laminar coflow methane-air diffusion flames was explored further by examining quantities derived from the computational results.

microgravity↗
DOE OSTI2025

Local chromatic correction optics for Future Circular Collider e + e −

Local chromatic correction optics are proposed for the Future Circular e + − e − Collider. These new optics assume an identical layout of the magnets at all operation energies and include a new design for the arcs, for the straight sections, and for the final focus. The arcs design is a step forward from the classic FODO lattice that achieves near cancellation up to the fourth order of chromatic and geometric aberrations. Straight sections and the final focus benefit from the application of transparency conditions and are thus minimally impacting the global beam dynamics. The final focus design includes both vertical and horizontal chromaticity correction sections, crab sextupoles, and additional specific sextupoles and decapoles for the optimization of the dynamics off-energy. Decapoles are very effective in mitigating the reduction of dynamic aperture due to synchrotron radiation. Nonlinear magnet settings are further optimized with multiobjective algorithms. The overall properties of the new optics proposed are analyzed in the presence of errors and compared to other optics designs.

43 PARTICLE ACCELERATORS↗
NASA NTRS1985

Nonlinear quasi-fracture behavior of polymers

This paper deals with the nonlinear quasi-fracture behavior of polymers. Using the finite element method, the quasi-fracture model composed of an isolated craze opening with a nonlinear boundary envelope forces represented by oriented fibril domains is analyzed in simply-stressed rectangular plate. The stress distribution, the displacement profile, as well as the yield or plastic or flow zone adjacent to the craze tips are obtained. It is found that the variation of the stress distribution near the craze tip is not as sharp as that for an elastic crack. The initial load required for yielding or flowing is greater for the quasifracture as compared with that for a classical fracture. Under the same load the yield or flow zone is somewhat smaller as compared with that obtained on the basis of linear elastic-plastic fracture mechanics theory.

Sun, B. N.↗