Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Hydrodynamic 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 37 records · Page 2

Multistream hydrodynamic modeling of interhemispheric plasma flow

Interhemispheric plasma flow was simulated using one-stream and two-stream hydrodymic models in order to test the suggestion of Banks et al. (1971) and others that the collision of high-speed flows originating from the conjugate hemispheres will cause the formation of a pair of shocks. The single-fluid hydrodynamic equations were modified to include multiple ion streams, allowing for the possibility of counterstreaming flow. It was found that a counterstreaming of ion streams from conjugate hemispheres does occur during the early stages of the refilling of plamaspheric flux tubes, and that a pair of reverse shocks does form. These shocks form away from the equator, and their subsequent motion creates conditions similar to those predicted by the single-stream hydrodynamic models. The findings support the conclusion of earlier studies that the refilling of the plasmasphere occurs from the equatorial region downward.

Rasmussen, C. E.

An Axisymmetric Hydrodynamical Model for the Torus Wind in AGN: X-ray Excited Funnel Flow - 2

We have calculated a series of models of outflows from the obscuring torus in active galactic nuclei (AGN). Our modeling assumes that the inner face of a rotationally supported torus is illuminated and heated by the intense X-rays from the inner accretion disk and black hole. As a result of such heating a strong biconical outflow is observed in our simulations. We calculate 3-dimensional hydrodynamical models, assuming axial symmetry, and including the effects of X-ray heating, ionization, and radiation pressure. We discuss the behavior of a large family of these models, their velocity fields, mass fluxes and temperature, as functions of the torus properties and X-ray flux. Synthetic warm absorber spectra are calculated, assuming pure absorption, for sample models at various inclination angles and observing times. We show that these models have mass fluxes and flow speeds which are comparable to those which have been inferred from observations of Seyfert 1 warm absorbers, and that they can produce rich absorption line spectra.

Dorodnitsyn, A.

A fluctuating hydrodynamics model for nanoscale surfactant-laden interfaces

A multispecies diffuse interface model is formulated in a fluctuating hydrodynamics framework for the purpose of simulating surfactant interfaces at the nanoscale. The model generalizes previous work to ternary mixtures, employing a Cahn-Hilliard free energy density combined with incompressible, isothermal fluctuating hydrodynamics where dissipative fluxes include both deterministic and stochastic terms. The intermolecular parameters in the free energy are chosen such that one species acts as a partially miscible surfactant. From Laplace pressure measurements, we show that in this model the surface tension decreases linearly with surfactant concentration, leading to Marangoni convection for interfaces with concentration gradients. In the capillary wave spectrum for interfaces with and without surfactant, we find that for the former, the spectrum deviates significantly from classical capillary wave theory, presumably due to Gibbs elasticity. In non-equilibrium simulations of the Rayleigh-Plateau instability, deterministic simulations showed that the surfactant delays pinching of a fluid cylinder into droplets. However, stochastic simulations indicate that thermal fluctuations disrupt the surfactant's stabilizing effect. Similarly, the spreading of a patch of surfactant, driven by Marangoni convection, was found to be partially suppressed by thermal fluctuations.

Capillary waves

CNO abundances and hydrodynamic models of the nova outburst.

We have used a fully implicit, Lagrangian, hydrodynamic computer code incorporating a nuclear reaction network to follow thermonuclear runaways in the hydrogen-rich envelopes of white dwarfs in order to produce a nova outburst. Because of the short time-scales and the high nuclear burning rates produced in our models, the nuclear reactions are far out of equilibrium and the beta-plus unstable nuclei become the most abundant nuclei in the envelope except for hydrogen and helium. Our models have ejected 1.00017 solar mass with kinetic energies of 8 times 10 to the 44-th power ergs, a value that agrees quite closely with the observed values for novae.

Starrfield, S.

Hydrodynamic models of a cepheid atmosphere

A method for including the solution of the transfer equation in a standard Henyey type hydrodynamic code was developed. This modified Henyey method was used in an implicit hydrodynamic code to compute deep envelope models of a classical Cepheid with a period of 12(d) including radiative transfer effects in the optically thin zones. It was found that the velocity gradients in the atmosphere are not responsible for the large microturbulent velocities observed in Cepheids but may be responsible for the occurrence of supersonic microturbulence. It was found that the splitting of the cores of the strong lines is due to shock induced temperature inversions in the line forming region. The adopted light, color, and velocity curves were used to study three methods frequently used to determine the mean radii of Cepheids. It is concluded that an accuracy of 10% is possible only if high quality observations are used.

Karp, A. H.

Sharp Eccentric Rings in Planetless Hydrodynamical Models of Debris Disks

