Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrodynamics modeling”

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 1,621 records · Page 90

A review of modelling techniques for floating offshore wind turbines

Modelling floating offshore wind turbines (FOWTs) is challenging due to the strong coupling between the aerodynamics of the turbine and the hydrodynamics of the floating platform. Physical testing at scale is faced with the additional challenge of the scaling mismatch between Froude number and Reynolds number due to working in the two fluid domains, air and water. In the drive for cost-reduction of floating wind energy, designers may be seeking to move towards high-fidelity numerical modelling as a substitute for physical testing. However, the numerical engineering tools typically used for FOWT modelling are considered as mid-fidelity to low-fidelity tools, and currently lack the level of accuracy required to do so. Furthermore, there is a lack of operational FOWT data available for further development and validation. High-fidelity tools, such as CFD, have greater accuracy but are cumbersome tools and still require validation. Physical scale model testing therefore continues to play an essential role in the development of FOWTs both as a source of validation data for numerical models and as an important development step along the path to commercialization of all platform concepts. The aim of this paper is to provide an overview of both numerical modelling and physical FOWT scale model testing approaches and to provide guidance on the selection of the most appropriate approach (or combination of approaches). The current state-of-the-art will be discussed along with current research trends and areas for further investigation.

17 WIND ENERGY↗

Xenon Fractionation, Hydrogen Escape, and the Oxidation of the Earth

Xenon in Earth's atmosphere is severely mass fractionated and depleted compared to any plausible solar system source material, yet Kr is unfractionated. These observations seem to imply that Xe has escaped from Earth. But to date no process has been identified that can cause Xe, which is heavier than Kr, to escape while Kr does not. Vigorous hydrodynamic hydrogen escape can produce mass fractionation in heavy gases. The required hydrogen flux is very high but within the possible range permitted by solar Extreme Ultraviolet radiation (EUV, which here means radiation at wavelengths short enough to be absorbed efficiently by hydrogen) heating when Earth was on the order of 100 Myrs old or younger. However this model cannot explain why Xe escapes but Kr does not. Recently, what appears to be ancient atmospheric xenon has been recovered from several very ancient (3-3.5 Ga) terrestrial hydrothermal barites and cherts. What is eye-catching about this ancient Xe is that it is less fractionated that Xe in modern air. In other words, it appears that a process was active on Earth some 3 to 3.5 billion years ago that caused xenon to fractionate. By this time the Sun was no longer the EUV source that it used to. If xenon was being fractionated by escape — currently the only viable hypothesis — it had to be in the less unfamiliar context of Earth’s Archean atmosphere and under rather modest levels of EUV forcing. This requires a new model. Here we address the circumstances in which Xe, but not Kr, could escape from Earth as an ion. In a hydrodynamically escaping hydrogen wind the hydrogen is partially photo-ionized. The key concepts are that ions are much more strongly coupled to the escaping flow than are neutrals (so that a relatively modest flow of H and H+ to space could carry Xe+ along with it), and that xenon alone among the noble gases is more easily ionized than hydrogen. This sort of escape is possible if not prevented by a planetary magnetic field. The best prospects for Earth are therefore escape along the polar field lines, although a very weak or absent magnetic field would likely work as well. As applied to the Archean Earth the discussion will be constrained by diffusion-limited hydrogen escape. The extended history of hydrogen escape implicit in Xe escape in the Archean is consistent with suggestions that hydrogen escape from the anoxic Archean atmosphere was considerable, because biogenic methane is expected to have been rather abundant. Hydrogen escape plausibly played the key role in creating oxidizing condition at the surface of the Earth and setting the stage for the creation of an O2 atmosphere.

photo-ionized gases↗

Theoretical spectra of circumstellar dust shells around carbon-rich asymptotic giant branch stars

