Search NASA⌕ Search

SEARCH · Search NASA

Results for “topological fluid dynamics”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Noether currents for Eulerian variational principles in non-barotropic magnetohydrodynamics and topological conservations laws

We derive a Noether current for the Eulerian variational principle of ideal non-barotropic magnetohydrodynamics (MHD). It was shown previously that ideal non-barotropic MHD is mathematically equivalent to a five function field theory with an induced geometrical structure in the case that field lines cover surfaces and this theory can be described using a variational principle. Here we use various symmetries of the flow to derive topological constants of motion through the derived Noether current and discuss their implication for non-barotropic MHD.

42 ENGINEERING↗
DOE OSTI2025

Machine learning models for PDE constrained optimization

Partial differential equation (PDE)-constrained optimization problems arise in a variety of scientific and engineering applications, such as topology optimization, electrodynamics, fluid dynamics, and structural dynamics. However, these problems are often challenging and computationally expensive to solve, due to the need to solve the PDEs within the optimization loop. One approach to reducing the computational cost of these methods while providing convergence guarantees is through inexact trust region methods; this method uses lower fidelity solutions of the PDE at early stages of the optimization and adjusts the required accuracy of inexact PDE solvers as the optimization progresses. In this work, we explore the use of machine learning based surrogate models with these inexact trust region methods. We first demonstrate the potential of this approach by using Gaussian processes as the surrogate model and test this on a simple PDE-constrained optimization problem. We then document explorations into improving the computational costs of evolutional deep neural network / neural Galerkin methods, with the eventual goal of using these methods with the inexact trust region algorithms. We are able to speed up these approaches, albeit at the cost of lower accuracy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Liquid crystal torons in Poiseuille-like flows

Three-dimensional (3D) simulations of the structure of liquid crystal (LC) torons, topologically protected distortions of the LC director field, under material flows are rare but essential in microfluidic applications. Here, we show that torons adopt a steady-state configuration at low flow velocity before disintegrating at higher velocities, in line with experimental results. Furthermore, we show that under partial slip conditions at the boundaries, the flow induces a reversible elongation of the torons, also consistent with the experimental observations. These results are in contrast with previous simulation results for 2D skyrmions under similar flow conditions, highlighting the need for a 3D description of this LC soliton in relation to its coupling to the material flow. These findings pave the way for future studies of other topological solitons, like hopfions and heliknotons, in flowing soft matter systems.

Science & Technology - Other Topics↗
DOE OSTI2026

Orbital Inverse Faraday and Cotton-Mouton Effects in Hall Fluids

We report two light-induced orbital magnetization effects in quantum Hall (QH) fluids, stemming from their transverse response. The first is a purely transverse contribution to the inverse Faraday effect (IFE), where circularly polarized light induces a dc magnetization by stirring the charged fluid. This contribution dominates the IFE in the QH regime. The second is the orbital inverse Cotton-Mouton effect (ICME), in which linearly polarized light generates a dc magnetization. Since the applied field in the ICME does not break time-reversal symmetry, the induced magnetization directly probes the chiral orbital response of the fluid at the driving frequency. We estimate that the resulting magnetization lies in the range of 0.5–10 Bohr magnetons per charge carrier in materials such as graphene and transition-metal dichalcogenides (TMDs) in the QH regime. Finally, we show that the induced magnetization is accompanied by a local correction to the static particle density, enabling optical quantum printing of density profiles into the QH fluid.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Material coherence from trajectories via Burau eigenanalysis of braids

In this paper, we provide a numerical tool to study a material’s coherence from a set of 2D Lagrangian trajectories sampling a dynamical system, i.e., from the motion of passive tracers. We show that eigenvectors of the Burau representation of a topological braid derived from the trajectories have levelsets corresponding to components of the Nielsen–Thurston decomposition of the dynamical system. One can thus detect and identify clusters of space–time trajectories corresponding to coherent regions of the dynamical system by solving an eigenvalue problem. Unlike previous methods, the scalable computational complexity of our braid-based approach allows the analysis of large amounts of trajectories.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Aero-thermal numerical characterization of blunt fin-induced shock wave–boundary layer interaction and its control through leading-edge cooling injection

