Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fluid structure interaction”

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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Three dimensional simulation of an underexpanded jet interacting with a supersonic cross flow

A new three-dimensional computational fluid dynamics (CFD) code has been developed to simulate the flow fields of an underexpanded jet transversely injected into supersonic air stream inside the combustors of ramjets and scramjets. The code employs an implicit finite volume, lower-upper (LU) time marching method to solve the complete three-dimensional Navier-Stokes equations in a fully-coupled and very efficient manner. Results clearly depict the flow characteristics, including the shock structure, separated flow regions around the injector, and the wake flow in the lee of the injector.

Yu, Sheng-Tao↗
DOE OSTI2020

Fate of Internal Waves on a Shallow Shelf

Internal waves strongly influence the physical and chemical environment of coastal ecosystems worldwide. We report novel observations from a distributed temperature sensing (DTS) system that tracked the transformation of internal waves from the shelf break to the surf zone over a narrow shelf slope region in the South China Sea. The spatially continuous view of temperature fields provides a perspective of physical processes commonly available only in laboratory settings or numerical models, including internal wave reflection off a natural slope, shoreward transport of dense fluid within trapped cores, and observations of internal rundown (near-bed, offshore-directed jets of water preceding a breaking internal wave). Analysis shows that the fate of internal waves on this shelf—whether transmitted into shallow waters or reflected back offshore—is mediated by local water column density structure and background currents set by the previous shoaling internal waves, highlighting the importance of wave-wave interactions in nearshore internal wave dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

Relativistic Modeling Capabilities in Perseus Extended MHD Simulation Code for HED Plasmas

This project is submitted in response to DE-FOA-0001153, and will be carried out by Dr. Charles Seyler and Dr. Nathaniel Hamlin as post doc under the direction of Professor Seyler. The objectives are to enable PERSEUS to self-consistently model relativistic high-energy-density (HED) laboratory and astrophysical phenomena. Electrons and ions in a plasma are generally modeled either as fluids using a two-fluid code, like PERSEUS, or as particles using a Particle-In- Cell (PIC) code. Nonrelativistic PERSEUS solves the two-fluid equations, formulated in terms of a generalized Ohm’s law (GOL), so as to model about nine orders of magnitude in density variation using a local semi-implicit method. Relativistic PERSEUS preserves this structure and therefore retains these advantageous properties, enabling it to model a broader range of relativistic HED phenomena than previous codes. We have also overcome a major technical challenge associated with solving a relativistic system of equations relating conserved quantities (momentum, energy, etc.) to primitive variables (density, velocity, and pressure). We have thus far developed a relativistic version of PERSEUS and demonstrated that it recovers expected nonrelativistic results. We are currently using relativistic PERSEUS to model two HED phenomena, namely laser-plasma interactions and X-pinches, both of which show spectroscopic evidence of relativistic electrons. In an X-pinch, a section of wire carrying sufficiently large current forms a plasma column that pinches down to a point-like source, generating X-rays for use in radiographic imaging. Laser-plasma interactions have been simulated using PIC codes, and we use their results as a benchmark for comparison. However, neither PIC codes nor existing two-fluid schemes can simulate X-pinches, due to the broad range of densities which have thus far made such a simulation computationally prohibitive. In both laser-plasma interaction and Xpinch simulations, we have observed relativistic phenomena that PERSEUS could not previously model. This includes relativistic channeling of a laser into an over-dense deuterium gas (i.e. induced transparency), consistent with PIC simulation results, along with energetic electrons and ions following an X-pinch and likely accompanied by strong X-ray emission, consistent with experimental observations of energetic X-rays. This proposal includes our latest simulation results. We will improve computational efficiency with techniques that allow the use of a less refined numerical grid. This will enable extensions to three-dimensional geometries. We will improve modeling of laser-plasma interactions and X-pinches, including radiative transport processes, in order to improve comparisons with PIC simulations and experimental diagnostics. We will also model power loss through electrode surface plasmas in magnetically insulated transmission lines (MITLs). We will model relativistic astrophysical phenomena in both astrophysical and HED regimes, which could shed light on the validity of laboratory astrophysics for certain problems. These phenomena are currently modeled using PIC codes and existing relativistic two-fluid codes, neither of which can model the range of length scales of which relativistic PERSEUS is capable. There is strong potential for this project to have a broad impact on both the astrophysical and laboratory plasma communities. In the laboratory community, our work is applicable to those concerned with laser-plasma interactions such as laser-driven fusion, and those in the dense Z pinch community concerned with X-pinches and electrode surface plasmas. An improved understanding of the mechanisms of hard (i.e. high-energy) X-ray generation will likely make important contributions toward the use of X-pinches for radiography through the harnessing or mitigation of hard X-rays, depending on the application. The modeling of electrode surface plasmas, including vacuum electron flow, will improve our understanding of power loss to various loads.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Turbulence theories and statistical closure approaches