Realistic modeling of circumstellar dust shells around evolved stars has to include a physical treatment of the interaction among hydrodynamics, thermodynamics, radiative transfer, chemistry and dust formation and -growth. A self-consistent solution of this problem is presented in the case of stationary, spherical symmetric dust-driven winds. The resulting shell structure and the mass-loss rate are completely determined by the three fundamental stellar parameters stellar mass M(stellar), stellar luminosity L(stellar) and effective temperature T(sub eff) and by the element abundances epsilon(sub i). A detailed calculation of the transport coefficients of the dust component by means of the particle size distribution function and the solution of the non-grey radiative transfer problem provide realistic synthetic spectra of the dust shell models. We discuss the dependence of the resulting spectra on the stellar parameters in terms of infrared two color diagrams. Application of these model calculations to the prominent infrared object IRC +10216 yields best agreement with the observed spectrum and the visibility data at maximum light for the stellar parameters M(stellar) = 0.7 solar mass, L(stellar) = 2.4 x 10(exp 4) solar luminosity, T(stellar) = 2010K and a carbon to oxygen ratio of epsilon(sub c)/epsilon(sub o) = 1.40, which corresponds to a mass-loss rate of M-dot = 8 x 10(exp -5) solar mass/yr. In this model only amorphous carbon grains are considered as the main opacity source. From this model a distance to IRC +10216 of d = 170pc is deduced. The total mass contained in the circumstellar dust shell implies and initial main sequence mass of M(sub ZAMS) greater than or = 1.3 solar mass for IRC +10216.

Winters, J. M.↗

Experimental and numerical studies of seismic fluid-structure interaction in a base-supported cylindrical vessel

Seismic design and qualification of advanced reactors will rely heavily on the use of verified and validated numerical models capable of capturing the interaction of the vessel, its contained fluid, and the internal equipment: fluid-structure interaction (FSI) analysis. Analytical solutions can be used for preliminary sizing and design of such vessels but their application is limited to simple geometries and boundary conditions, and small amplitude, translational (and rotational) inputs. To validate numerical models for seismic FSI analysis in finite element codes, a comprehensive set of experiments was performed on a liquid-filled cylindrical vessel, using a 6 degree-of-freedom earthquake simulator. Results in terms of sloshing frequency, damping ratio in sloshing modes, and hydrodynamic responses (wave height, hydrodynamic pressure, base shear, and base moment) for multi-directional earthquake simulator inputs are reported and compared with analytical solutions for liquid-filled vessels. The impact of seismic (base) isolation on hydrodynamic responses was studied using earthquake simulator inputs generated using a virtual isolation system. Data from the experiments are used to validate a numerical model of the fluid-filled vessel using the Arbitrary Lagrangian Eulerian (ALE) solver in the commercial finite element program LS-DYNA. Validation studies are presented for multi-directional seismic inputs, including rocking motions. Lagrangian modeling approaches using an elastic material formulation for the fluid are also investigated and their limitations and possible applications are identified. Here, the results are broadly applicable to the seismic response of base supported, liquid-filled vessels.

58 GEOSCIENCES↗

Trailing edge noise prediction from measured surface pressures

Trailing edge (TE) noise is investigated for the case of a two-dimensional airfoil embedded in a uniform low Mach number flow, and the usefulness of several TE noise theories is examined by applying them to the measured data. The TE noise spectra and directivity are quantitatively determined for the case of a high Reynolds number and a fully turbulent boundary layer. Parameters include angle of attack, flow velocity and TE bluntness. Evanescent wave theories by Chase (1975) and Chandiramani (1974) are compared to the present results and show good agreement. Agreement of the near field pressure scatter phenomenon analysis with measurements implies that the basic assumptions used in the analysis are correct, i.e., the turbulent boundary layer (TBL) flow passes the trailing edge into the wake region. No hydrodynamic wake shedding activity is confirmed for the two-sided TBL flow, and a method incorporating the principles of the coherent output power method is used to determine the sound field. The near field edge scatter model is found to not only establish optimum sizing of edge treatment for noise control, but also to separate and identify the scattered field from the incident hydrodynamic field.

