Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrodynamics simulations”

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,801 records · Page 100

Investigation of the Role of the Role of Nuclear Physics in Heavy Element Nucleosynthesis, through the Study of Key Reactions, and the Improvement of Theoretical Reaction Rates (Final Report)

Neutrino-driven winds in core-collapse supernovae have been identified as a possible site for the production of elements heavier than iron. Traditionally, these neutrino-driven winds have been proposed as the site of the main r-process. Recent simulations fail to reproduce the conditions required for the main r-process. while they remain a promising site for producing the lightest elements beyond iron, e.g., Sr, Y, and Zr through the νp process. The efficiency of the νp process depends on the hydrodynamical conditions, the electron fraction (which is related to the neutrino properties), and the nuclear reactions on many short-lived nuclei with limited (if at all) experimental information. The reaction rates on these nuclei are based on theoretical predictions using the Hauser-Feshbach model. Recent sensitivity studies have highlighted the importance of neutron-induced reactions on these nuclei along the νp process path. This work aimed to experimentally constrain reaction rates that are known to play a key role in the neutrino-p process nucleosynthesis. A secondary subsequently-added objective was to start the implementation of techniques that improve the description of nuclear properties in the Hauser-Feshbach model by extending the microscopic nuclear level density description offered via the shell model to high excitation energies without using experiment-based renormalizations. The main objective of this work was the experimental constraint of the 56 Ni(n,p) 56 Co reaction rate via a measurement of the inverse reaction 56 Co(p,n) 56 Ni at the National Superconducting Cyclotron Laboratory (NSCL) and later the Facility for Rare Isotope Beams (FRIB). This reaction is considered the key one for determining the yields possible by the neutrino-p process. A technique for this type of measurement in inverse kinematics at low energies did not exist before this work. The work also had two secondary objectives. First, to contribute to efforts to measure the same reaction in direct kinematics using a radioactive target at Los Alamos National Laboratory (LANL), and second, to advance work to implement shell-model-deduced microscopic level densities in Hauser-Feshbach calculations. The project has resulted in the development of the first technique to perform (p,n) cross-section measurements in relevant-for-astrophysics low energies in inverse kinematics using a magnetic spectrometer or separator, and neutron detectors for neutron tagging. It has also resulted in the precise measurement of the cross-section of the 40 Ar(p,n) 40 K reaction in a proof-of-principle experiment realized by using a beam-line quadrupole of the ReA3 accelerator of NSCL/FRIB. As part of this project the technique was successfully adapted to make use of the superior acceptance of the Separator for Capture Reactions (SECAR) at FRIB. In this project, the required experimental setup simulations and beam optics were developed and tested with the measurement of the 58 Fe(p,n) 58 Cu reaction cross-section. Additionally, this project contributed with simulation work to the development of a technique to measure (n,p) reactions with radioactive targets at LANL, and which resulted in the measurement of the key 56 Ni(n,p) 56 Co reaction cross-section at neutron energies above ≈1 MeV. Last, the project initiated work in the development of shell model based level densities using the moments method.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

OC6 Phase Ia: CFD Simulations of the Free-Decay Motion of the DeepCwind Semisubmersible

Currently, the design of floating offshore wind systems is primarily based on mid-fidelity models with empirical drag forces. The tuning of the model coefficients requires data from either experiments or high-fidelity simulations. As part of the OC6 (Offshore Code Comparison Collaboration, Continued, with Correlation, and unCertainty (OC6) is a project under the International Energy Agency Wind Task 30 framework) project, the present investigation explores the latter option. A verification and validation study of computational fluid dynamics (CFD) models of the DeepCwind semisubmersible undergoing free-decay motion is performed. Several institutions provided CFD results for validation against the OC6 experimental campaign. The objective is to evaluate whether the CFD setups of the participants can provide valid estimates of the hydrodynamic damping coefficients needed by mid-fidelity models. The linear and quadratic damping coefficients and the equivalent damping ratio are chosen as metrics for validation. Large numerical uncertainties are estimated for the linear and quadratic damping coefficients; however, the equivalent damping ratios are more consistently predicted with lower uncertainty. Some difference is observed between the experimental and CFD surge-decay motion, which is caused by mechanical damping not considered in the simulations that likely originated from the mooring setup, including a Coulomb-friction-type force. Overall, the simulations and the experiment show reasonable agreement, thus demonstrating the feasibility of using CFD simulations to tune mid-fidelity models.