When discussing research in physics and in science more generally, it is common to ascribe equal importance to the three components of the scientific trinity: theoretical, experimental, and computational studies. This review will explore the future of modern turbulence theory by tracing its history, which began in earnest with Kolmogorov’s 1941 analysis of turbulence cascade and inertial range [A.N. Kolmogorov, Dokl. Akad. Nauk SSSR, 30, 299, (1941); 32, 19, (1941)]. The 80th Anniversary of Kolmogorov’s landmark study is a welcome opportunity to survey the achievements and evaluate the future of the theoretical approach of turbulence research. Over the years, turbulence theories have been critically important in laying the foundation of our understanding of the nature of turbulent flows. In particular, the Direct Interaction Approximation (DIA) [R.H. Kraichnan, J. Fluid Mech., 5, 497 (1959)] and its subsequent development, known as the statistical closure approach, can be identified as perhaps the most profound single advancement. The remarkable success of the statistical closure has furnished a platform to study such essential concepts as the energy transfer process and interacting scales, and the roles of the straining and sweeping motions. More recently, the quasi-Lagrangian formulation of V. L’vov & I. Procaccia and Kraichnan’s solvable passive scalar model provided powerful ways to explore another fundamental aspect of turbulent flows, the phenomena of intermittency, and the associated anomalous scaling exponents. In the meantime, the theory of fluid equilibria has been developed to describe the large-scale structures that can emerge from turbulent cascades of two-dimensional and geophysical flows at a later time. And yet, despite all these successes, analytical treatments suffer from mathematical complexities. As a result, the utility of theoretical approaches has been limited to relatively idealized flows. On the other hand, in recent decades, computational abilities and experimental facilities have reached an unprecedented scale. Looking beyond the horizon, the imminent deployment of exascale supercomputers will generate complete datasets of the entire flow field of key benchmark flows, allowing researchers to extract additional measurements concerning fully developed, complex turbulent flow fields far beyond those available from the statistical closure theories. Some other developments that could potentially influence the future course of turbulence theories include the advancement of machine learning, artificial intelligence, and data science; likely disruptions arising from the advent of quantum computation; and the increasingly prominent role of turbulence research in providing more accurate climate scientific data. Finally, turbulence theorists can leverage these developments by asking the right questions and developing advanced, sophisticated frameworks that will be able to predict and correlate vast amounts of data from the other two components of the trinity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Dynamic properties of different liquid states in systems with competing interactions studied with lysozyme solutions

Recent studies of colloidal systems with a short-range attraction and long-range repulsion (SALR) have been demonstrated to have a generalized phase diagram with multiple liquid states defined by their structures. In this paper, we identify the different liquid states of previous experimentally studied lysozyme samples within this proposed generalized state diagram and explore the dynamic properties of each liquid state. We show that most lysozyme samples studied here and previously at low and intermediate concentrations are dispersed fluids while a few high concentration samples are randomly percolated liquids. In the dispersed fluid region, the short-time diffusion coefficient measured by neutron spin echo agrees well with the long time diffusion coefficient estimated with the solution viscosity. This dynamic feature is maintained even for some samples in the random percolated region. However, the short-time and long-time diffusion coefficients of random percolated fluids deviate at larger concentration and attraction strength. At high enough concentrations, the mean square displacement can be as slow as those of many glassy colloidal systems at time scales near the characteristic diffusion time even though these lysozyme samples remain in liquid states at the long-time limit. Here, we thus identify the region in the generalized phase diagram where these equilibrium states with extremely slow local dynamics exist relative to bulk percolation and kinetic arrest (gel and glassy) transitions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Unsupervised discovery of extreme weather events using universal representations of emergent organization