Brooks, T. F.↗

Water Impact Test and Simulation of a Composite Energy Absorbing Fuselage Section

In March 2002, a 25-ft/s vertical drop test of a composite fuselage section was conducted onto water. The purpose of the test was to obtain experimental data characterizing the structural response of the fuselage section during water impact for comparison with two previous drop tests that were performed onto a rigid surface and soft soil. For the drop test, the fuselage section was configured with ten 100-lb. lead masses, five per side, that were attached to seat rails mounted to the floor. The fuselage section was raised to a height of 10-ft. and dropped vertically into a 15-ft. diameter pool filled to a depth of 3.5-ft. with water. Approximately 70 channels of data were collected during the drop test at a 10-kHz sampling rate. The test data were used to validate crash simulations of the water impact that were developed using the nonlinear, explicit transient dynamic codes, MSC.Dytran and LS-DYNA. The fuselage structure was modeled using shell and solid elements with a Lagrangian mesh, and the water was modeled with both Eulerian and Lagrangian techniques. The fluid-structure interactions were executed using the fast general coupling in MSC.Dytran and the Arbitrary Lagrange-Euler (ALE) coupling in LS-DYNA. Additionally, the smooth particle hydrodynamics (SPH) meshless Lagrangian technique was used in LS-DYNA to represent the fluid. The simulation results were correlated with the test data to validate the modeling approach. Additional simulation studies were performed to determine how changes in mesh density, mesh uniformity, fluid viscosity, and failure strain influence the test-analysis correlation.

Fasanella, Edwin L.↗

The Atmospheric Response to High Nonthermal Electron Beam Fluxes in Solar Flares. I. Modeling the Brightest NUV Footpoints in the X1 Solar Flare of 2014 March 29

The 2014 March 29 X1 solar flare (SOL20140329T17:48) produced bright continuum emission in the far- and near-ultraviolet (NUV) and highly asymmetric chromospheric emission lines, providing long-sought constraints on the heating mechanisms of the lower atmosphere in solar flares. We analyze the continuum and emission line data from the Interface Region Imaging Spectrograph (IRIS) of the brightest flaring magnetic footpoints in this flare. We compare the NUV spectra of the brightest pixels to new radiative-hydrodynamic predictions calculated with the RADYN code using constraints on a nonthermal electron beam inferred from the collisional thick-target modeling of hard X-ray data from Reuven Ramaty High Energy Solar Spectroscopic Imager. We show that the atmospheric response to a high beam flux density satisfactorily achieves the observed continuum brightness in the NUV. The NUV continuum emission in this flare is consistent with hydrogen (Balmer) recombination radiation that originates from low optical depth in a dense chromospheric condensation and from the stationary beam-heated layers just below the condensation. A model producing two flaring regions (a condensation and stationary layers) in the lower atmosphere is also consistent with the asymmetric Fe II chromospheric emission line profiles observed in the impulsive phase.

Kowalski, Adam F.↗

Examining the effects of soil entrainment during nuclear cloud rise on fallout predictions using a multiscale atmospheric modeling framework

Current operational models for nuclear cloud rise over land were developed and validated using observations from shallow-buried or surface detonations, where lofted soil quickly mixed with fission products from the detonation. These models poorly predict fallout from elevated detonations near the fallout-free height of burst (FFHOB), where interactions with the ground are limited and the mixing of fission products and lofted soil is incomplete. Fallout-free is a misnomer at this HOB, as fallout was observed in these cases, but was below the levels of concern, especially off-grounds of the nuclear test site. To correctly characterize and model fallout from detonations near the FFHOB, models must be developed which can capture the stratified nature of the particle and activity-size distributions within the cloud. Previously, it was shown that the Weather Research and Forecasting (WRF) model can accurately simulate nuclear cloud rise for airbursts with little to no ground interactions (Arthur et al., 2021). That work is expanded here by (1) using a radiation-hydrodynamics code to improve the fireball initialization in WRF, (2) further developing an aerosol package from WRF-Chem to simulate lofted soil, and (3) combining the WRF cloud rise simulations with the operational models used at the National Atmospheric Release Advisory Center (NARAC) for fallout modeling. Using this combination of codes, the Upshot-Knothole Grable detonation, which was just below the FFHOB, is simulated from seconds after detonation through cloud rise and fallout, and results are compared to historical test data. Here, the results show improved prediction of dose rate and highlight the need to correctly characterize the entrainment of material into the cloud and the subsequent mixing of fission products with entrained material.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Active Galaxy Unification in the Era of X-Ray Polarimetry