17 WIND ENERGY↗

Second release of the CoRe database of binary neutron star merger waveforms

Abstract We present the second data release of gravitational waveforms from binary neutron star (BNS) merger simulations performed by the Computational Relativity ( CoRe ) collaboration. The current database consists of 254 different BNS configurations and a total of 590 individual numerical-relativity simulations using various grid resolutions. The released waveform data contain the strain and the Weyl curvature multipoles up to ℓ = m = 4 . They span a significant portion of the mass, mass-ratio, spin and eccentricity parameter space and include targeted configurations to the events GW170817 and GW190425. CoRe simulations are performed with 18 different equations of state, seven of which are finite temperature models, and three of which account for non-hadronic degrees of freedom. About half of the released data are computed with high-order hydrodynamics schemes for tens of orbits to merger; the other half is computed with advanced microphysics. We showcase a standard waveform error analysis and discuss the accuracy of the database in terms of faithfulness. We present ready-to-use fitting formulas for equation of state-insensitive relations at merger (e.g. merger frequency), luminosity peak, and post-merger spectrum.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

O 16 O 16 collisions at energies available at the BNL Relativistic Heavy Ion Collider and at the CERN Large Hadron Collider comparing α clustering versus substructure

Collisions of . light and heavy nuclei in relativistic heavy-ion collisions have been shown to be sensitive to nuclear structure. With a proposed 16 O 16 O run at the CERN Large Hadron Collider (LHC) and at the BNL Relativistic Heavy Ion Collider (RHIC) we study the potential for finding α clustering in 16 O. Here we use the state-of-the-art iEBE-VISHNU package with 16 O nucleonic configurations from ab initio nuclear lattice simulations. This setup was tuned using a Bayesian analysis on p Pb and PbPb systems. We find that the 16 O 16 O system always begins far from equilibrium and that at LHC and RHIC it approaches the regime of hydrodynamic applicability only at very late times. Finally, by taking ratios of flow harmonics we are able to find measurable differences between α-clustering, nucleonic, and subnucleonic degrees of freedom in the initial state.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Validation Framework for Post-Combustion Carbon Capture CFD Simulations

First-principles based computational fluid dynamics (CFD) simulations are proposed as a fundamental tool for investigating solvent-based CO2 absorption in packed columns, due to the ability to accurately represent the underlying non-linear multiscale dynamics. In this work, we employ such models to investigate hydrodynamics of columns with structured by assessing the key hydrodynamic metrics, such as pressure drop and liquid holdup. Our models are validated with experimental data from a specifically designed column for this line of work. The test cases of gas and liquid flowrates and operating conditions were selected through a comprehensive sequential design of experiments approach offered by the CCSI2 toolset.

Panagakos, Grigorios↗

Accretion in wind-driven X-ray sources

The hydrodynamics of axisymmetric accretion flow in stellar wind-fed X-ray sources is investigated, including momentum deposition by radiation, as well as radiative heating and cooling by Compton and bremsstrahlung processes. The results of two-dimensional numerical simulations for mass accretion rates ranging from 0.02 to 0.4 of the Eddington value reveal a variety of behaviors. At low rates radiative effects are unimportant, and the accretion flow is nearly adiabatic. In this case the flow is steady. For intermediate accretion rates, radiative heating and cooling effects become important. The flow remains steady, but the accretion rate is significantly reduced from the Hoyle-Lyttleton estimate. At the highest mass accretion rates, radiative momentum transfer and energy exchange are important and lead to nonsteady behavior. For the parameters relevant to massive X-ray binary systems accreting at high rates, the results reveal that the accretion flows can be complex and time dependent.