This study represents a novel evaluation of active flow control to alleviate the aerothermal penalties created by the blunt fin-induced shock wave–boundary layer interaction. The manuscript analyzes the effect of flow injection on a blunt fin-induced shock wave–boundary layer interaction via computational fluid dynamics simulations with various degrees of resolution. The impact on the mean flow topology and wall variables was investigated utilizing Reynolds-averaged Navier–Stokes simulations. Detached-eddy simulations revealed the low-frequency shock motion, shock wave–boundary layer, and horseshoe vortex interaction. The test article was exposed to two different incoming boundary layer thicknesses; the thicker boundary layer led to the appearance of larger turbulent scales. The Detached-eddy simulations revealed the time history of the shock wave–boundary layer interaction, focusing on the inception and development of the recirculated flow regions. Ultimately, spectral proper orthogonal decomposition was employed to identify the structures associated with the low-frequency shock motion caused by the shock wave–boundary layer interaction.

Lozano, Francisco (ORCID:0000000277871125)↗
DOE OSTI2022

Topology optimization of an airfoil fin microchannel heat exchanger using artificial intelligence

High-performance microchannel heat exchangers are needed to supply heat for power conversion for nuclear microreactors. An airfoil fin microchannel design, constructed of Alloy 617 with helium as the working fluid, is analyzed and optimized using a design of experiments with artificial intelligence techniques. The use of airfoil fins offers the potential to reduce pressure drop across the heat exchanger, as compared to other types of channel configurations. A framework for topology optimization of airfoil fin printed circuit heat exchangers (PCHEs) has been developed that can be readily extended to different fin sizes and shapes, as well as different inlet and operating conditions, materials of construction, and working fluids. An optimization procedure is developed that employs computational fluid dynamics for a set of design points identified using Latin hypercube sampling. Computational fluid dynamics is used to analyze a simplified two-channel configuration where five design parameters are varied – inlet angle, fin scale, extent of staggering, transverse and longitudinal pitches. Two methods (a 5D polynomial and a regression neural network) are compared for generating surrogate models and the resulting response surface approximation is input to a genetic algorithm that is used to identify a set of optimal parameters. The optimal geometries are found across six channel Reynolds numbers ranging from 1000 to 5000, since inlet conditions affect flow through the heat exchanger. Additionally, a set of optimal designs that maximizes heat transfer and minimizes pressure drop is identified, and a thermal stress analysis is performed on the optimal design. Correlations for the Nusselt number and Darcy friction factor are developed that can be useful for thermal hydraulic analyses using system codes. Thermal stresses are analyzed and a brief discussion of the status of code cases of PCHEs for nuclear applications is given. Testing and thermomechanical modeling is needed to facilitate future code compliance of PCHEs for high pressure and high temperature applications.

42 ENGINEERING↗
DOE OSTI2022

Multi-scale and multi-physics analysis, design optimization, and experimental validation of heat exchangers utilizing high performance, non-round tubes

Air-to-refrigerant heat exchangers (HXs) are fundamental components in HVAC&R systems, and considerable research has been dedicated designing continually smaller, lighter, and more efficient HX designs. In recent years, researchers have leveraged advancements in Computational Fluid Dynamics (CFD), Finite Element Analysis (FEA), and optimization algorithms to consider primary tube shape and topology optimization to design highly compact, high performance HXs for a multitude of applications. In this research, we present a computationally efficient, comprehensive, multi-scale, and multi-physics analysis and optimization method for air-to-refrigerant HXs featuring automated CFD and FEA simulations and approximation-assisted optimization. This methodology was utilized to design HXs with shape-optimized, non-round tubes which outperform current state-of-the-art tube-fin HXs without compromising structural integrity. The optimal HXs were shown to deliver similar thermal performance to the baseline HXs while also achieving more than 20% reductions in airside pressure drop and core envelope volume and more than 30% reduction in internal volume. Comprehensive experimental validation of the optimization methodology was conducted through the testing of two prototypes in a standardized wind tunnel facility under multiple operating conditions. For prototype #1 under dry evaporator conditions, the predicted heat load agreed within ± 10% of measured values and the predicted airside pressure drop agreed within ± 30%, while for dehumidifying conditions, the predicted sensible and latent heat loads agreed within ± 10% and ± 20% of the measured values, respectively. For prototype #2, the predicted condenser heat load agreed within ± 3.0% of measured values, and the predicted airside pressure drop agreed within ± 27%. The acceptable agreement between simulation and experimental results for the present application highlights the flexibility of the novel optimization methodology to design next generation HXs with improved performance and reduced volume, weight, and environmental impact.

42 ENGINEERING↗
DOE OSTI2020

Deformation of Lie–Poisson algebras and chirality