Active galactic nuclei (AGNs), Seyfert galaxies, and quasars are powered by luminous accretion and often accompanied by winds that are powerful enough to affect the AGN mass budget, and whose observational appearance bears an imprint of processes that are happening within the central parsec around the black hole (BH). One example of such a wind is the partially ionized gas responsible for X-ray and UV absorption (warm absorbers). Here, we show that such gas will have a distinct signature when viewed in polarized X-rays. Observations of such polarization can test models for the geometry of the flow and the gas responsible for launching and collimating it. We present calculations that show that the polarization depends on the hydrodynamics of the flow, the quantum mechanics of resonance-line scattering, and the transfer of polarized X-ray light in the highly ionized moving gas. The results emphasize the three-dimensional nature of the wind for modeling spectra. We show that the polarization in the 0.1-10 keV energy range is dominated by the effects of resonance lines. We predict a 5%-25% X-ray polarization signature of type-2 objects in this energy range. These results are generalized to flows that originate from a cold torus-like structure, located approximately 1 pc from the BH, which wraps the BH and is ultimately responsible for the apparent dichotomy between type 1 and type 2 AGNs. Such signals will be detectable by future dedicated X-ray polarimetry space missions, such as the NASA Gravity and Extreme Magnetism Small Explorer SMEX, "GEMS" Swank et al. (2008).

Dorodnitsyn, A.↗

Pre-equilibrium photons from the early stages of heavy-ion collisions

Abstract We use QCD kinetic theory to compute photon production in the chemically equilibrating Quark-Gluon Plasma created in the early stages of high-energy heavy-ion collisions. We do a detailed comparison of pre-equilibrium photon rates to the thermal photon production. We show that the photon spectrum radiated from a hydrodynamic attractor evolution satisfies a simple scaling form in terms of the specific shear viscosityη/sand entropy densitydS/dζ∼ (Tτ 1/3 ) 3/2 ∞ . We confirm the analytical predictions with numerical kinetic theory simulations. We use the extracted scaling function to compute the pre-equilibrium photon contribution in$$ \sqrt{s_{NN}} $$ s NN = 2.76 TeV 0–20% PbPb collisions. We demonstrate that our matching procedure allows for a smooth switching from pre-equilibrium kinetic to thermal hydrodynamic photon production. Finally, our publicly available implementation can be straightforwardly added to existing heavy ion models.

Physics↗

Four-dimensional QCD equation of state with multiple chemical potentials

Here, we construct a four-dimensional version of the equation of state (EoS) model neos, neos-4d, as a function of the temperature and chemical potentials of baryon, electric charge, and strangeness for the hot and dense quantum chromodynamics (QCD) matter created in relativistic nuclear collisions. This EoS enables multiple conserved charge current evolution in a relativistic fluid. Input from lattice QCD simulations and a hadron resonance gas model is considered for constructing the equation of state. We investigate its applicability to the relativistic hydrodynamic description of nuclear collisions and present a method for efficient numerical implementation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Impact of late-time neutrino emission on the diffuse supernova neutrino background

