Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nonlinear Dynamics and Pattern Formation”

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.

12 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Surface Engineering by Simultaneous Action of Multiple External Fields (Final Technical Report)

This research project aims at enabling surface engineering strategies based on the fundamental understanding and predictive modeling of the surface morphological response of solid materials subjected to the combined action of multiple external forces. Surface morphological response to mismatch strain in strained-layer heterostructures also is analyzed. Our research focuses on the surface morphological evolution and stability of electrically conducting (metallic, such as Cu) and semiconducting (such as Si) solids under the simultaneous application of mechanical stresses, electric fields, and temperature gradients. We emphasize on identifying the conditions under which the multiply driven surface morphology is stable and on exploring the complexity of the corresponding various morphologically stable surface patterns. We also emphasize on the driven evolution of nanoscale features aiming at stabilizing and controlling nanoscale patterns on surfaces by their manipulation through simultaneously applied multiple external forces. Specifically, we address systematically the morphological response to the combined action of multiple external forces of bulk solid surfaces, surfaces of thin films grown epitaxially on thick or thin substrates, as well as nanoscale surface features such as coherently strained islands grown epitaxially on substrate surfaces. The research is based on a modeling approach that combines theoretical analyses of surface morphological stability with self-consistent dynamical simulations of surface evolution based on properly parameterized continuum and multiscale surface transport models.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Black and gray box learning of amplitude equations: Application to phase field systems

Here, we present a data-driven approach to learning surrogate models for amplitude equations and illustrate its application to interfacial dynamics of phase field systems. In particular, we demonstrate learning effective partial differential equations describing the evolution of phase field interfaces from full phase field data. We illustrate this on a model phase field system, where analytical approximate equations for the dynamics of the phase field interface (a higher-order eikonal equation and its approximation, the Kardar-Parisi-Zhang equation) are known. For this system, we discuss data-driven approaches for the identification of equations that accurately describe the front interface dynamics. When the analytical approximate models mentioned above become inaccurate, as we move beyond the region of validity of the underlying assumptions, the data-driven equations outperform them. In these regimes, going beyond black box identification, we explore different approaches to learning data-driven corrections to the analytically approximate models, leading to effective gray box partial differential equations.

42 ENGINEERING↗
DOE OSTI2023

Exploring regular and turbulent flow states in active nematic channel flow via Exact Coherent Structures and their invariant manifolds

This work is a unified study of stable and unstable steady states of 2D active nematic channel flow using the framework of Exact Coherent Structures (ECS). ECS are stationary, periodic, quasiperiodic, or traveling wave solutions of the governing equations that, together with their invariant manifolds, organize the dynamics of nonlinear continuum systems. We extend our earlier work on ECS in the preturbulent regime by performing a comprehensive study of stable and unstable ECS for a wide range of activity values spanning the preturbulent and turbulent regimes. In the weakly turbulent regime, we compute more than 200 unstable ECS that co-exist at a single set of parameters, and uncover the role of symmetries in organizing the phase space geometry. We provide conclusive numerical evidence that in the preturbulent regime, generic trajectories shadow a series of unstable ECS before settling onto an attractor. Lastly, our studies hint at shadowing of quasiperiodic type ECS in the turbulent regime.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Beyond Ohm's law: Bernoulli effect and streaming in electron hydrodynamics

Recent observations of nonlocal transport in ultraclean two-dimensional materials raised the tantalizing possibility of accessing hydrodynamic correlated transport of a many-electron state. However, it has been pointed out that nonlocal transport can also arise from impurity scattering rather than interaction. At the crux of the ambiguity is the focus on linear effects, i.e., Ohm's law, which cannot easily differentiate among different modes of transport. Furthermore we propose experiments that can reveal rich hydrodynamic features in the system by tapping into the nonlinearity of the Navier-Stokes equation. Three experiments we propose will each manifest a unique phenomenon that is well known in classical fluids: the Bernoulli effect, Eckart streaming, and Rayleigh streaming. Analysis of known parameters confirms that the proposed experiments are feasible and the hydrodynamic signatures are within reach of graphene-based devices. Experimental realization of any one of the three phenomena will provide a stepping stone to formulating and exploring the notions of nonlinear electron fluid dynamics with an eye to celebrated examples from classical nonlaminar flows, e.g., pattern formation and turbulence.

