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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 253 records · Page 14

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Nonlinear Model Predictive Control Based on Real-Time Iteration Scheme for Wave Energy Converters Using WEC-Sim: Preprint

One of several challenges that wave energy technologies face is their inability to generate electricity cost-competitively with other grid-scale energy generation sources. Several studies have identified two approaches to lower the levelised cost of electricity: reduce the cost over the device's lifetime or increase its overall electrical energy production. Several advanced control strategies have been developed to address the latter. However, only a few take into account the overall efficiency of the power take-off (PTO) system, and none of them solve the optimisation problem that arises at each sampling time on real-time. In this paper, a detailed Nonlinear model predictive control (NMPC) approach based on the real-time iteration (RTI) scheme is presented, and the controller performance is evaluated using a time-domain hydrodynamics model (WEC-Sim). The proposed control law incorporates the PTO system's efficiency in a control law to maximise the energy extracted. The study also revealed that RTI-NMPC clearly outperforms a simple resistive controller.

model predictive control↗

DOE OSTI2022

Pulsed power to control growth of silicon nanoparticles in low temperature flowing plasmas

Low-temperature plasmas have seen increasing use for synthesizing high-quality, mono-disperse nanoparticles (NPs). Recent work has highlighted that an important process in NP growth in plasmas is particle trapping—small, negatively charged nanoparticles become trapped by the positive electrostatic potential in the plasma, even if only momentarily charged. In this article, results are discussed from a computational investigation into how pulsing the power applied to an inductively coupled plasma (ICP) reactor may be used for controlling the size of NPs synthesized in the plasma. The model system is an ICP at 1 Torr to grow silicon NPs from an Ar/SiH 4 gas mixture. This system was simulated using a two-dimensional plasma hydrodynamics model coupled to a three-dimensional kinetic NP growth and trajectory tracking model. The effects of pulse frequency and pulse duty cycle are discussed. We identified separate regimes of pulsing where particles become trapped for one pulsed cycle, a few cycles, and many cycles—each having noticeable effects on particle size distributions. For the same average power, pulsing can produce a stronger trapping potential for particles when compared to continuous wave power, potentially increasing particle mono-dispersity. Pulsing may also offer a larger degree of control over particle size for the same average power. Experimental confirmation of predicted trends is discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Assessing Effects of Climate Change on Legacy Waste at the Enewetak Atoll

The Republic of the Marshall Islands (RMI) is in the central Pacific Ocean ~4,500 km west of Hawaii. The Enewetak Atoll, located in the northwest part of the RMI, was the site for 43 nuclear weapon tests between 1948 and 1958. Fallout and deposition from the tests contaminated the island surfaces, lagoon waters and sediment, and nearby ocean waters at the atoll. In the 1970s, a cleanup effort collected radioactive waste and placed it in the Cactus Crater on Runit Island (also called the Runit Dome). In December 2021, Congress directed the U.S. Department of Energy to study the impacts of climate change on the Runit Dome nuclear waste disposal site. Pacific Northwest National Laboratory (PNNL) assembled a multidisciplinary team of climate scientists, ocean modelers, environmental scientists, and health physicists to assess the likely effects of remaining radionuclides at the Enewetak Atoll. PNNL’s approach focused on effects of tropical cyclones that were postulated to mobilize and transport contaminated lagoon sediments and result in human and biota exposure. PNNL’s study estimated (1) the radionuclide source term, (2) the effects of climate change on severe storms, (3) mobilization and transport of radionuclides, and (4) radiation dose to humans and biota. Radionuclides in the lagoon and/or ocean waters of the Enewetak Atoll were characterized by the U.S. Atomic Energy Commission (AEC) in 1972, Woods Hole Oceanographic Institution in 2015, and Lawrence Livermore National Laboratory in 2018. The RMI Nationwide Radiological Study was conducted in the early 1990s for radionuclides remaining in island soils. The 1972 AEC survey remains the most comprehensive source of radionuclide data on lagoon sediments. Climate change modeling at a regional scale in the central Pacific Ocean is limited. PNNL climate scientists simulated severe historical storms postulated to occur both in a recent climate (2015) and in the future (2090) using the Advanced Research Weather Research and Forecasting (WRF-ARW) model, employing a pseudo-global-warming technique. A postulated complete, future failure of the Runit Dome was also considered. PNNL developed a high-resolution regional ocean hydrodynamics model covering the entire RMI extended economic zone using the Finite Volume Coastal Ocean Model (FVCOM). The FVCOM model was run using global reanalysis data for current climate and WRF-ARW simulation for the future climate. PNNL also developed a radionuclide fate and transport model using the FVCOM Integrated Compartment Model (FVCOM-ICM) to simulate the current and future mobilization and transport of radionuclides sorbed to lagoon sediments and the exchange of radionuclides between the water and sediment. FVCOM-ICM-predicted radionuclide concentrations were then used to estimate radiation dose to humans and biota at all islands of the Enewetak Atoll. Under current climate conditions, annual radiation exposures for the southern islands including Enewetak (Fred) and Medren (Elmer) were below the current U.S. standards. Radiation doses were somewhat elevated starting at Runit Island northward and westward to Enjebi Island (Janet). The islands in the northwest quadrant, particularly Bokoluo (Alice) and Bokombako (Belle), remain relatively contaminated. The islands in the southwestern quadrant have low contamination. The highest contribution to radiation doses comes from consumption of locally grown foods. Two radionuclides, 90Sr and 137Cs, contributed the greatest fraction for most terrestrial foods. In current climate conditions, the storms temporarily increased radionuclide concentrations in the lagoon waters, increasing the radiation dose slightly. In future conditions, doses are expected to be smaller, primarily because of the radioactive decay of the shorter-lived radioisotopes of 90Sr and 137Cs. This could make all islands in the far northwest of the atoll – except Bokombako (Belle) and perhaps Bokoluo (Alice) – suitable for residency. For the f