Linearization of a Hamiltonian system around an equilibrium point yields a set of Hamiltonian symmetric spectra: If λ is an eigenvalue of the linearized generator, -λ and $\barλ$ (hence, -$\barλ$) are also eigenvalues—the former implies a time-reversal symmetry, while the latter guarantees the reality of the solution. However, linearization around a singular equilibrium point (which commonly exists in noncanonical Hamiltonian systems) works out differently, resulting in breaking of the Hamiltonian symmetry of spectra; time-reversal asymmetry causes chirality. This interesting phenomenon was first found in analyzing the chiral motion of the rattleback, a boat-shaped top having misaligned axes of inertia and geometry [Z. Yoshida et al., Phys. Lett. A 381, 2772–2777 (2017)]. To elucidate how chiral spectra are generated, we study the three-dimensional Lie–Poisson systems and classify the prototypes of singularities that cause symmetry breaking. The central idea is the deformation of the underlying Lie algebra; invoking Bianchi’s list of all three-dimensional Lie algebras, we show that the so-called class-B algebras, which are produced by asymmetric deformations of the simple algebra s o ( 3 ) , yield chiral spectra when linearized around their singularities. The theory of deformation is generalized to higher dimensions, including the infinite-dimensional Poisson manifolds relevant to fluid mechanics.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

The kinetic origin of the fluid helicity—A symmetry in the kinetic phase space

Helicity, a topological degree that measures the winding and linking of vortex lines, is preserved by ideal (barotropic) fluid dynamics. In the context of the Hamiltonian description, the helicity is a Casimir invariant characterizing a foliation of the associated Poisson manifold. Casimir invariants are special invariants that depend on the Poisson bracket, not on the particular choice of the Hamiltonian. The total mass (or particle number) is another Casimir invariant, whose invariance guarantees the mass (particle) conservation (independent of any specific choice of the Hamiltonian). In a kinetic description (e.g., that of the Vlasov equation), the helicity is no longer an invariant (although the total mass remains a Casimir of the Vlasov’s Poisson algebra). The implication is that some “kinetic effect” can violate the constancy of the helicity. To elucidate how the helicity constraint emerges or submerges, we examine the fluid reduction of the Vlasov system; the fluid (macroscopic) system is a “sub-algebra” of the kinetic (microscopic) Vlasov system. In the Vlasov system, the helicity can be conserved if a special helicity symmetry condition holds. To put it another way, breaking helicity symmetry induces a change in the helicity. We delineate the geometrical meaning of helicity symmetry and show that for a special class of flows (the so-called epi-two-dimensional flows), the helicity symmetry is written as ∂γ = 0 for a coordinate γ of the configuration space.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

The role of breakup and coalescence in fine-scale bubble-induced turbulence. I. Dynamics

We study the effect of bubble breakup and coalescence on fine-scale dynamics of bubbly turbulent flows using direct numerical simulations. We perform two different simulations of dilute bubbly flows of void fraction 0.5%: one with bubbles breaking up and coalescing and the other without these physical processes. The volume of the fluid method is used for simulating bubbles undergoing breakup and coalescence while the bubbles are treated as rigid spheres in the immersed boundary method simulation. The energy spectrum in both types of simulation, consistent with previous studies, exhibits a -3 slope. We follow a single infinitesimal fluid element as it evolves to understand velocity gradient dynamics using conditional mean trajectories. We note finite-time divergence when the fluid element evolves under the action of inertial and pressure dynamics. The inertial, pressure, and viscous velocity gradient dynamics, when considered individually, produce the same results for bubble-induced turbulence (BIT) as with the classical homogeneous isotropic turbulence (HIT). Yet when the overall velocity gradient dynamics is considered, BIT results in non-cyclic trajectories moving toward stable node and unstable saddle while classical HIT shows cyclic behavior in their trajectories that move toward the origin. Interestingly, both the volume of fluid and immersed boundary simulations produce similar results. Therefore, there are two main takeaways from this research. First, new velocity gradient models are needed for BIT as their velocity gradients behavior is entirely different from the HIT. Second, we can neglect the bubble topology, breakup, and coalescence while studying or modeling the fine-scale dynamics of BIT.

42 ENGINEERING↗
DOE OSTI2023

Helfrich-Hurault elastic instabilities driven by geometrical frustration