Exoplanets are often associated with disks of dust and debris, analogs of the Kuiper Belt in our solar system. These "debris disks" show a variety of non-trivial structures attributed to planetary perturbations and utilized to constrain the properties of the planets. However, analyses of these systems have largely ignored the fact that, increasingly, debris disks are found to contain small quantities of gas, a component all debris disks should contain at some level. Several debris disks have been measured with a dust-to-gas ratio around unity where the effect of hydrodynamics on the structure of the disk cannot be ignored. Here we report that dust-gas interactions can produce some of the key patterns seen in debris disks that were previously attributed to planets. Through linear and nonlinear modeling of the hydrodynamical problem, we find that a robust clumping instability exists in this configuration, organizing the dust into narrow, eccentric rings, similar to the Fomalhaut debris disk. The hypothesis that these disks might contain planets, though thrilling, is not necessarily required to explain these systems.

debris disks

3D Radiative Hydrodynamics Modeling of Convection of Stars to Probe Their Interiors and Photospheric Properties

The dramatic flow of data from the Kepler and K2 missions opens the opportunity to significantly improve our knowledge of stellar interiors, surface dynamics, and structure. However, interpretation of these observations is a challenging task because it depends on tiny effects that can be studied only with advanced first-principles modeling. We present results of 3D time-dependent radiative hydrodynamic simulations of stellar outer convection zones and atmospheres taking into account chemical composition, radiative transfer, turbulence effects, and a realistic equation of state for main sequence stars. We will discuss properties of convective structure and dynamics, convective overshoot, effects of magnetic fields and rotation, as well as the potential influence of turbulent surface dynamics on high-precision RV measurements.

Photospheric

Hydrodynamic Modelling for Tidal Energy Kites: Cooperative Research and Development Final Report, CRADA Number CRD-24-30123

The primary goal of the work is to improve accuracy of hydrodynamic loads acting on the lifting surface of the kite turbine due to complex air-water-lifting body interaction in a low-order numerical model using high-fidelity numerical input. The result of this study will help the Participant better design control algorithms for the kite turbine so that the turbine can harvest more energy while maintaining stability of the platform as well as its reliability. The use of hydrodynamic coefficients/load components from computational fluid dynamic simulation will inform better inputs for a medium-fidelity dynamic model (specifically an OrcaFlex model) which is intended for use in rapid design iterations in the design process.

16 TIDAL AND WAVE POWER

A hydrodynamic model of an outer hair cell

On the model it is possible to measure the force and the force direction for each individual hair as a function of the flow direction and velocity. Measurements were made at the man flow velocity .01 m/s, which is equivalent to a flow velocity in the real ear of about 1 micrometer/s. The kinematic viscosity of the liquid used in the model was 10,000 times higher than the viscosity of perilymph to attain hydrodynamic equality. Two different geometries for the sterocilia pattern were tested. First the force distribution for a W-shaped sterocilia pattern was recorded. This is the sterocilia pattern found in all real ears. It is found that the forces acting on the hairs are very regular and perpendicular to the legs of the W when the flow is directed from the outside of the W. When the flow is reversed, the forces are not reversed, but are much more irregular. This can eventually explain the half wave rectification of the nerve signals. As a second experiment, the force distribution for a V-shaped sterocilia pattern was recorded. Here the forces were irregular both when the flow was directed into the V and when it was directed against the edge of the V.

Jacobson, B. O.

Hydrodynamic Modeling of Oxidizer-Rich Staged Combustion Injector Flow

The main objective of this work is to determine the unsteady hydrodynamic characteristics of coaxial swirl atomizers of interest in oxidizer-rich staged combustion (ORSC) liquid rocket engines. To this end, the pseudo-density (homogeneous flow) treatment combined with the Marker-and-Cell (MAC) numerical algorithm has been used to develop an axisymmetric with swirl, two-phase, unsteady model. The numerical model is capable of assessing the time-dependent orifice exit conditions and internal mixing for arbitrary fuel and oxidizer gas injection conditions. Parametric studies have been conducted to determine the effect of geometry, gas properties, and liquid properties on the exit massflow rate and velocity. It has been found that the frequency at which the liquid film oscillates increases as the density ratio and thickness increase, decreases as film thickness and liquid swirl velocity increase, and is unaffected by the mixing length. Additionally, it has been determined that the variation in the massflow rate increases as the liquid swirl velocity and liquid film thickness increase, and decreases as the density ratio, collar thickness, and mixing length increase.

Canino, J. V.

Constraining hydrodynamic models of inertial confinement fusion implosions using capsule surrogate experiments