Prasad, Rajiv↗

DOE OSTI2022

Off-of-equilibrium effects on Kurtosis Along Strangeness-Neutral Trajectories

The Beam Energy Scan program at the Relativistic Heavy Ion Collider (RHIC) is searching for the QCD critical point. The main signal for the critical point is the kurtosis of the distribution of proton yields obtained on an event by event basis where one expects a peak at the critical point. However, its exact behavior is still an open question due to out-of-equilibrium effects and uncertainty in the equation of state. Here we use a simplistic hydrodynamic model that enforces strangeness-neutrality, selecting trajectories that pass close to the critical point. We vary the initial conditions to estimate the effect of out-of-equilibrium hydrodynamics on the kurtosis signal.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS2004

An Ellipsoidal Particle-Finite Element Method for Hypervelocity Impact Simulation

A number of coupled particle-element and hybrid particle-element methods have been developed for the simulation of hypervelocity impact problems, to avoid certain disadvantages associated with the use of pure continuum based or pure particle based methods. To date these methods have employed spherical particles. In recent work a hybrid formulation has been extended to the ellipsoidal particle case. A model formulation approach based on Lagrange's equations, with particles entropies serving as generalized coordinates, avoids the angular momentum conservation problems which have been reported with ellipsoidal smooth particle hydrodynamics models.

Shivarama, Ravishankar↗

DOE OSTI2021

Skewness of mean transverse momentum fluctuations in heavy-ion collisions

We propose the skewness of mean transverse momentum, $\langle{p_t}\rangle$, fluctuations as a fine probe of hydrodynamic behavior in relativistic nuclear collisions. We describe how the skewness of the $\langle{p_t}\rangle$ distribution can be analyzed experimentally, and we use hydrodynamic simulations to predict its value. We predict in particular that $\langle{p_t}\rangle$ fluctuations have positive skew, which is significantly larger than if particles were emitted independently. We elucidate the origin of this result by deriving generic formulas relating the fluctuations of $\langle{p_t}\rangle$ to the fluctuations of the early-time thermodynamic quantities. We postulate that the large positive skewness of $\langle{p_t}\rangle$ fluctuations is a generic prediction of hydrodynamic models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1989

Pressure-confined Lyman-alpha clouds

Results are presented of numerical models of pressure-confined spherical gas clouds which produce absorption resembling the low to intermediate atomic column density lines found in high-redshift QSO spectra. One-dimensional hydrodynamical models including electron conduction are described, and the rate equations are solved to find ionization and excitation states. Results are presented for both static and adiabatically expanding confining media covering a range of initial pressures. It is found that Ly-alpha lines are very similar over a wide range of conditions and that the most promising diagnostic of pressure is to compare the column density in H I to that in He I and He II. No single-pressure model can explain the wide range of observed H I column densities.