Taam, Ronald E.↗

The Micromechanics of the Moving Contact Line

A transient moving contact line is investigated experimentally. The dynamic interface shape between 20 and 800 microns from the contact line is compared with theory. A novel experiment is devised, in which the contact line is set into motion by electrically altering the solid-liquid surface tension gamma(sub SL). The contact line motion simulates that of spontaneous wetting along a vertical plate with a maximum capillary number Ca approx. = 4 x 10(exp -2). The images of the dynamic meniscus are analyzed as a funtion of Ca. For comparison, the steady-state hydrodynamic equation based on the creeping flow model in a wedge geometry and the three-region uniform perturbation expansion of Cox (1986) is adopted. The interface shape is well depicted by the uniform solutions for Ca <= 10(exp -3). However, for Ca > 10(exp -3), the uniform solution over-predicts the viscous bending. This over-prediction can be accounted for by modifying the slip coefficient within the intermediate solution. With this correction, the measured interface shape is seen to match the theoretical prediction for all capillary numbers. The amount of slip needed to fit the measurements does not scale with the capillary number.

Lichter, Seth↗

Safety Applications Project FY20-21 Highlights, FY2022 Goals [Slides]

Collaborations leverage the investment and expertise of others. Code capabilities determine the questions that can be addressed using SAP tools. Infrastructure acts as a force multiplier for software developer skills. Improved code performance allows simulations to finish faster or with greater fidelity. Code performance has dramatically improved over the past ~4 years.

42 ENGINEERING↗

Use of hydrodynamic theory to estimate electrical current redistribution in metals

Using the analogy between hydrodynamic and electrical current flow, we study herein how electrical current density j redistributes and amplifies due to two commonly encountered inhomogeneities in metals. First, we consider flow around a spherical resistive inclusion and find significant j amplification, independent of inclusion size. Hence, even μ m-scale inclusions can affect performance in applications by creating localized regions of enhanced Joule heating. Next, we investigate j redistribution due to surface roughness, idealized as a sinusoidal perturbation with amplitude A and wavelength λ . Theory predicts that j amplification is determined by the ratio A / λ , so that even “smooth” surface finishes (i.e., small A ) can generate significant amplification, if λ is correspondingly small. We compare theory with magnetohydrodynamic simulation to illustrate both the utility and limitations of the steady-state theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

X-ray spectroscopy evidence for plasma shell formation in experiments modeling accretion columns in young stars

Recent achievements in laboratory astrophysics experiments with high-power lasers have allowed progress in our understanding of the early stages of star formation. In particular, we have recently demonstrated the possibility of simulating in the laboratory the process of the accretion of matter on young stars [G. Revet et al., Sci. Adv. 3, e1700982 (2017)]. The present paper focuses on x-ray spectroscopy methods that allow us to investigate the complex plasma hydrodynamics involved in such experiments. We demonstrate that we can infer the formation of a plasma shell, surrounding the accretion column at the location of impact with the stellar surface, and thus resolve the present discrepancies between mass accretion rates derived from x-ray and optical-radiation astronomical observations originating from the same object. In our experiments, the accretion column is modeled by having a collimated narrow (1 mm diameter) plasma stream first propagate along the lines of a large-scale external magnetic field and then impact onto an obstacle, mimicking the high-density region of the stellar chromosphere. A combined approach using steady-state and quasi-stationary models was successfully applied to measure the parameters of the plasma all along its propagation, at the impact site, and in the structure surrounding the impact region. The formation of a hot plasma shell, surrounding the denser and colder core, formed by the incoming stream of matter is observed near the obstacle using x-ray spatially resolved spectroscopy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Evidence of non-Maxwellian ion velocity distributions in spherical shock-driven implosions