In the absence of high-statistics supernova neutrino measurements, estimates of the diffuse supernova neutrino background (DSNB) hinge on the precision of simulations of core-collapse supernovae. Understanding the cooling phase of protoneutron star (PNS) evolution (≳ 1 s after core bounce) is crucial, since approximately 50% of the energy liberated by neutrinos is emitted during the cooling phase. We model the cooling phase with a hybrid method by combining the neutrino emission predicted by 3D hydrodynamic simulations with several cooling-phase estimates, including a novel two-parameter correlation depending on the final baryonic PNS mass and the time of shock revival. We find that the predicted DSNB event rate at Super-Kamiokande can vary by a factor of ~2–3 depending on the cooling-phase treatment. Here we also find that except for one cooling estimate, the range in predicted DSNB events is largely driven by the uncertainty in the neutrino mean energy. With a good understanding of the late-time neutrino emission, more precise DSNB estimates can be made for the next generation of DSNB searches.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Precision linear shaped charge severance of graphite-epoxy materials

This paper presents Precision Linear Shaped Charge (PLSC) components designed to sever a variety of target materials. Recent data for the severance of graphite-epoxy panels or targets with PLSC's are presented. A brief history of the requirement to originate the development of PLSC's for weapon components at Sandia National Laboratories is presented. The Department of Energy's (DOE) nuclear weapon systems have continually decreased in size. Today's relatively small weapons require the design of much more efficient, lighter, and smaller explosive components because fragments, air shocks, and pyro-shocks associated with the function of these components can damage electrical and other sensitive components located nearby. The DOE requirements for PLSC's are listed. Therefore, linear shaped charge (LSC) components for weapon systems can no longer be empirically or experimentally designed for a given application. Many of today's designs require severing concentric cylinders, for example, where the LSC jet is designed to sever only one of the two cylinders as was the case for the B90/Nuclear Depth Strike Bomb. Therefore, code modeling and simulation technology must be utilized to obtain a better understanding of the LSC jet hydrodynamic penetration, fracture, shear, and spall mechanisms associated with the severance of metallic as well as composite targets.

Vigil, Manuel G.↗

Giant impacts on a primitive Uranus

Simulations of collisions are conducted between a model of the primitive Uranus and 1-3 earth-mass impactors, using smooth-particle hydrodynamics. A series of collisions was simulated for each impactor while varying the total angular momentum of the system. Most of the simulation runs left ices in orbit; a subset of the runs also left rock or iron (from the impactor). It is concluded on the basis of these results that there is a wide range of giant impacts which could have produced the current period and inclination of the spin axis relative to the plane of the ecliptic. A subset of these could have deposited the material in orbit from which the regular satellites of Uranus were assembled.

Slattery, Wayne L.↗

Controlling the Mobility of the Fluid Interface of Moving Gas Bubbles or Liquid Drops by Using Micellar Solutions of Surfactants

Microgravity processes must rely on mechanisms other than buoyancy to move bubbles or droplets from one region to another in a continuous liquid phase. One suggested method is thermocapillary migration in which a temperature gradient is applied to the continuous phase. A significant and as yet unresolved impediment to the use of thermocapillary migration to direct bubble or drop motion is that these migrations can be significantly retarded by the adsorption onto the fluid particle surface of surface active impurities unavoidably present in the continuous or (if the particle is a liquid) droplet phases. The focus of our research was to develop a theory for remobilizing fluid particle interfaces retarded by a surfactant impurity in an effort to make more viable the use of thermocapillary migrations for the management of bubbles and drops in microgravity. We postulated that a surfactant at high bulk concentration which kinetically exchanges rapidly with the surface can restore interface mobility. The scaling arguments along with a discussion of the previous literature is reviewed in the context of the scaling framework. The specific objectives of the research were twofold. The first was to prove the remobilization theory by studying a model problem. As the mechanism for remobilization is independent of the force which drives the particle, the fluid particle shape and the presence of fluid inertia, we chose the simplest model consisting of a spherical bubble rising steadily by buoyancy in creeping flow. We solved the hydrodynamic and surfactant transport equations for rapid kinetic exchange to demonstrate that as the concentration increases, the Marangoni retardation at first increases (the low k behavior) and then decreases (the high k behavior). The second objective was to develop a method to determine the kinetic rate constants of a surfactant molecule, since this information is necessary to select surfactants which will exchange rapidly enough relative to the convective rate in the thermocapillary process of interest. To measure the kinetic rate, we measure the dynamic tension change accompanying adsorption onto an initially clean interface, or the re-equilibration in tension when an equilibrium interface is compressed. The dynamic tension measurements are made by a pendant bubble method, in which surfactant adsorbs onto a pendant bubble, and the tension is measured by analyzing the shape change in the bubble. We conclude this report by detailing the publications, presentations and doctoral thesis completed under the auspices of this grant.

