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 307 records · Page 17

The role of hydrogen thermalization in shaping the Lyman-alpha coma

Two interrelated aspects of the thermalization of cometary hydrogen were addressed using a three-dimensional Monte Carlo particle trajectory model: photochemical heating of the coma and the velocity distribution of cometary hydrogen leaving the inner coma that produces the extended Lyman-alpha coma. Observed cometary conditions were examined, ranging from productive comets like Kohoutek at very small heliocentric distances to intrinsically small comets like Giacobini-Zinner at 1 AU. The collisional decoupling of hydrogen yields bulk flow speeds and gas temperatures which are smaller than generally predicted in hydrodynamic models. By properly including these conditions the hydrogen velocity distribution required to produce the observed Lyman-alpha coma can be derived naturally. The two-dimensional morphology of the H coma is diagnostic of the physical conditions in the inner coma.

Combi, Michael R.↗

Analysis of Broad-Lined Type Ic Supernovae From the (Intermediate) Palomar Transient Factory

We study 34 Type Ic supernovae that have broad spectral features (SNe Ic-BL). This is the only SN type found in association with long duration gamma-ray bursts (GRBs).We obtained our photometric data with the Palomar Transient Factory (PTF) and its continuation, the intermediate PTF (iPTF). This is the first large, homogeneous sample of SNe Ic-BL from an untargeted survey. Furthermore, given the high observational cadence of iPTF, most of these SNe Ic-BL were discovered soon after explosion. We present K-corrected Bgriz light curves of these SNe, obtained through photometry on template-subtracted images. We analyzed the shape of the r-band light curves, finding a correlation between the decline parameter delta m15 and the rise parameter delta m-10. We studied the SN colors and, based on g - r, we estimated the host-galaxy extinction for each event. Peak r-band absolute magnitudes have an average of -18:6 +- 0:5 mag. We fit each r-band light curve with that of SN 1998bw (scaled and stretched) to derive the explosion epochs. We computed the bolometric light curves using bolometric corrections, r-band data, and g - r colors. Expansion velocities from Fe ii were obtained by fitting spectral templates of SNe Ic. Bolometric light curves and velocities at peak were fitted using the semianalytic Arnett model to estimate ejecta mass Mej, explosion energy EK and 56Ni mass M(56Ni) for each SN. We find average values of Mej = 4 +- 3 M sun, EK = (7 +- 6) x 1051 erg, and M(56Ni) = 0:31 +- 0:16 M . The parameter distributions were compared to those presented in the literature and are overall in agreement with them. We also estimated the degree of 56Ni mixing using scaling relations derived from hydrodynamical models and we find that all the SNe are strongly mixed. The derived explosion parameters imply that at least 21% of the progenitors of SNe Ic-BL are compatible with massive (>28 M ), possibly single stars, whereas at least 64% might come from less massive stars in close binary systems.

Supernovae↗

Data Assimilation and Transport Modeling in Terrestrial and Planetary Atmospheres

Data assimilation is a blanket term used to describe a number of techniques for retrieving important physical parameters from observational data, subject to constraints imposed by prior knowledge (such as, in the case of meteorology, the primitive equations that govern atmospheric motion). Since these newly developed methods make efficient use of computational resources, they are of great importance in the interpretation of the voluminous datasets that are now produced by satellite missions. As proposed, these techniques have been applied to the study of the Martian and terrestrial atmospheres based on available satellite observations. In addition, a sophisticated hydrodynamic model (non-hydrostatic, and therefore applicable to the study of the interiors of the giant planets) has also been developed and successfully applied to the study of tidally induced motions in Jupiter.

Houben, Howard C.↗

Spectral Eclipse Timing

We utilize multi-dimensional simulations of varying equatorial jet strength to predict wavelength-dependent variations in the eclipse times of gas-giant planets. A displaced hot spot introduces an asymmetry in the secondary eclipse light curve that manifests itself as a measured offset in the timing of the center of eclipse. A multiwavelength observation of secondary eclipse, one probing the timing of barycentric eclipse at short wavelengths and another probing at longer wavelengths, will reveal the longitudinal displacement of the hot spot and break the degeneracy between this effect and that associated with the asymmetry due to an eccentric orbit. The effect of time offsets was first explored in the IRAC wavebands by Williams et al. Here we improve upon their methodology, extend to a broad range of wavelengths, and demonstrate our technique on a series of multi-dimensional radiative-hydrodynamical simulations of HD 209458b with varying equatorial jet strength and hot-spot displacement. Simulations with the largest hot-spot displacement result in timing offsets of up to 100 s in the infrared. Though we utilize a particular radiative hydrodynamical model to demonstrate this effect, the technique is model independent. This technique should allow a much larger survey of hot-spot displacements with the James Webb Space Telescope than currently accessible with time-intensive phase curves, hopefully shedding light on the physical mechanisms associated with thermal energy advection in irradiated gas giants.