Here, the Helfrich-Hurault (HH) elastic instability is a well-known mechanism behind patterns that form as a result of strain upon liquid crystal systems with periodic ground states. In the HH model, layered structures undulate and buckle in response to local, geometric incompatibilities in order to maintain the preferred layer spacing. Classic HH systems include cholesteric liquid crystals under electromagnetic field distortions and smectic liquid crystals under mechanical strains, where both materials are confined between rigid substrates. However, richer phenomena are observed when undulation instabilities occur in the presence of deformable interfaces and variable boundary conditions. Understanding how the HH instability is affected by deformable surfaces is imperative for applying the instability to a broader range of materials. In this review, the HH mechanism is reexamined and special focus is given to how the boundary conditions influence the response of lamellar systems to geometrical frustration. Lamellar liquid crystals confined within a spherical shell geometry are used as the model system. Made possible by the relatively recent advances in microfluidics within the past 15 years, liquid crystal shells are composed entirely of fluid interfaces and have boundary conditions that can be dynamically controlled at will. Past and recent work that exemplifies how topological constraints, molecular anchoring conditions, and boundary curvature can trigger the HH mechanism in liquid crystals with periodic ground states is examined. The review ends by identifying similar phenomena across a wide variety of materials, both biological and synthetic. The fact that the HH mechanism is a generic and often overlooked response of periodic materials to geometrical frustration is highlighted.

42 ENGINEERING↗
DOE OSTI2025

Wave topology in Hall magnetohydrodynamics

Hall magnetohydrodynamics (HMHD) extends ideal MHD by incorporating the Hall effect via the induction equation, making it more accurate for describing plasma behavior at length scales below the ion skin depth. Despite its importance, a comprehensive description of the eigenmodes in HMHD has been lacking. In this work, we derive the complete spectrum and eigenvectors of HMHD waves and identify their underlying topological structure. We prove that the HMHD wave spectrum is homotopic to that of ideal MHD, consisting of three distinct branches: the slow magnetosonic-Hall waves, the shear Alfvén-Hall waves, and the fast magnetosonic-Hall waves, which continuously reduce to their ideal MHD counterparts in the limit of vanishing Hall parameter. Contrary to a recent claim [Mahajan, Sharma, and Lingam, Phys. Plasmas 31, 090701 (2024)], we find that HMHD does not admit any additional wave branches beyond those in ideal MHD. In conclusion, the key qualitative difference lies in the topological nature of the HMHD wave structure: it exhibits nontrivial topology characterized by a Weyl point—an isolated eigenmode degeneracy point—and associated nonzero Chern numbers of the eigenmode bundles over a 2-sphere in 𝐤-space surrounding the Weyl point.

Alfvén waves↗
DOE OSTI2025

Using principal component analysis to distinguish different dynamic phases in superconducting vortex matter

Vortices in type-II superconductors driven over random disorder are known to exhibit a remarkable variety of distinct nonequilibrium dynamical phases that arise owing to the competition between vortex-vortex interactions, the quenched disorder, and the drive. These include pinned states, elastic flows, plastic or disordered flows, and dynamically reordered moving crystal or moving smectic states. The plastic flow phases can be particularly difficult to characterize since the flows are strongly disordered. Here, we perform principal component analysis (PCA) on the positions and velocities of vortex matter moving over random disorder for different disorder strengths and drives. We find that PCA can distinguish the known dynamic phases as well as or better than previous measures based on transport signatures or topological defect densities. In addition, PCA recognizes distinct plastic flow regimes, a slowly changing channel flow and a moving amorphous fluid flow, that do not produce distinct signatures in the standard measurements. In conclusion, our results suggest that this position and velocity-based PCA approach could be used to characterize dynamic phases in a broader class of systems that exhibit depinning and nonequilibrium phase transitions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Nematic Colloidal Micro‐Robots as Physically Intelligent Systems

Abstract Physically intelligent micro‐robotic systems exploit information embedded in micro‐robots, their colloidal cargo, and their milieu to interact, assemble, and form functional structures. Nonlinear anisotropic fluids such as nematic liquid crystals (NLCs) provide untapped opportunities to embed interactions via their topological defects, complex elastic responses, and ability to dramatically restructure in dynamic settings. Here a four‐armed ferromagnetic micro‐robot is designed and fabricated to embed and dynamically reconfigure information in the nematic director field, generating a suite of physical interactions for cargo manipulation. The micro‐robot shape and surface chemistry are designed to generate a nemato‐elastic energy landscape in the domain that defines multiple modes of emergent, bottom‐up interactions with passive colloids. Micro‐robot rotation expands the ability to sculpt interactions; the energy landscape around a rotating micro‐robot is dynamically reconfigured by complex far‐from‐equilibrium dynamics of the micro‐robot's companion topological defect. These defect dynamics allow transient information to be programmed into the domain and exploited. Robust micro‐robotic manipulation strategies are demonstrated that exploit these diverse modes of nemato‐elastic interaction to achieve cargo docking, transport, release, and assembly of complex reconfigurable structures at multi‐stable sites. Such structures are of great interest to future developments of LC‐based advanced optical device and micro‐manufacturing in anisotropic environments.