The ion velocity distribution functions of thermonuclear plasmas generated by spherical laser direct drive implosions are studied using deuterium-tritium (DT) and deuterium-deuterium (DD) fusion neutron energy spectrum measurements. A hydrodynamic Maxwellian plasma model accurately describes measurements made from lower temperature (< 10 keV), hydrodynamiclike plasmas, but is insufficient to describe measurements made from higher temperature more kineticlike plasmas. The high temperature measurements are more consistent with Vlasov-Fokker-Planck (VFP) simulation results which predict the presence of a bimodal plasma ion velocity distribution near peak neutron production. Furthermore, these measurements provide direct experimental evidence of non-Maxwellian ion velocity distributions in spherical shock driven implosions and provide useful data for benchmarking kinetic VFP simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

TEAMER: Sandia and CalWave Torque Spring Assessment

Power-take-off (PTO) systems for wave energy converters (WEC) require restoring forces for efficient hydrodynamic as well as mechanical to electric power transfer. Implementation of an effective spring mechanism that can provide a restoring torque on e.g. a rotational, winch type PTO shaft is challenging. SNL in collaboration with CalWave, will use numerical finite element simulations to study application of mechanical torsion springs as a highly efficient restoring force element for CalWave's rotational PTO. Material studies will include common steel as a baseline and potentially advanced materials such as fiber composite materials with different cross-sectional profiles. Test report including description and results of FEA (finite element analysis) models. Excel sheet summarizes which scripts were used for which tables and figures in the report. Part of TEAMER request for technical support (RFTS 1) award to CalWave and Sandia National Lab (SNL).

16 TIDAL AND WAVE POWER↗

Raytracing

The raytracing (spectrally specific radiative transfer that does not couple back to the hydrodynamic evolution) code was developed because it helps improve our understanding of the physics of explosive events, validates computational physics codes, and estimates optical signals that are cheaper and safer than experiments involving high explosives. The raytracing code can also produce simulations that are useful to analyze what we can understand with sensors, such as determining what we can discriminate with sensors, like the presence of certain metals or the amount of soot. In addition, simulations can provide extra validations to computational physics to help us better understand discrepancies.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Crosslink V.0.11.x User Manual

CrossLink is a novel two-dimensional and three-dimensional geometry and mesh generation software package developed by the Simulation Tools team at Los Alamos National Laboratory. This software represents the third generation of topology-based mesh generation technology developed by the Department of Defense and the Department of Energy with a special focus on complex multi-material hydrodynamic applications, mesh scalability, and high-order element mesh generation. The topology-based meshing approach offered by CrossLink enables users to quickly and easily mesh complex geometries in a repeatable and robust manner. CrossLink’s topology-based meshing approach is well-suited for parametric design studies, parametric design optimization, damage scenario assessment, and iterative design modification (i.e. feature addition and/or removal). CrossLink’s python API allows workflow scripting of the geometry creation and mesh generation process for traceability, repeatability, data provenance, and version control. CrossLink consists of three main components: a graphical user interface (GUI), a geometry creation and mesh generation engine, and a python API that provides a workflow scripting interface to the geometry and meshing functions.

97 MATHEMATICS AND COMPUTING↗

National direct-drive program on OMEGA and the National Ignition Facility