Maldarelli, Charles↗

Tank Tests of 1/5.5-Scale Forward Dynamic Model of the Columbia XJL-1 Amphibian - Langley Tank Model 208, TED No. NACA 2336

Tests of a powered dynamic model of the Columbia XJL-1 amphibian were made in Langley tank no.1 to determine the hydrodynamic stability and spray characteristics of the basic hull and to investigate the effects of modifications on these characteristics. Modifications to the forebody chime flare, the step, and the afterbody, and an increase in the angle of incidence of the wing were included in the test program. The seaworthiness and spray characteristics were studied from simulated taxi runs in smooth and rough water. The trim limits of stability, the range of stable positions of the enter of gravity for take-off, and the landing stability were determined in smooth water. The aerodynamic lift, pitching moment, and thrust were determined at speeds up to take-off speed.

Havens, Robert F.↗

Formation of Sharp Eccentric Rings in Debris Disks with Gas but Without Planets

'Debris disks' around young stars (analogues of the Kuiper Belt in our Solar System) show a variety of non-trivial structures attributed to planetary perturbations and used to constrain the properties of those planets. However, these analyses have largely ignored the fact that some debris disks are found to contain small quantities of gas, a component that all such disks should contain at some level. Several debris disks have been measured with a dust-to-gas ratio of about unity, at which the effect of hydrodynamics on the structure of the disk cannot be ignored. Here we report linear and nonlinear modelling that shows that dust-gas interactions can produce some of the key patterns attributed to planets. We find a robust clumping instability that organizes the dust into narrow, eccentric rings, similar to the Fomalhaut debris disk. The conclusion that such disks might contain planets is not necessarily required to explain these systems.

Fomalhaut debris disk↗

An Extension of the Truncated Ingard-Myers Impedance Boundary Condition for High Mach Number Grazing Flows

The Ingard-Myers impedance boundary condition is widely recognized to be prone to hydrodynamic Kelvin-Helmholtz-type instability, primarily due to its use of a vortex sheet to model the flow at the boundary. Recently, a stabilized approximation of the Ingard-Myers condition, called the Truncated Ingard-Myers Impedance Boundary Condition (TIMIBC), has been introduced. This approximation offers a good representation of the Ingard-Myers condition for grazing flows characterized by low to mid subsonic Mach numbers. This paper explores an extension of the TIMIBC, referred to as TIMIBC-ext, designed for liners in grazing flows with high subsonic Mach numbers. The TIMIBC-ext introduces a tunable parameter and demonstrates that by judiciously selecting the parameter value, the accuracy of the TIMIBC can be improved for either upstream or downstream propagating waves in grazing flows with high subsonic Mach numbers. Consequently, the TIMIBC-ext represents an enhancement to the TIMIBC when the location of the source relative to the liner (upstream or downstream) is known. Time domain implementations of the TIMIBC-ext are also provided, employing a multipole expansion model for the impedance function. A numerical example is provided to assess the performance of the TIMIBC-ext, comparing computational results with measurements from a recent set of measurements in the NASA Langley Grazing Flow Impedance Tube.

time domain acoustic scattering↗