We conduct capsule surrogate experiments at the National Ignition Facility to calibrate radiation hydrodynamic simulations to infer hydrodynamic conditions that are not observable in indirect drive ignition implosions. We tune the simulations by applying laser power and cross beam energy transfer (CBET) saturation multipliers to match the observables from capsule surrogate experiments. Shock timing, velocity, and symmetry are measured in liquid D 2 filled Keyhole capsule surrogate experiments and implosion trajectory, stagnation time, and shape time history are measured in in-flight 2D backlit x-ray radiography experiments (“2DConA”) of D 2 gas filled capsule implosions. Calibrated simulations suggest that the N210808 ignition implosion (fusion target gain = 0.7) had a shell mass remaining at stagnation of less than the nominal %5 (3.8%) and resulted in less confinement. For N221204, the shell was made 5.75 μm thicker to trade implosion velocity for increased confinement and resulted in a target gain = 1.5 with a shell mass remaining of 5.7%. Furthermore, a single adjusted model can reproduce all shock timing data as changes are made to shell thickness (79–85 μm) and laser wavelength separation (1.8–4.0 Å). However, for the 2DConA implosions, a 5% variation in the peak power laser multipliers and a 30% variation in late-time CBET between experiments are needed to match the observed stagnation times, in-flight $P_2$ shape, and hot-spot $P_2$ shape. While progress is being made to improve the models in simulations using focused experiments, capsule surrogate experiments will continue to be needed to optimize future ignition designs.

Lasers

Data-Driven Radiative Hydrodynamic Modeling of the 2014 March 29 X1.0 Solar Flare

Spectroscopic observations of solar flares provide critical diagnostics of the physical conditions in the flaring atmosphere. Some key features in observed spectra have not yet been accounted for in existing flare models. Here we report a data-driven simulation of the well-observed X1.0 flare on 2014 March 29 that can reconcile some well-known spectral discrepancies. We analyzed spectra of the flaring region from the Interface Region Imaging Spectrograph (IRIS) in Mg II hk, the Interferometric BIdimensional Spectropolarimeter at the Dunn Solar Telescope (DSTIBIS) in H(alpha) 6563A and Ca II 8542A, and the Reuven Ramaty High Energy Solar Spectroscope Imager (RHESSI) in hard X-rays. We constructed a multithreaded flare loop model and used the electron flux inferred from RHESSI data as the input to the radiative hydrodynamic code RADYN to simulate the atmospheric response. We then synthesized various chromospheric emission lines and compared them with the IRIS and IBIS observations. In general, the synthetic intensities agree with the observed ones, especially near the northern footpoint of the flare. The simulated Mg II line profile has narrower wings than the observed one. This discrepancy can be reduced by using a higher microturbulent velocity (27 km/s) in a narrow chromospheric layer. In addition, we found that an increase of electron density in the upper chromosphere within a narrow height range of approx. 800 km below the transition region can turn the simulated Mg II line core into emission and thus reproduce the single peaked profile, which is a common feature in all IRIS flares.

Da Costa, Fatima Rubio

Hydrodynamic modeling of plasma channel systems for laser plasma accelerators

Structured plasma channels are an essential technology for driving high-gradient, plasma-based acceleration and control of electron and positron beams for advanced concepts accelerators. Laser and gas technologies can permit the generation of long plasma columns known as hydrodynamic, optically-field-ionized (HOFI) channels, which feature low on-axis densities and steep walls. By carefully selecting the background gas and laser properties, one can generate narrow, tunable plasma channels for guiding high intensity laser pulses. Here, we present on the development of simulations of HOFI channels using the FLASH code, a publicly available radiation hydrodynamics code. We explore sensitivities of the channel evolution to laser profile, intensity, and background gas conditions, and identify relevant scalings with laser intensity through a range of practical channel delays.

Accelerators

AGN Obscuration Through Dusty Infrared Dominated Flows. II. Multidimensional, Radiation-Hydrodynamics Modeling

We explore a detailed model in which the active galactic nucleus (AGN) obscuration results from the extinction of AGN radiation in a global ow driven by the pressure of infrared radiation on dust grains. We assume that external illumination by UV and soft X-rays of the dusty gas located at approximately 1pc away from the supermassive black hole is followed by a conversion of such radiation into IR. Using 2.5D, time-dependent radiation hydrodynamics simulations in a ux-limited di usion approximation we nd that the external illumination can support a geometrically thick obscuration via out ows driven by infrared radiation pressure in AGN with luminosities greater than 0:05 L(sub edd) and Compton optical depth, Tau(sub T) approx > & 1.

Dorodnitsyn, Anton