Spontaneous self-organization is ubiquitous in systems far from thermodynamic equilibrium. While organized structures that emerge dominate transport properties, universal representations that identify and describe these key objects remain elusive. Here, we introduce a theoretically grounded framework for describing emergent organization that, via data-driven algorithms, is constructive in practice. Its building blocks are spacetime lightcones that embody how information propagates across a system through local interactions. We show that predictive equivalence classes of lightcones—local causal states—capture organized behaviors in complex spatiotemporal systems. Employing an unsupervised physics-informed machine learning algorithm and a high-performance computing implementation, we demonstrate automatically discovering organized structures in two real-world domain science problems. We show that local causal states identify vortices and track their power-law decay behavior in two-dimensional fluid turbulence. We then show how to detect and track familiar extreme weather events—hurricanes and atmospheric rivers—and discover other novel structures associated with precipitation extremes in high-resolution climate data at the grid-cell level.

Rupe, Adam [Pacific Northwest National Laboratory ↗
DOE OSTI2015

Structure of retrograde circumbinary accretion disks

We evaluate the interactions between a central equal mass binary and a surrounding retrograde circumbinary accretion disk in a three-dimensional, MHD simulation. We find (as widely expected) that the net torque expressed on a retrograde circumbinary disk by the central binary is much (more than two orders of magnitude) smaller than when the disk orbits in the prograde direction. For this reason, unlike the prograde case, there is no “gap” carved out of the retrograde disk, and accreting fluid can travel directly from the disk to the domain of the binary. For the same reason, retrograde disks do not develop the large-amplitude “lumps” generically seen near the inner edges of prograde disks. Although not directly treated in this simulation, the finite extent of accretion disks around the individual masses of the binary can lead to shocks between their outer edges and the circumbinary disk, potentially creating a photon signal modulated at the binary orbital period. The spectra of these shocks should peak in the hard X-ray band as is expected for prograde systems, but retrograde systems would not be expected to display a “notch” in the optical/UV band. Because the disk’s angular momentum is directed oppositely to that of the binary, accretion onto the binary causes a rapid increase in eccentricity. Lastly, we discuss the comparative rates of orbital evolution associated with accretion in prograde and retrograde circumbinary disks, finding that the principal contrast is that e-dot should be rather greater in the retrograde case.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Going beyond BEM with BEM: an insight into dynamic inflow effects on floating wind turbines

Blade element momentum (BEM) theory is the backbone of many industry-standard wind turbine aerodynamic models. To be applied to a broader set of engineering problems, BEM models have been extended since their inception and now include several empirical corrections. These models have benefitted from decades of development and refinement and have been extensively used and validated, proving their adequacy in predicting aerodynamic forces of horizontal-axis wind turbine rotors in most scenarios. However, the analysis of floating offshore wind turbines (FOWTs) introduces new sets of challenges, especially if new-generation large and flexible machines are considered. In fact, due to the combined action of wind and waves and their interaction with the turbine structure and control system, these machines are subject to unsteady motion and thus unsteady inflow on the wind turbine's blades, which could put BEM models to the test. Consensus has not been reached on the accuracy limits of BEM in these conditions. This study contributes to the ongoing research on the topic by systematically comparing four different aerodynamic models, ranging from BEM to computational fluid dynamics, in an attempt to shed light on the unsteady aerodynamic phenomena that are at stake in FOWTs and whether BEM is able to model them appropriately. Simulations are performed on the UNAFLOW 1:75 scale rotor during imposed harmonic surge and pitch motion. Experimental results are available for these conditions and are used for baseline validation. The rotor is analyzed in both rated operating conditions and low wind speeds, where unsteady aerodynamic effects are expected to be more pronounced. Results show that BEM, despite its simplicity, can adequately model the aerodynamics of FOWTs in most conditions if augmented with a dynamic inflow model.

17 WIND ENERGY↗
NASA NTRS1996

Relativistic Navigation: A Theoretical Foundation

We present a theoretical foundation for relativistic astronomical measurements in curved space-time. In particular, we discuss a new iterative approach for describing the dynamics of an isolated astronomical N-body system in metric theories of gravity. To do this, we generalize the Fock-Chandrasekhar method of the weak-field and slow-motion approximation (WFSMA) and develop a theory of relativistic reference frames (RF's) for a gravitationally bounded many-extended-body problem. In any proper RF constructed in the immediate vicinity of an arbitrary body, the N-body solutions of the gravitational field equations are formally presented as a sum of the Riemann-flat inertial space-time, the gravitational field generated by the body itself, the unperturbed solutions for each body in the system transformed to the coordinates of this proper RF, and the gravitational interaction term. We develop the basic concept of a general WFSMA theory of the celestial RF's applicable to a wide class of metric theories of gravity and an arbitrary model of matter distribution. We apply the proposed method to general relativity. Celestial bodies are described using a perfect fluid model; as such, they possess any number of internal mass and current multipole moments that explicitly characterize their internal structures. The obtained relativistic corrections to the geodetic equations of motion arise because of a coupling of the bodies' multiple moments to the surrounding gravitational field. The resulting relativistic transformations between the different RF's extend the Poincare group to the motion of deformable self-gravitating bodies. Within the present accuracy of astronomical measurements we discuss the properties of the Fermi-normal-like proper RF that is defined in the immediate vicinity of the extended compact bodies. We further generalize the proposed approximation method and include two Eddington parameters (gamma, Beta). This generalized approach was used to derive the relativistic equations of satellite motion in the vicinity of the extended bodies. Anticipating improvements in radio and laser tracking technologies over the next few decades, we apply this method to spacecraft orbit determination. We emphasize the number of feasible relativistic gravity tests that may be performed within the context of the parameterized WFSMA. Based on the planeto-centric equations of motion of a spacecraft around the planet, we suggested a new null test of the Strong Equivalence Principle (SEP). The experiment to measure the corresponding SEP violation effect could be performed with the future Mercury Orbiter mission. We discuss other relativistic effects, including the perihelion advance and the redshift and geodetic precession of the orbiter's orbital plane about Mercury, as well as the possible future implementation of the proposed formalism in software codes developed for solar-system orbit determination. All the important calculations are completely documented, and the references contain an extensive list of cited literature.

Turyshev, Slava G.↗
DOE OSTI2019

LBNF Hadron Absorber: Updated Mechanical Design and Analysis for 2.4 MW Operation

The Long-Baseline Neutrino Facility (LBNF) Hadron Absorber is located downstream of the decay pipe. It consists of actively cooled aluminum and steel blocks surrounded by steel and concrete shielding. Majority of the beam power is deposited in the absorber core which is water cooled. The surrounding steel and concrete shielding are air-cooled. The absorber provides radiation protection to personnel and keeps soil and ground activation levels below allowable limits. It is designed for 2.4 MW beam operations. The total heat load deposited into the absorber is approximately 400 kW. The current design considers the longer 4-interaction length target of the optimized beam design. In addition, the ‘bafflette’ around the target reduces the energy deposited into the absorber. For this reason, the sculpting in the aluminum core blocks, which was in the previous design, was removed, making the design uniform and less complicated. In addition, the uniformity of the absorber makes it easier to understand the muon monitor data. Steady state thermal, structural, and Computational Fluid Dynamics (CFD) analysis of critical absorber aluminum and steel components during steady state operations is discussed herein. A similar analysis for a 120 GeV, 10 μs pulse, accident condition is also discussed. A preliminary design for the accident pulse prevention system that protects the absorber is also described.

43 PARTICLE ACCELERATORS↗
NASA NTRS1993

Double simple-harmonic-oscillator formulation of the thermal equilibrium of a fluid interacting with a coherent source of phonons

A formulation is given for a collection of phonons (sound) in a fluid at a non-zero temperature which uses the simple harmonic oscillator twice; one to give a stochastic thermal 'noise' process and the other which generates a coherent Glauber state of phonons. Simple thermodynamic observables are calculated and the acoustic two point function, 'contrast' is presented. The role of 'coherence' in an equilibrium system is clarified by these results and the simple harmonic oscillator is a key structure in both the formulation and the calculations.

Defacio, B.↗
DOE OSTI2023

Wetting transition and fluid trapping in a microfluidic fracture

Immiscible fluid–fluid displacement in confined geometries is a fundamental process occurring in many natural phenomena and technological applications, from geological CO 2 sequestration to microfluidics. Due to the interactions between the fluids and the solid walls, fluid invasion undergoes a wetting transition from complete displacement at low displacement rates to leaving a film of the defending fluid on the confining surfaces at high displacement rates. While most real surfaces are rough, fundamental questions remain about the type of fluid–fluid displacement that can emerge in a confined, rough geometry. Here, we study immiscible displacement in a microfluidic device with a precisely controlled structured surface as an analogue for a rough fracture. We analyze the influence of the degree of surface roughness on the wetting transition and the formation of thin films of the defending liquid. We show experimentally, and rationalize theoretically, that roughness affects both the stability and dewetting dynamics of thin films, leading to distinct late-time morphologies of the undisplaced (trapped) fluid. Finally, we discuss the implications of our observations for geologic and technological applications.

42 ENGINEERING↗
DOE OSTI2020

Solvation Energy of Ions in a Stockmayer Fluid

We calculate the solvation energy of monovalent and divalent ions in various liquids with coarse-grained molecular dynamics simulations. Our theory treats the solvent as a Stockmayer fluid, which accounts for the intrinsic dipole moment of molecules and the rotational dynamics of the dipoles. Despite the simplicity of the model, we obtain qualitative agreement between the simulations and experimental data for the free energy and enthalpy of ion solvation, which indicates that the primary contribution to the solvation energy arises mainly from the first and possibly second solvation shells near the ions. Our results suggest that a Stockmayer fluid can serve as a reference model that enables direct comparison between theory and experiment and may be invoked to scale up electrostatic interactions from the atomic to the molecular length scale.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2012

Molecular Structures of Fluid Phase Phosphatidylglycerol Bilayers as Determined by Small Angle Neutron and X-ray Scattering

We have determined the molecular structures of commonly used phosphatidylglycerols (PGs) in the commonly accepted biologically relevant fluid phase. This was done by simultaneously analyzing small-angle neutron and X-ray scattering data, with the constraint of a measured volume per lipid. We report the temperature dependence of bilayer parameters obtained using the one dimensional scattering density profile model derived from molecular dynamics simulations, including the area per lipid, the overall bilayer thickness, as well as other intrabilayer parameters (e.g., hydrocarbon thickness). Lipid areas are found to be larger than their phosphatidylcholine (PC) counterparts, a result likely due to repulsive electrostatic interactions taking place between the charged PG headgroups even in the presence of Na+ counterions. In general, PG and PC bilayers show a similar response to changes in temperature and chain length, but differ in their response to chain unsaturation. For example, compared to their PC counterparts, the inclusion of a double bond in PG lipids results in a smaller incremental change in lipid area and bilayer thickness. However, the extrapolated lipid area of saturated PG lipids at infinite chain length is found to be similar to that of PCs, an indication of the glycerol-carbonyl backbone's pivotal role in influencing the lipid-water interface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2013

THE ROLE OF SUPERLUMINAL ELECTROMAGNETIC WAVES IN PULSAR WIND TERMINATION SHOCKS

The dynamics of a standing shock front in a Poynting-flux-dominated relativistic flow is investigated by using a one-dimensional, relativistic, two-fluid simulation. An upstream flow containing a circularly polarized, sinusoidal magnetic shear wave is considered, mimicking a wave driven by an obliquely rotating pulsar. It is demonstrated that this wave is converted into large-amplitude electromagnetic waves with superluminal phase speeds by interacting with the shock when the shock-frame frequency of the wave exceeds the proper plasma frequency. The superluminal waves propagate in the upstream, modify the shock structure substantially, and form a well-developed precursor region ahead of a subshock. Dissipation of Poynting flux occurs in the precursor as well as in the downstream region through a parametric instability driven by the superluminal waves. The Poynting flux remaining in the downstream region is carried entirely by the superluminal waves. The downstream plasma is therefore an essentially unmagnetized, relativistically hot plasma with a non-relativistic flow speed, as suggested by observations of pulsar wind nebulae.

79 ASTRONOMY AND ASTROPHYSICS↗