2-dimensional systems↗
DOE OSTI2026

Emergence of complex-regular eutectic patterns in Al–Ge: Observations from correlative nano-imaging

Eutectic solidification exemplifies nonequilibrium pattern formation, making it a well-studied moving boundary problem. Yet the mechanisms behind the formation of complex-regular microstructures – particularly in highly anisotropic systems with a significant volume fraction of a faceted phase – remain poorly understood. Our understanding of such systems is made complicated by the nonlinear interface kinetics and unique growth dynamics characteristic of faceted phases. To address these challenges, we investigate a model Al–Ge eutectic system, where the faceted Ge phase constitutes a substantial volume fraction (~0.35) and where the two solid phases arrange into so-called “fishbone” or “feather” complex-regular patterns. Using synchrotron-based x-ray nano-imaging and nanotomography with high spatial resolution (22 nm per pixel), we capture in real-time the evolution of the solid–liquid interfaces and the resulting three-dimensional microstructures in this faceted/non-faceted eutectic system. By integrating these observations with electron backscattered diffraction, we elucidate the crystallographic biases on the solidification process and the mechanisms driving the formation of such complex-regular microstructures. These findings inform a new growth model for irregular eutectics in (near-)symmetrical phase diagrams, offering insight on advanced microstructural design and processing strategies. More broadly, we demonstrate how interfacial curvature is generated in irregular eutectic alloys and how it depends on the volume fraction of the faceted phase.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Pattern formation in a coupled membrane-bulk reaction-diffusion model for intracellular polarization and oscillations

Reaction-diffusion systems have been widely used to study spatio-temporal phenomena in cell biology, such as cell polarization. Coupled bulk-surface models naturally include compartmentalization of cytosolic and membrane-bound polarity molecules. Here we study the distribution of the polarity protein Cdc42 in a mass-conserved membrane-bulk model, and explore the effects of diffusion and spatial dimensionality on spatio-temporal pattern formation. We first analyze a one-dimensional (1-D) model for Cdc42 oscillations in fission yeast, consisting of two diffusion equations in the bulk domain coupled to nonlinear ODEs for binding kinetics at each end of the cell. In 1-D, our analysis reveals the existence of symmetric and asymmetric steady states, as well as anti-phase relaxation oscillations typical of slow-fast systems. We then extend our analysis to a two-dimensional (2-D) model with circular bulk geometry, for which species can either diffuse inside the cell or become bound to the membrane and undergo a nonlinear reaction-diffusion process. We also consider a nonlocal system of PDEs approximating the dynamics of the 2-D membrane-bulk model in the limit of fast bulk diffusion. In all three model variants we find that mass conservation selects perturbations of spatial modes that simply redistribute mass. In 1-D, only anti-phase oscillations between the two ends of the cell can occur, and in-phase oscillations are excluded. In higher dimensions, no radially symmetric oscillations are observed. Instead, the only instabilities are symmetry-breaking, either corresponding to stationary Turing instabilities, leading to the formation of stationary patterns, or to oscillatory Turing instabilities, leading to traveling and standing waves. Codimension-two Bogdanov—Takens bifurcations occur when the two distinct instabilities coincide, causing traveling waves to slow down and to eventually become stationary patterns. Our work clarifies the effect of geometry and dimensionality on behaviors observed in ma.ss-conserved cell polarity models.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Roadmap for Photonics with 2D Materials