Yao, Tianyi↗
DOE OSTI2024

ExaWind: Open‐source CFD for hybrid‐RANS/LES geometry‐resolved wind turbine simulations in atmospheric flows

Abstract Predictive high‐fidelity modeling of wind turbines with computational fluid dynamics, wherein turbine geometry is resolved in an atmospheric boundary layer, is important to understanding complex flow accounting for design strategies and operational phenomena such as blade erosion, pitch‐control, stall/vortex‐induced vibrations, and aftermarket add‐ons. The biggest challenge with high‐fidelity modeling is the realization of numerical algorithms that can capture the relevant physics in detail through effective use of high‐performance computing. For modern supercomputers, that means relying on GPUs for acceleration. In this paper, we present ExaWind, a GPU‐enabled open‐source incompressible‐flow hybrid‐computational fluid dynamics framework, comprising the near‐body unstructured grid solver Nalu‐Wind, and the off‐body block‐structured‐grid solver AMR‐Wind, which are coupled using the Topology Independent Overset Grid Assembler. Turbine simulations employ either a pure Reynolds‐averaged Navier–Stokes turbulence model or hybrid turbulence modeling wherein Reynolds‐averaged Navier–Stokes is used for near‐body flow and large eddy simulation is used for off‐body flow. Being two‐way coupled through overset grids, the two solvers enable simulation of flows across a huge range of length scales, for example, 10 orders of magnitude going from O(μm) boundary layers along the blades to O(10 km) across a wind farm. In this paper, we describe the numerical algorithms for geometry‐resolved turbine simulations in atmospheric boundary layers using ExaWind. We present verification studies using canonical flow problems. Validation studies are presented using megawatt‐scale turbines established in literature. Additionally presented are demonstration simulations of a small wind farm under atmospheric inflow with different stability states.

17 WIND ENERGY↗
DOE OSTI2020

Control of the Folding Dynamics of Self-Reconfiguring Magnetic Microbots Using Liquid Crystallinity

Reconfigurable microdevices are being explored in a range of contexts where their life-like abilities to move and change shape are important. While much work has been done to control the motion of these microdevices by engineering their geometry and composition, little is known about their dynamics in complex fluid environments with non-Newtonian rheology. Here in this paper, it is shown how the actuation dynamics of reconfigurable microdevices made by assembly of patchy magnetic microcubes, which are referred to as “microbots,” can be modulated by their interactions with the anisotropic viscoelastic environment of a liquid crystal (LC). The free energy arising from the elastic strain of LC and formation of topological defects around the microbots influences their folding dynamics, which can be tuned by tailoring both the far-field orientation of the LC and the local ordering of the LC at the microbot surfaces. These findings represent a first step toward establishing a general set of design rules to control the dynamics of microbots using complex anisotropic fluids.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Influence of Heave Plate Topology on Reaction Force

In this study, we explore the role of heave plate topology on fluid reaction forces using three, representative shapes: a hexagonal flat plate, a hexagonal conic with an open top, and the same with a closed top that encloses a fluid mass. We force each test article sinusoidally in a quiescent tank and decompose the reaction force using forms of the Morison equation for phase-invariant and phase-dependent parameterizations. We find that a flat plate generates 5.3 % more fluid reaction force than the open conic topology, and 21.4% more than the enclosed conic. Asymmetric topologies generate asymmetric reaction forces, but the magnitude of asymmetry is limited by nearly symmetric fluid inertia forces, which dominate over drag for these test articles. Additionally, we observe asymmetric vortex dynamics for the flat plate when the Keulegan-Carpenter number is between 1 and 2, accompanied by a shift in the phase of the peak force by of the oscillation period. As a consequence of this shift, the hydrodynamic coefficients estimated from the phase-dependent Morison equation decomposition are asymmetric, suggesting that phase-dependent representations may not provide physical insight in some hydrodynamic regimes.

added mass↗