A major advantage of the laser direct-drive (DD) approach to ignition is the increased fraction of laser drive energy coupled to the hot spot and relaxed hot-spot requirements for the peak pressure and convergence ratios relative to the indirect-drive approach at equivalent laser energy. With the goal of a successful ignition demonstration using DD, the recently established national strategy has several elements and involves multiple national and international institutions. These elements include the experimental demonstration on OMEGA cryogenic implosions of hot-spot conditions relevant for ignition at MJ-scale energies available at the National Ignition Facility (NIF) and developing an understanding of laser-plasma interactions and laser coupling using DD experiments on the NIF. DD designs require reaching central stagnation pressures in excess of 100 Gbar. The current experiments on OMEGA have achieved inferred peak pressures of 56 Gbar (Regan et al 2016 Phys. Rev. Lett. 117 025001). Extensive analysis of the cryogenic target experiments and two- and three-dimensional simulations suggest that power balance, target offset, and target quality are the main limiting factors in target performance. In addition, cross-beam energy transfer (CBET) has been identified as the main mechanism reducing laser coupling. Reaching the goal of demonstrating hydrodynamic equivalence on OMEGA includes improving laser power balance, target position, and target quality at shot time. CBET must also be significantly reduced and several strategies have been identified to address this issue. (2016-01-01)

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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.↗

The simulation of the alternate turbopump development high pressure oxygen and fuel turbopumps for the space shuttle main engine using the Shaberth computer program

The Space Shuttle Main Engine (SSME) is basically comprised of a combustion chamber and nozzle, high and low pressure oxygen turbopumps and high and low pressure fuel turbopumps. In the current configuration, the high pressure fuel (HPTFP) and high pressure oxygen turbopumps (HPOTP) have experienced a history of ball bearing wear. The wear problem can be attributed to numerous factors including the hydrodynamic axial and radial loads caused by the flow of liquid oxygen and liquid hydrogen through the turbopump impellers and turbine. Also, friction effects between the rolling elements, races, and cage can create thermally induced bearing geometry changes. To alleviate some of the current configuration problems, an alternate turbopump development (ATD) was proposed. However, the ATD HPOTP and HPTFP are constrained to operate interchangeably with the current turbopumps, thus, the operation conditions must be similar. The ATD configuration features a major change in bearings used to support the integrated shaft, impeller, and turbine system. A single ball and single roller will replace the pump-end and turbine and duplex ball bearings. The Shaft-Bearing-Thermal (SHABERTH) computer code was used to model the ATD HPOTP and ATD HPFTP configurations. A two bearing model was used to simulate the HPOTP and HPFTP bearings and shaft geometry. From SHABERTH, a comparison of bearing reaction loads, frictional heat generation rates, and Hertz contact stresses will be attempted with analysis at the 109 percent and 65 percent power levels.

Mcdonald, Gary H.↗

One-sweep moment-based semi-implicit-explicit integration for gray thermal radiation transport

Thermal radiation transport (TRT) is a time dependent, high dimensional partial integro-differential equation. In practical applications such as inertial confinement fusion, TRT is coupled to other physics such as hydrodynamics, plasmas, etc., and the timescales one is interested in capturing are often much slower than the radiation timescale. As a result, TRT is treated implicitly, and due to its stiffness and high dimensionality, is often a dominant computational cost in multiphysics simulations. Here we develop a new approach for implicit-explicit (IMEX) integration of gray TRT in the deterministic SN setting, which requires only one sweep per stage, with the simplest first-order method requiring only one sweep per time step. The partitioning of equations is done via a moment-based high-order low-order formulation of TRT, where the streaming operator and first two moments are used to capture the asymptotic stiff regimes of the streaming limit and diffusion limit. Absorption-reemission is treated explicitly, and although stiff, is sufficiently damped by the implicit solve that we achieve stable accurate time integration without incorporating the coupling of the high order and low order equations implicitly. Due to nonlinear coupling of the high-order and low-order equations through temperature-dependent opacities, to facilitate IMEX partitioning and higher-order methods, we use a semi-implicit integration approach amenable to nonlinear partitions. In conclusion, results are demonstrated on thick Marshak and crooked pipe benchmark problems, demonstrating orders of magnitude improvement in accuracy and wallclock compared with the standard first-order implicit integration typically used.

97 MATHEMATICS AND COMPUTING↗

Refine your search to explore more results.