Hydrodynamics↗

Tide Corrections for Coastal Altimetry: Status and Prospects

Knowledge of global oceanic tides has markedly advanced over the last two decades, in no small part because of the near-global measurements provided by satellite altimeters, and especially the long and precise Topex/Poseidon time series e.g. [2]. Satellite altimetry in turn places very severe demands on the accuracy of tidal models. The reason is clear: tides are by far the largest contributor to the variance of sea-surface elevation, so any study of non-tidal ocean signals requires removal of this dominant tidal component. Efforts toward improving models for altimetric tide corrections have understandably focused on deep-water, open-ocean regions. These efforts have produced models thought to be generally accurate to about 2 cm rms. Corresponding tide predictions in shelf and near-coastal regions, however, are far less accurate. This paper discusses the status of our current abilities to provide near-global tidal predictions in shelf and near-coastal waters, highlights some of the difficulties that must be overcome, and attempts to divine a path toward some degree of progress. There are, of course, many groups worldwide who model tides over fairly localized shallow-water regions, and such work is extremely valuable for any altimeter study limited to those regions, but this paper considers the more global models necessary for the general user. There have indeed been efforts to patch local and global models together, but such work is difficult to maintain over many updates and can often encounter problems of proprietary or political nature. Such a path, however, might yet prove the most fruitful, and there are now new plans afoot to try again. As is well known, tides in shallow waters tend to be large, possibly nonlinear, and high wavenumber. The short spatial scales mean that current mapping capabilities with (multiple) nadir-oriented altimeters often yield inadequate coverage. This necessitates added reliance on numerical hydrodynamic models and data assimilation, which in turn necessitates very accurate bathymetry with high spatial resolution. Nonlinearity means that many additional compound tides and overtides must be accounted for in our predictions, which increases the degree of modeling effort and increases the amounts of data required to disentangle closely aliased tides.

Ray, Richard D.↗

A Survey of Multimaterial Treatments for Thermal Radiative Transfer

Arbitrary Lagrangian-Eulerian methods are a popular choice for hydrodynamic modeling in radiation (rad-hydro) simulations. Because these methods involve a relaxation step that moves the mesh relative to material boundaries, multimaterial spatial zones are generally present. Accurate treatments of these zones are needed to resolve various physical phenomena of interest for inertial confinement fusion applications. However, these codes are often paired with single-material, deterministic thermal radiative transfer (TRT) codes that are oblivious to the material compositions of each zone. These single-material TRT codes can only accept homogenized material properties (opacities, specific heats, etc.) from the hydrodynamic code and output homogenized solutions. After each TRT time step, the multimaterial hydrodynamic code must dehomogenize the quantities computed by the TRT package in order to update subzonal material temperatures. The process by which hydrodynamic codes perform this dehomogenization has not been well documented in previous literature, and the methods can vary significantly from code to code. The purpose of this paper is to document, study, and compare existing techniques used for rad-hydro simulations as well as present a new method with potentially promising results. We summarize several methods and give comparisons on infinite-medium problems as well a finite-medium problem for two of the methods.

42 ENGINEERING↗

Understanding spatial and temporal drivers of variation in tree hydraulic processes and their consequences for climate feedbacks (Final Technical Report)