Triggered by advances in atomic-layer exfoliation and growth techniques, along with the identification of a wide range of extraordinary physical properties in self-standing films consisting of one or a few atomic layers, two-dimensional (2D) materials such as graphene, transition metal dichalcogenides (TMDs), and other van der Waals (vdW) crystals now constitute a broad research field expanding in multiple directions through the combination of layer stacking and twisting, nanofabrication, surface-science methods, and integration into nanostructured environments. Photonics encompasses a multidisciplinary subset of those directions, where 2D materials contribute remarkable nonlinearities, long-lived and ultraconfined polaritons, strong excitons, topological and chiral effects, susceptibility to external stimuli, accessibility, robustness, and a completely new range of photonic materials based on layer stacking, gating, and the formation of moiré patterns. These properties are being leveraged to develop applications in electro-optical modulation, light emission and detection, imaging and metasurfaces, integrated optics, sensing, and quantum physics across a broad spectral range extending from the far-infrared to the ultraviolet, as well as enabling hybridization with spin and momentum textures of electronic band structures and magnetic degrees of freedom. The rapid expansion of photonics with 2D materials as a dynamic research arena is yielding breakthroughs, which this Roadmap summarizes while identifying challenges and opportunities for future goals and how to meet them through a wide collection of topical sections prepared by leading practitioners.

2D materials↗
DOE OSTI2025

Second harmonic Bessel-Gauss beam shaping with elliptic axicon aberrations

The second harmonic (SH) of an axicon generated Bessel-Gauss beam is created through the nonlinear interaction of photons with a crystal, resulting in the energy doubling of the output photons. In this work, we show experimentally that in addition to frequency doubling, the SH of Bessel-Gauss beams under asymmetric aberrations from an elliptic axicon exhibit intriguing beam formation. Particularly, the central region of the SH beam profile is composed of two central spots of various geometries surrounded by nested ellipses; one of which is the configuration of two central gamma dots with similar radius knotted by nested ellipses for a zeroth-order Bessel-Gauss pump. These SH beams consistently maintain their spatial profile throughout propagation, reminiscing the behavior of screw dislocations in wave patterns. Our numerical simulations produce beam dynamics consistent with that of experiments and further implicate the remarkable interweaving of bright spots with dark vortices. This is especially noticeable when the beams dynamically oscillate along the optical axis, resulting in the genesis of spatially polarized beams with a knotted framework. The insights gained from our study establish a novel paradigm for exploring interactions of Bessel-like beams with vortex dynamics. This, in turn, has the potential to spark innovations in optical applications, fostering new methodologies to harness and manipulate complex light structures. Our experimental findings could spur new methods of generating logical states of light and new opportunities for material processing control. Published by the American Physical Society 2025

Wang, Tianhong (ORCID:0000000169104201)↗
DOE OSTI2022

Scaling Behavior of Thermally Driven Fractures in Deep Low‐Permeability Formations: A Plane Strain Model With 1‐D Heat Conduction

Abstract Injection of cold fluids through/into deep formations may cause significant cooling, thermal stress, and possible thermal fracturing. In this study, the thermal fracturing of low‐permeability formations under one‐dimensional heat conduction was investigated using a plane strain model. Dimensionless governing equations, with dimensionless fracture length , aperture , spacing , time , and effective confining stress , were derived. Solution of single thermal fracture was derived analytically, while solution of multiple fractures with constant (or dynamic) spacing were obtained using the displacement discontinuity method (and stability analysis). For single fracture, increases nonlinearly with and then transitions to scaling law , indicating that late‐time fracture length increases linearly with the square root of cooling time. For constantly spaced fractures, deviates from the single‐fracture solution at a later for a larger , showing slower propagation under inter‐fracture stress interaction . For dynamically spaced fractures, fracture arrest induced by stress interaction was determined by the stability analysis; the fully transient solution provides evolution of dimensionless fracture length, spacing, aperture, and pattern; a similar scaling law, with , obtained shows the effect of both stress interaction and fracture arrest. The solution and scaling law provide fast predictions for all reservoir and cooling conditions using (single) model parameter . Application to a geothermal site with demonstrates that thermal fractures reach 0.67, 6.25, and 78.00 m in length, 0.49, 2.30, and 13.00 m in spacing, and 0.43, 2.09, and 12.19 mm in aperture at 1, 100, and 10,000 days.

Chen, Bin↗
DOE OSTI2026

Dynamic evolution of ultracold electron beams from nanostructured cathodes