Baron, E.↗

NASA NTRS1991

Physics of heavy ions (1989-1990)

The results from studies on polar wind ion heating due to kinetic ion beam instabilities and the effects of such ion heating on the outflow of O(+) in the polar wind are presented and discussed. First, the linear instabilities associated with an O(+) and H(+) polar wind plasma in the presence of O(+) and H(+) beams for a range of O(+)/H(+) beam densities, T(sub e)/T(sub i), and ion beam speeds were examined. Then, nonlinear heating of the polar wind ions was studied, using numerical simulations. The O(+) and H(+) polar wind ions were modeled by isotropic Maxwellian distributions, and the electrons, O(+) beams, and H(+) beams were modeled by drifting Maxwellian distributions. The effects of the kinetic ion heating on the outflow of the polar wind ions were examined from the ionosphere, using a time-dependent hydrodynamic model. A numerical code to solve the O(+) and H(+) continuity and momentum equations in a flux tube from ionospheric to magnetospheric altitudes were developed. The effects of ion heating were included by allowing for the altitudinal variation of the ion temperatures in the momentum equation. The ion temperature profiles were specified based on the ion heating characteristics found from previous kinetic simulations. It was assumed that heating occurred above 1500 km and increased to a saturated value of temperature that was obtained directly from the kinetic simulation study. The characteristics of the dynamical polar wind without ion heating were studied, and a flux tube on closed field lines that suddenly became open at t = 0 was simulated. Then, the effects of ion heating were included. To gain some physical insight, two limiting cases were considered: preferential H(+) heating and preferential O(+) heating. How O(+) heating can lead to enhanced polar wind O(+) fluxes in the polar magnetosphere is shown.

Ashour-Abdalla, M.↗

NASA NTRS1975

Shock waves in a beta Doradus model

The accumulation of observational evidence and interpretation supports the conclusion that beta Doradus has running shock waves in its atmosphere. It is the work by Hutchinson with its predictions for the phases of the shock waves that prompted this author to construct a numerical hydrodynamical model for the atmosphere of beta Doradus. A short description of the procedure of the calculation follows. The results for the beta Doradus model contain a comparison with the observations. It is of particular interest to verify that the atmosphere contains multiple shock waves as suggested by Hutchinson for the ultraviolet observations.

Hill, S. J.↗

NASA NTRS2014

A High-Resolution X-Ray and Optical Study of SN1006: Asymmetric Expansion and Small-Scale Structure in a Type Ia Supernova Remnant

We introduce a deep (670 ks) X-ray survey of the entire SN 1006 remnant from the Chandra X-Ray Observatory, together with a deep Hα image of SN 1006 from the 4 m Blanco telescope at CTIO. Comparison with Chandra images from 2003 gives the first measurement of the X-ray proper motions around the entire periphery, carried out over a 9 yr baseline. We find that the expansion velocity varies significantly with azimuth. The highest velocity of approx.7400 km/s (almost 2.5 times that in the northwest (NW)) is found along the southeast (SE) periphery, where both the kinematics and the spectra indicate that most of the X-ray emission stems from ejecta that have been decelerated little, if at all. Asymmetries in the distribution of ejecta are seen on a variety of spatial scales. Si-rich ejecta are especially prominent in the SE quadrant, while O and Mg are more uniformly distributed, indicating large-scale asymmetries arising from the explosion itself. Neon emission is strongest in a sharp filament just behind the primary shock along the NWrim, where the pre-shock density is highest. Here the Ne is likely interstellar, while Ne within the shell may include a contribution from ejecta. Within the interior of the projected shell we find a few isolated "bullets" of what appear to be supernova ejecta that are immediately preceded by bowshocks seen in Hα, features that we interpret as ejecta knots that have reached relatively dense regions of the surrounding interstellar medium, but that appear in the interior in projection. Recent three-dimensional hydrodynamic models for Type Ia supernovae display small-scale features that strongly resemble the ones seen in X-rays in SN 1006; an origin in the explosion itself or from subsequent hydrodynamic instabilities both remain viable options. We have expanded the search for precursor X-ray emission ahead of a synchrotron-dominated shock front, as expected from diffusive shock acceleration theory, to numerous regions along both the northeast and southwest rims of the shell. Our data require that a precursor be thinner than about 3, and fainter than about 5% of the post-shock peak. These limits suggest that the magnetic field is amplified by a factor of seven or more in a narrow precursor region, promoting diffusive particle acceleration.