This is the final technical report from the first phase of a project that changed institutions. The grant was titled “Understanding spatial and temporal drivers of variation in tree hydraulic processes and their consequences for climate feedbacks.” The overall objectives of this project were to (1) provide model‐compatible datasets of key plant hydraulic traits and status for model evaluation, parameterization and validation and (2) use these data to pinpoint ecosystem responses to a changing hydroclimate by addressing both long‐term climatic drying and episodic extreme droughts. We planned to address the objectives with three research activities to quantify plant responses to chronic water stress and episodic drought: (1) generate high frequency observations of soil and plant hydraulic data across different landscape positions at multiple sites, (2) quantify plant hydraulic trait plasticity in response to experimental soil moisture reduction in situ in two central hardwood forests, and (3) simulate the carbon consequences of incorporating plant hydrodynamics and plant acclimation to water stress in the DOE‐sponsored plant hydrodynamics model FATES‐HYDRO. As of the transfer of this project to another institution, we had made substantial progress on activities 1 and 2, and started activity 3.

54 ENVIRONMENTAL SCIENCES↗

Integrated Data Analysis to Improve Te and Zeff Measurement on NSTX-U

The multi-energy soft x-ray (ME-SXR) camera is a new concept for a fusion plasma diagnostic based on a pixelated detector with adjustable photon cutoff thresholds. It is a single, compact, solid-state detector which combines spatial and spectral sensitivity and can be configured remotely at run time, making it ideal for a fusion environment. The detector can be configured for application to a wide range of physics problems including impurity transport, electron temperature profiles, and non-thermal effects. The detector has been installed on the Madison Symmetric Torus (MST) as part of a collaboration between Princeton Plasma Physics Laboratory and the University of Wisconsin–Madison to develop the technology and analysis techniques. A full diagnostic forward model of the ME-SXR diagnostic has been developed to enable analysis of the raw photon count data. This model incorporates the atomic physics underlying the plasma emission, diagnostic geometry, and detector response characteristics. Employing this model, a Bayesian analysis approach is used to extract electron temperature from the raw data. This analysis approach uses all available information, clearly quantifies prior information, and handles uncertainties transparently. The ME-SXR camera is being used to study the 3D (helical, stellarator-like) state that forms in the MST device. This helical core structure has improved energy confinement. The transition to the helical state is accompanied by a rise in electron temperature in the central part of the helical structure, but then this enhancement fluctuates. The characteristics of this fluctuation are being compared to a magneto-hydrodynamics model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Time-dependent MHD modeling of the global structure of the heliosphere

We present results from time-dependent modeling of the global structure of the heliosphere with neutral and magnetic field effects included. The magnetic field is assumed parallel to the interstellar flow in this two-dimensional axisymmetric model; the neutrals are treated as a fluid. The effects of interstellar neutrals and the interplanetary magnetic field on the location of the termination shock are studied using the most recent estimate of the interstellar medium parameters, results will be compared to those of Baranov and Zaitsev. The effect of the solar wind - VLISM interaction on the density and velocity of interstellar neutrals within the heliosphere will also be presented and related to observations. The response of the termination shock to the solar cycle variation in the solar wind will be compared to the response found previously using an axisymmetric hydrodynamic model without neutrals.

Liewer, P. C.↗

Pros and Cons of ID vs. 3D Modeling

Advances in computing capability have led to tremendous improvements in 3D modeling. Entire active regions are being simulated in what might be described as a first principles way, in which plasma heating is treated self consistently rather than through the specification of heating functions. There are limitations to this approach, however, as actual heating mechanisms on the Sun involve spatial scales orders of magnitude smaller than what these simulations can resolve. Other simulations begin to resolve these scales, but they only treat a tiny volume and do not include the all important coupling with larger scales or with other parts of the atmosphere, and so cannot be readily compared with observations. Finally, ID hydrodynamic models capture the field-aligned evolution of the plasma extremely well and are ideally suited for data comparison, but they treat the heating in a totally ad hoc manner. All of these approaches have important contributions to make, but we must be aware of their limitations. I will highlight some of the strengths. and weaknesses of each.

Klimchuk, James A.↗

Elliptic flow of charged particles at midrapidity relative to the spectator plane in Pb–Pb and Xe–Xe collisions