Late-Time Evolution of Composite Supernova Remnants: Deep Chandra Observations and Hydrodynamical Modeling of a Crushed Pulsar Wind Nebula in SNR G327.1-1.1

In an effort to better understand the evolution of composite supernova remnants (SNRs) and the eventual fate of relativistic particles injected by their pulsars, we present a multifaceted investigation of the interaction between a pulsar wind nebula (PWN) and its host SNR G327.1-1.1. Our 350 ks Chandra X-ray observations of SNR G327.1-1.1 reveal a highly complex morphology; a cometary structure resembling a bow shock, prong-like features extending into large arcs in the SNR interior, and thermal emission from the SNR shell. Spectral analysis of the non-thermal emission offers clues about the origin of the PWN structures, while enhanced abundances in the PWN region provide evidence for mixing of supernova ejecta with PWN material. The overall morphology and spectral properties of the SNR suggest that the PWN has undergone an asymmetric interaction with the SNR reverse shock(RS) that can occur as a result of a density gradient in the ambient medium and or a moving pulsar that displaces the PWN from the center of the remnant. We present hydrodynamical simulations of G327.1-1.1 that show that its morphology and evolution can be described by a approx. 17,000 yr old composite SNR that expanded into a density gradient with an orientation perpendicular to the pulsar's motion. We also show that the RSPWN interaction scenario can reproduce the broadband spectrum of the PWN from radio to gamma-ray wavelengths. The analysis and modeling presented in this work have important implications for our general understanding of the structure and evolution of composite SNRs.

Evolution

A Transient Hydrodynamic Model of Screen Channel Liquid Acquisition Devices for In-Space Cryogenic Propellant Management

Screen channel liquid acquisition devices (LADs) will play a crucial role in future deep space travel. It is essential that vapor-free delivery of propellants during tank-to-tank transfer is ensured to maximize yield from storage tanks and prevent potential combustion instabilities. The screen channel LAD utilizes a fine screen wire mesh that can separate phases in a low Bond number (i.e. microgravity) environment using surface tension forces. This study presents the development and verification of a new model for transient screen compliance, one of the influential factors for screen channel LAD design. Screen compliance is crucial during LAD channel outflow transients because the slight deflection of the screen can provide needed mass to satisfy rapid outflow demands and reduce the pressure difference across the screen. The model is successfully verified against CFD simulations. In addition, the characteristic speed for the governing screen compliance equations is derived which allows for numerical stability criteria to be established. As shown in this study, the transient maximum pressure difference across the screen can greatly exceed the steady state maximum pressure difference across the screen in many cases.

Hydrodynamics Simulations

Hydrodynamic Modeling of Solar Acoustic Wave Propagation

Acoustic waves in the sun are separated into two categories: global and local. Global oscillations are excited through convection, and resonate in the solar interior. As with all resonant media, propagation of these acoustic waves is limited to specific frequencies for a given wavenumber. These modes are easily distinguished in 1-v diagrams, such as in figure 1. Global modes dominate the spectrum in the 2-5 mHz range. Local oscillations, so-called "sunquakes", are usually excited by particularly strong flares and can be observed as ring-shaped waves propagating away from the point of impact. The propagation is best displayed in a time-distance diagram, where a slit passing from the source to some distance away. Distance is displayed on the x-axis, and the time-series is plotted along the y-axis; an example is shown in figure 1, showing the sunquake induced by the Sept 2017 X-class flare.

Hydrodynamic

exoALMA. XI. ALMA Observations and Hydrodynamic Models of LkCa 15: Implications for Planetary Mass Companions in the Dust Continuum Cavity

In the past decade, the Atacama Large Millimeter/submillimeter Array (ALMA) has revealed a plethora of substructures in the disks surrounding young stars. These substructures have several proposed formation mechanisms, with one leading theory being the interaction between the disk and newly formed planets. In this Letter, we present high angular resolution ALMA observations of LkCa 15’s disk that reveal a striking difference in dust and CO emission morphology. The dust continuum emission shows a ringlike structure characterized by a dust-depleted inner region of ~40 au in radius. Conversely, the CO emission is radially smoother and shows no sign of gas depletion within the dust cavity. We compare the observations with models for the disk–planet interaction, including radiative transfer calculation in the dust and CO emission. This source is particularly interesting, as the presence of massive planets within the dust cavity has been suggested based on previous near-IR observations. We find that the level of CO emission observed within the dust cavity is inconsistent with the presence of planets more massive than Jupiter orbiting between 10 and 40 au. Instead, we argue that the LkCa 15 innermost dust cavity might be created either by a chain of low-mass planets or by other processes that do not require the presence of planets.

79 ASTRONOMY AND ASTROPHYSICS