Dynamic Coulomb expansion of dense, nonuniform particle beamlets can lead to causticlike formations in the charge density profile of the resultant beam. Such features only occur for cold beams where thermal diffusion of the particles is minimized. These phenomena have been observed experimentally for select cases with ion beams but not yet for electron beams as the electron dynamics occur over far faster timescales. We simulate the interaction of electron bunches while varying the initial transverse temperature and density profiles to determine the thresholds that characterize this pattern formation. Furthermore, we extend the simulations to geometric arrangements that are more pertinent to the nanotip cathode application and consider the effects of asymmetries and the impact of a low-density halo on the overall process. The simulations are conducted using a previously documented high-fidelity algorithm for collisional particle dynamics. These studies identify the temperature and density thresholds required for caustic formation and provide guidance for optimizing the design of nanoscale structured cathodes. We also demonstrate that the correlation between the final density patterns and the underlying beam parameters enables a virtual diagnostic tool for measuring ultralow emittance beams solely based on a transverse charge density profile.

43 PARTICLE ACCELERATORS↗
DOE OSTI2021

Edge-localized-mode simulation in CFETR steady-state scenario

We report the EPED1 model and self-consistent core-pedestal coupling in integrated modeling are used to design the pedestal structure of the China Fusion Engineering Testing Reactor (CFETR) steady-state scenario. The key parameters, such as β p and q 95 , are based on the grassy edge-localized-mode (ELM) experimental database. In this work, we use the BOUT++ six-field two-fluid code to simulate the onset of the ELM in the CFETR steady-state scenario. The ELM size is around 0.2% in nonlinear simulations, which is in the experimental range of the grassy ELM discharges, 0.1%–1% observed in multiple tokamak devices. Linear and nonlinear simulations show that the dominant high-n ballooning modes peak around n = 40. Compared to type-I ELM crashing dynamics, grassy ELM crashing has a smaller initial crash and is then followed by three phases of turbulence spreading, which are dominated by multi-modes, a high-n mode of n = 45 and low-n mode of n = 5, respectively. In contras to type-I ELM, the perturbation of the high-n mode has a narrow width around ψ = 0.95, and magnetic island formation and reconnection occur only beyond ψ = 0.95, leading to a small initial crash. Mode–mode interaction in the multi-mode coexistence stage stops the growth of individual modes and reduces the transport of particles and heat, and these are the two reasons why the ELM size is small. In–out asymmetry of transient heat flux with a ratio of E out /E in = 3.5 is found during grassy ELM crash. The rise and delay times of the heat flux match the calculation from the free-streaming model. To evaluate the erosion of the divertor target, the energy fluence at the outer divertor target is calculated, which is 0.029 MJ m -2 , 5.5 times smaller than the tungsten melting limit 0.16 MJ m -2 . The calculated energy fluency still follows the experimental scaling law from type-I ELM experiments. The fluctuation eddies in the toroidal direction show a filament structure at the outer mid-plane. Parallel heat flux patterns with a toroidal mode number n = 10 are found at the outer divertor with an amplitude of 680 MW m -2 .

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Soft interface instability and gas flow channeling in low-permeability deformable media

Understanding gas percolation through a clay layer or a shale formation is of great importance for the development of a geologic repository for nuclear waste disposal, a subsurface system for gas storage, and an engineering approach for hydrocarbon extraction from unconventional reservoirs. Gas injection experiments have revealed complex dynamic behaviours of gas percolation through water saturated compacted bentonite, characterized by a high breakthrough pressure, rapid breakthrough, a pressure/stress decay after the breakthrough, a relatively high migration rate, high-frequency periodic/nonperiodic variations in flow rate, stepwise rate reductions during relaxation, and low gas saturation over the whole process, all indicating channelling nature of the processes. Using linear stability analyses, we show that this channelling can autonomously emerge from the instability of the deformable interface between the injected gas and the compacted bentonite matrix driven by local stress concentration, pore dilation, and hydrologic gradient. Channel patterns formed would possess a fractal geometry. We further show that, once a percolating channel is established, the gas injected would percolate through the channel in a chain of gas bubbles, also due to the interface instability, resulting in periodic/chaotic variations in gas flow rate. Our work provides a unified explanation for key features observed for gas percolation in low-permeability deformable media. The work also suggests a possibility of designing an engineered barrier system for a nuclear waste repository that can have controllable gas release while limit water transport.

Bentonite↗