Measurements of the elliptic flow coefficient relative to the collision plane defined by the spectator neutrons v 2 {$Ψ_{SP}$} in collisions of Pb ions at center-of-mass energy per nucleon–nucleon pair $\sqrt{s_{NN}}=$ 2.76 TeV and Xe ions at $\sqrt{s_{NN}}=$ 5.44 TeV are reported. The results are presented for charged particles produced at midrapidity as a function of centrality and transverse momentum for the 5–70% and 0.2–6 GeV/c ranges, respectively. The ratio between v 2 {$Ψ_{SP}$} and the elliptic flow coefficient relative to the participant plane v 2 {4}, estimated using four-particle correlations, deviates by up to 20% from unity depending on centrality. This observation differs strongly from the magnitude of the corresponding eccentricity ratios predicted by the T R ENTo and the elliptic power models of initial state fluctuations that are tuned to describe the participant plane anisotropies. The differences can be interpreted as a decorrelation of the neutron spectator plane and the reaction plane because of fragmentation of the remnants from the colliding nuclei, which points to an incompleteness of current models describing the initial state fluctuations. A significant transverse momentum dependence of the ratio v 2 {$Ψ_{SP}$} /v 2 {4} is observed in all but the most central collisions, which may help to understand whether momentum anisotropies at low and intermediate transverse momentum have a common origin in initial state fluctuations. The ratios of v 2 {$Ψ_{SP}$} and v 2 {4} to the corresponding initial state eccentricities for Xe–Xe and Pb–Pb collisions at similar initial entropy density show a difference of (7.0 ± 0.9)% with an additional variation of +1.8% when including RHIC data in the T R ENTo parameter extraction. These observations provide new experimental constraints for viscous effects in the hydrodynamic modeling of the expanding quark–gluon plasma produced in heavy-ion collisions at the LHC.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Using Remotely Sensed Data and Hydrologic Models to Evaluate the Effects of Climate Change on Shallow Aquatic Ecosystems in the Mobile Bay, AL Estuary

Coastal systems in the northern Gulf of Mexico, including the Mobile Bay, AL estuary, are subject to increasing pressure from a variety of activities including climate change. Climate changes have a direct effect on the discharge of rivers that drain into Mobile Bay and adjacent coastal water bodies. The outflows change water quality (temperature, salinity, and sediment concentrations) in the shallow aquatic areas and affect ecosystem functioning. Mobile Bay is a vital ecosystem that provides habitat for many species of fauna and flora. Historically, submerged aquatic vegetation (SAV) and seagrasses were found in this area of the northern Gulf of Mexico; however the extent of vegetation has significantly decreased over the last 60 years. The objectives of this research are to determine: how climate changes affect runoff and water quality in the estuary and how these changes will affect habitat suitability for SAV and seagrasses. Our approach is to use watershed and hydrodynamic modeling to evaluate the impact of climate change on shallow water aquatic ecosystems in Mobile Bay and adjacent coastal areas. Remotely sensed Landsat data were used for current land cover land use (LCLU) model input and the data provided by Intergovernmental Panel on Climate Change (IPCC) of the future changes in temperature, precipitation, and sea level rise were used to create the climate scenarios for the 2025 and 2050 model simulations. Project results are being shared with Gulf coast stakeholders through the Gulf of Mexico Data Atlas to benefit coastal policy and climate change adaptation strategies.

Estes, M. G.↗

Design and modeling of indirectly driven magnetized implosions on the NIF

The use of magnetic fields to improve the performance of hohlraum-driven implosions on the National Ignition Facility (NIF) is discussed. The focus is on magnetically insulated inertial confinement fusion, where the primary field effect is to reduce electron-thermal and alpha-particle loss from the compressed hotspot (magnetic pressure is of secondary importance). We summarize the requirements to achieve this state. The design of recent NIF magnetized hohlraum experiments is presented. These are close to earlier shots in the three-shock, high-adiabat (BigFoot) campaign, subject to the constraints that magnetized NIF targets must be fielded at room-temperature, and use ≲1 MJ of laser energy to avoid the risk of optics damage from stimulated Brillouin scattering. We present results from the original magnetized hohlraum platform, as well as a later variant that gives a higher hotspot temperature. In both platforms, imposed fields (at the capsule center) of up to 28 T increase the fusion yield and hotspot temperature. Integrated radiation-magneto-hydrodynamic modeling with the Lasnex code of these shots is shown, where laser power multipliers and a saturation clamp on cross-beam energy transfer are developed to match the time of peak capsule emission and the P2 Legendre moment of the hotspot x-ray image. The resulting fusion yield and ion temperature agree decently with the measured relative effects of the field, although the absolute simulated yields are higher than the data by 2.0−2.7×. The tuned parameters and yield discrepancy are comparable for experiments with and without an imposed field, indicating the model adequately captures the field effects. Self-generated and imposed fields are added sequentially to simulations of one BigFoot NIF shot to understand how they alter target dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Stochastic hydrodynamics and hydro-kinetics: Similarities and differences