galaxies: active↗

DOE OSTI2025

Supercooled Goldstone Bosons at the QCD Chiral Phase Transition

We discuss a universal nonequilibrium enhancement of long-wavelength Goldstone bosons induced by quenches to the broken phase in Model G—the dynamical universality class of an 𝑂⁡(4) antiferromagnet and the chiral phase transition in QCD. Scaling arguments for the coarsening dynamics describing the formation of the chiral condensate predict a parametric enhancement in the infrared spectra of Goldstone bosons, a prediction confirmed by stochastic simulations of the transition. The details of the enhancement are determined by the nonlinear dynamics of a superfluid effective theory, which is a limit of Model G reflecting the broken 𝑂⁡(4) symmetry. Our results translate to a parametric enhancement of low-momentum pions in heavy-ion collisions at the LHC, which are underpredicted in current hydrodynamic models without critical dynamics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS2000

Modeling of dense reactive granular flows

Because of the lack of fundamental knowledge of both the fluid mechanical and thermal behavior of multiphase flow systems, the design of these processes has been necessarily based on empirical correlations and experiments performed in laboratory or pilot scale units. Detailed models and simulations of these systems would allow for the optimization of the processes and their scaling, without the need for expensive and time consuming testing. In this paper, we present the proposed hydrodynamic model in detail.

granular↗

DOE OSTI2022

Self-mixing in microtubule-kinesin active fluid from nonuniform to uniform distribution of activity

Active fluids have applications in micromixing, but little is known about the mixing kinematics of systems with spatiotemporally-varying activity. To investigate, UV-activated caged ATP is used to activate controlled regions of microtubule-kinesin active fluid and the mixing process is observed with fluorescent tracers and molecular dyes. At low Péclet numbers (diffusive transport), the active-inactive interface progresses toward the inactive area in a diffusion-like manner that is described by a simple model combining diffusion with Michaelis-Menten kinetics. At high Péclet numbers (convective transport), the active-inactive interface progresses in a superdiffusion-like manner that is qualitatively captured by an active-fluid hydrodynamic model coupled to ATP transport. Results show that active fluid mixing involves complex coupling between distribution of active stress and active transport of ATP and reduces mixing time for suspended components with decreased impact of initial component distribution. This work will inform application of active fluids to promote micromixing in microfluidic devices.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Modeling Astrophysical Explosions and Nucleosynthesis (Final Technical Report)

The goal of this project was to further our theoretical understanding of the extreme physics and heavy element nucleosynthesis occurring in stellar explosions, in particular supernovae and neutron star mergers. To study these systems, we further developed and applied numerical simulations (radiation hydrodynamical models and nuclear reaction networks) run on high performance computers. The numerical calculations and companion theoretical analysis helped to better define the sites of astrophysical nucleosynthesis, and make predictions of the associated observables (nucleosynthetic yields and electromagnetic signals) that can be used to anticipate and interpret experimental data. This project enhances the scientific value of the nuclear data obtained from current and future experimental facilities (e.g., rare isotope beams) by improving our understanding of the astrophysical context and specific physical conditions under which we expect such reactions to take place in Nature.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRSDate not supplied

Constraining Nanoflare Frequency Through Observations and Models

Nanoflares are considered prime candidates for heating the non-flaring solar corona. However, direct detection of individual nanoflares with present-generation instruments remains challenging. Understanding the frequency and magnitude of nanoflares is crucial for comprehending their role in coronal heating. In this study, we employ field-aligned hydrodynamic model to simulate the emission from a non-flaring X-ray bright point (XBP) observed by the Marshall Grazing Incidence X-ray Spectrometer (MaGIXS). The length and magnetic field strength of the coronal loops associated with this XBP are derived from magnetic field extrapolation of observed photospheric magnetograms by SDO/HMI. Each loop is assumed to be heated by random nanoflares, whose frequency and magnitude are determined by the loop length and their field strength. The simulation results are then compared and matched against the observation to determine the nanoflare heating frequency. We will discuss the findings of this study. Furthermore, we will discuss the evolution of nanoflare heating frequency during the formation and decay of an active region (AR) on the solar disk.