The hydro-kinetic formalism has been used as a complementary approach to solving the Stochastic Differential Equations (SDE) corresponding to noisy hydrodynamics. The hydro-kinetic formalism consists of a deterministic set of relaxation type equations that tracks the evolution of 2-point correlation functions of stochastic hydrodynamic quantities. Hence they are comparatively easier to solve than the SDEs, which are computationally intensive and need to deal with arbitrarily large gradients. This work compares the two approaches for the propagation and diffusion of conserved charge fluctuations in the Bjorken hydrodynamic model. For white noise, the two approaches agree. For colored Catteneo noise, which is causal, the two approaches diverge. This is because white noise only induces two-point correlations, while Catteneo noise also induces higher-order correlations. Furthermore, this difference is quantified from the effects of causal evolution and influence from higher-order correlations induced by the Catteneo noise.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A Chandra X-Ray Survey of Ejecta in the Cassiopeia A Supernova Remnant

We present a survey of the X-ray emitting ejecta in the Cassiopeia A supernova remnant based on an extensive analysis of over 6000 spectral regions extracted on 2.5-10" angular scales using the Chandra 1 Ms observation. We interpret these results in the context of hydrodynamical models for the evolution of the remnant. The distributions of fitted temperature and ionization age are highly peaked and suggest that the ejecta were subjected to multiple secondary shocks. Based on the fitted emission measure and element abundances, and an estimate of the emitting volume, we derive masses for the X-ray emitting ejecta as well as showing the distribution of the mass of various elements over the remnant. The total shocked Fe mass appears to be roughly 0.14 Solar Mass, which accounts for nearly all of the mass expected in Fe ejecta. We find two populations of Fe ejecta, that associated with normal Si-burning and that associated with alpha-rich freeze-out, with a mass ratio of approximately 2:1. Surprisingly, essentially all of this Fe (both components) is well outside the central regions of the SNR, presumably having been ejected by hydrodynamic instabilities during the explosion. We discuss this, and its implications for the neutron star kick.

Hwang, Una↗

Correlations of conserved quantities at finite baryon density

Correlations involving the seven conserved quantities, namely energy, baryon number, electric charge, strangeness, and the three components of momentum, give rise to correlations in heavy-ion collisions. Through the utilization of a simple one-dimensional hydrodynamic model, we calculate the evolution of the entire 7×7 matrix of correlations as a function of relative spatial rapidity. This comprehensive analysis accounts for finite baryon density, which results in off-diagonal correlations between the charge-related quantities and the energy-momentum quantities. These correlations in coordinate space are subsequently transformed into correlations in momentum space using statistical weighting. As a result, the entire matrix of correlations is revealed to be highly sensitive to the equation of state, viscosity, and diffusivity.

Equations of state of nuclear matter↗

Density measurements for the National Ignition Facility (NIF) opacity platform

The Opacity Platform on the National Ignition Facility (NIF) has been developed to measure opacities at varying densities and temperatures relevant to the solar interior and thermal cooling rates in white dwarf stars. The typical temperatures reached at NIF range between 150 and 210 eV, which allow these measurements to be performed experimentally. The captured opacities are crucial to validating radiation-hydrodynamic models that are used in astrophysics. The NIF opacity platform has a unique new capability that allows in situ measurement of the sample expansion. The sample expansion data are used to better understand the plasma conditions in our experiments by inferring the sample density throughout the duration of the laser drive. Here, we present the details of the density measurement technique, data analysis, and recent results for Fe and MgO.

47 OTHER INSTRUMENTATION↗

Effects of dissipative baryon current in heavy-ion collisions at RHIC-BES energies

The CLVisc (3+1)D viscous hydrodynamic model is extended to include the equation of net baryon conservation and the Israel-Stewart-like equations for dissipative baryon current. Using the NEOSB equation of state, we simulate the dynamical evolution and collectivity of the quark-gluon plasma with finite chemical potential, assuming smooth energy density and net baryon density distributions at the initial proper time. Numerical results are shown for the impact of net-baryon dissipation on particle yields and p T spectra in heavy-ion collisions at beam energy scan energies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