coronalheating↗

DOE OSTI2021

2021 Summer SPE Project

SPE is a project in order to develop new, more physics-based, seismic models of explosions (see Nelson et al. 2013). One key component of this effort is numerical modeling enabled by modern state-of-art code/software. Accurate modeling of the shape and amplitude of seismic waves from their generation to their propagation to remote monitoring seismic stations is important to our ability to determine the origin and strength of the source from remote recording. In this project, the modeling is performed by coupling two codes HOSS and SPECFEM3D. HOSS models the dynamic nonlinear processes happening near the explosion. SPECFEM3D computes the propagation of seismic waves as they travel through 3D complex Earth models. SPECFEM3D was modified in order to be driven by a set of time-series calculated by HOSS in lieu of a seismic source. The goal of this summer project is the investigation of several questions pertinent to the establishment of a full end-to-end modeling capability from the high-rate strain regime area to remote distances where seismic station record seismic waves generated by explosions. The investigated questions are: (1) How to perform proper filtering? Direct-solution modeling only sustains a limited range of frequency depending on the grid size. As we go from one modeling domain to the other via coupling, the mesh size is getting coarser to allow modeling at large scale but also to account for the fact that high-frequency waves do not physically travel to large distances. So filtering of the time-series generated by the near-field hydrodynamic modeling is a current practice often employed but its effects on the modeled waveforms has to be investigated. (2) Quantitative assessment of the efficiency of attenuation to remove high-frequency content of the wavefield. This assessment will allow to create meshes with a grid size appropriate to the actual physics of wave propagation for a given explosion. (3) Checking that the coupling process and the two codes respect the cylindrical symmetry that is expected in the case of a pure explosion in a half-space. (4) The effect of the state-of-stress in the near-source area on the modeled seismic waveforms. These questions will be investigated through the modeling of SPE-4P, the fourth explosion of this series because it was designed to have little interaction with the geologic surrounding and the free surface of the Earth so that it is the most explosion-like experiment, with the most symmetries to be verified.

58 GEOSCIENCES↗

NASA NTRS1974

Laser Doppler velocimeter system simulation for sensing aircraft wake vortices

A hydrodynamic model of aircraft vortex wakes in an irregular wind shear field near the ground is developed and used as a basis for modeling the characteristics of a laser Doppler detection and vortex location system. The trailing vortex sheet and the wind shear are represented by discrete free vortices distributed over a two-dimensional grid. The time dependent hydrodynamic equations are solved by direct numerical integration in the Boussinesq approximation. The ground boundary is simulated by images, and fast Fourier Transform techniques are used to evaluate the vorticity stream function. The atmospheric turbulence was simulated by constructing specific realizations at time equal to zero, assuming that Kolmogoroff's law applies, and that the dissipation rate is constant throughout the flow field. The response of a simulated laser Doppler velocimeter is analyzed by simulating the signal return from the flow field as sensed by a simulation of the optical/electronic system.

Thomson, J. A. L.↗

DOE OSTI2021

Observation of ionization trends in a laboratory photoionized plasma experiment at Z

We report experimental and modeling results for the charge state distribution of laboratory photoionized neon plasmas in the first systematic study over nearly an order of magnitude range of ionization parameter ξ ∝ F/N e . The range of ξ is achieved by flexibility in the experimental platform to adjust either the x-ray drive flux F at the sample or the electron number density N e or both. Experimental measurements of photoionized plasma conditions over such a range of parameters enable a stringent test of atomic kinetics models used within codes that are applied to photoionized plasmas in the laboratory and astrophysics. From experimental transmission data, ion areal densities are extracted by spectroscopic analysis that is independent of atomic kinetics modeling. The measurements reveal the net result of the competition between photon-driven ionization and electron-driven recombination atomic processes as a function of ξ as it affects the charge state distribution. Results from radiation-hydrodynamics modeling calculations with detailed inline atomic kinetics modeling are compared with the experimental results. There is good agreement in the mean charge and overall qualitative similarities in the trends observed with ξ but significant quantitative differences in the fractional populations of individual ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