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 955 records · Page 53

On the mass of the Cepheid X Cygni

Using a new non-Lagrangian hydrodynamic radiation transport code, models of the Cepheid X Cygni have been calculated and compared with new observations. Adopting an effective temperature of 5300 K, the luminosity was varied for stellar masses of 6, 8, and 9 solar masses in order to obtain a period of 16.4d. The results indicate that in order to find reasonable agreement between the light curves as calculated and observed, the mass used in the models must be nearly the evolutionary mass, not some 60% less as found in earlier work for Cepheids with periods from 7d to 10d. In addition, the new observations support the physical reality of the so-called 'artificial viscosity dip'.

Davis, C. G.↗

Statistical mechanics of dilute liquid mixtures of He3 in He4

The Woo-Tan-Massey quantum-mechanical theory of dilute mixtures of He-3 in liquid He-4 is extended to incorporate superfluid flow. Then statistical-mechanical principles are employed to find explicit formulas for thermodyanmic and hydrodynamic properties of mixtures, such as internal energy, free energy, entropy, pressure, chemical potential, and normal-fluid and superfluid densities. This paper is the first in a planned series that applies methods originally developed for pure liquid He-4 to dilute mixtures as a means of establishing a comprehensive theory that includes a theoretical bridge between a microscopic description of the mixtures and a hydrodynamic two-fluid model. Ultimately effects of applied electric and magnetic fields on dilute solutions are also planned to be treated.

Jackson, H. W.↗

On radial oscillations in viscous accretion discs surrounding neutron stars

Radial oscillations resulting from axisymmetric perturbations in viscous accretion disks surrounding neutron stars in X-ray binary systems have been investigated. Within the framework of the alpha-viscosity model a series of hydrodynamic calculations demonstrates that the oscillations are global for alpha of about 1. On the other hand, for alpha of 0.4 or less, the oscillations are local and confined to the disk boundaries. If viscous stresses acting in the radial direction are included, however, it is found that the disk can be stabilized. The application of such instabilities in accretion disks, without reference to the boundary layer region between the neutron star (or magnetosphere) and the inner edge of the disk, to the phenomenology of quasi-periodic oscillations is brought into question.

Chen, Xingming↗

Enrichment and heating of the intracluster medium by ejection from galaxies

Results of N-body + hydrodynamic simulations designed to model the formation and evolution of clusters of galaxies and intracluster gas are presented. Clusters of galaxies are the largest bound, relaxed objects in the universe. They are strong x-ray emitters; this radiation originates through thermal bremsstrahlung from a diffuse plasma filling the space between cluster galaxies, the intracluster medium or ICM. From observations, one can infer that the mass of the ICM is comparable to or greater than the mass of all the galaxies in the cluster, and that the ratio of mass in hot gas to mass in galaxies, M(sub ICM)/M(sub STARS), increases with the richness of the cluster. Spectroscopic studies of cluster x-ray emission show heavy element emission lines. While M(sub ICM)/M(sub STARS) is greater than or equal to 1 implies that most of the ICM is primordial in nature, the discovery of heavy elements indicates that some of the gas must have been processed through galaxies. Galaxy evolution thus directly impacts cluster evolution.

Metzler, Chris↗

A method of smoothed particle hydrodynamics using spheroidal kernels

We present a new method of three-dimensional smoothed particle hydrodynamics (SPH) designed to model systems dominated by deformation along a preferential axis. These systems cause severe problems for SPH codes using spherical kernels, which are best suited for modeling systems which retain rough spherical symmetry. Our method allows the smoothing length in the direction of the deformation to evolve independently of the smoothing length in the perpendicular plane, resulting in a kernel with a spheroidal shape. As a result the spatial resolution in the direction of deformation is significantly improved. As a test case we present the one-dimensional homologous collapse of a zero-temperature, uniform-density cloud, which serves to demonstrate the advantages of spheroidal kernels. We also present new results on the problem of the tidal disruption of a star by a massive black hole.

Fulbright, Michael S.↗

System and Method for Modeling the Flow Performance Features of an Object

The method and apparatus includes a neural network for generating a model of an object in a wind tunnel from performance data on the object. The network is trained from test input signals (e.g., leading edge flap position, trailing edge flap position, angle of attack, and other geometric configurations, and power settings) and test output signals (e.g., lift, drag, pitching moment, or other performance features). In one embodiment, the neural network training method employs a modified Levenberg-Marquardt optimization technique. The model can be generated 'real time' as wind tunnel testing proceeds. Once trained, the model is used to estimate performance features associated with the aircraft given geometric configuration and/or power setting input. The invention can also be applied in other similar static flow modeling applications in aerodynamics, hydrodynamics, fluid dynamics, and other such disciplines. For example, the static testing of cars, sails, and foils, propellers, keels, rudders, turbines, fins, and the like, in a wind tunnel, water trough, or other flowing medium.

Jorgensen, Charles↗

Simulation of Hypervelocity Impact on Aluminum-Nextel-Kevlar Orbital Debris Shields

An improved hybrid particle-finite element method has been developed for hypervelocity impact simulation. The method combines the general contact-impact capabilities of particle codes with the true Lagrangian kinematics of large strain finite element formulations. Unlike some alternative schemes which couple Lagrangian finite element models with smooth particle hydrodynamics, the present formulation makes no use of slidelines or penalty forces. The method has been implemented in a parallel, three dimensional computer code. Simulations of three dimensional orbital debris impact problems using this parallel hybrid particle-finite element code, show good agreement with experiment and good speedup in parallel computation. The simulations included single and multi-plate shields as well as aluminum and composite shielding materials. at an impact velocity of eleven kilometers per second.

Fahrenthold, Eric P.↗

NED_Asteroid_Energy_Deposition

In the event of a potentially catastrophic asteroid impact, with sufficient warning time, deploying a nuclear device remains a powerful option for planetary defense if a kinetic impactor proves insufficient. Predicting the effectiveness of a potential nuclear deflection or disruption mission depends on accurate multiphysics simulations of the device's x-ray energy deposition into the asteroid and the resulting material ablation. These simulations span many orders of magnitude, require a variety of different complex physics packages, and are computationally expensive. Having an efficient and accurate way of modeling this system is necessary for exploring a mission's sensitivity to the asteroid's range of physical properties. To expedite future simulations, we present a completed x-ray energy deposition model developed using the radiation-hydrodynamics code Kull which can be used to initiate a nuclear mitigation mission hydrocode. The model spans a wide variety of possible mission initial conditions: four different asteroid-like materials (Silicon Dioxide, Forsterite, Iron, and Ice), two different source spectra (1 and 2 keV blackbodies), and then a broad range of radiation fluences (0.0001 to 1 kt per square meter), source durations (10 to 100 ns), and asteroid porosities (0 to 80 percent). Using blowoff momentum as the primary metric, the modelinitiated simulation results match the full radiation-hydrodynamics results to within 10 percent. Please reference the journal article: Burkey et al., X-Ray Energy Deposition Model for Simulating Asteroid Response to a Nuclear Planetary Defense Mitigation Mission, Planetary Science Journal, (2023) for more information.

Managan, RobertA↗

Verification study of xRAGE’s multi-ion viscosity model

This report on the adequacy of the plasma viscosity model implemented in LANL’s ASC program code, xRAGE is organized as follows. Section I provides a brief introduction on the kinetic theory which underpins plasma viscosity, and Section II describes how plasma kinetic theory can be used to calculate plasma viscosity coefficients accurate within the hydrodynamic limit. Section III describes the existing model for calculating the plasma viscosity coefficient in xRAGE. Next, Section IV shows tests of the xRAGE viscosity model against established hydrodynamic theory in the case of binary and trinary plasma mixtures.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Experimental platforms for investigating feature-driven jets for HED mix model validation

High-energy-density (HED) systems, such as inertial confinement fusion (ICF), are susceptible to hydrodynamic instabilities that can significantly affect both experimental results and modeling predictions. Isolated features, such as fill tubes or divots in the capsule, can cause material to jet as a result of the compressive shock exciting the Richtmyer–Meshkov instability, and serve as one of the primary degradation mechanisms in ICF yield. Simulations of feature-driven jets and how they mix require extensive experimental validation, particularly for understanding to what degree the initial size and shape of a feature influence jet dynamics, and how much instability feeds through downstream layers. A better understanding of feature-driven jetting can improve our mix modeling capabilities and increase hydrodynamic simulation accuracy. This manuscript describes a series of experimental platforms fielded by Los Alamos National Laboratory as a part of the Mshock Omega 60 and ModCons Omega EP campaigns to explore feature-driven jetting. These platforms are designed to benchmark jet evolution and growth as a function of initial feature size and shape, investigate jet-layer interactions leading to instability feedthrough, and will be used to characterize jet-jet interactions resulting from clusters of features. In conclusion, preliminary results for both platforms are shown. The ModCons experiments are on-going, and a discussion of future work directions is included.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Hard Life of Floatovoltaics: Modeling Wind-Driven Oscillations of Floating Solar Panels

Modern, thin photovoltaic (PV) panels for solar power are susceptible to high stress loads in windy conditions. Manufacturers are eager to determine optimal installation practices to reduce these loads, including in the relatively new practice of installing PV on floating structures located on artificial and natural lakes. As part of our effort to develop simulation capabilities for stress on such floating PV systems, we here present our combined model-simulation approach, which simulates the dynamic wind loading and uses modeled elements to capture both the hydrodynamic and mooring-line forces. We discuss the forces important for our model and the challenges inherent in our simulation. Of particular interest is our model of the hysteresis response displayed by the attached mooring lines, which are engineered to damp motion and oscillation of the floating panel system. We validate our approach against the benchmark problem of vortex-induced vibration of a cylinder, which is driven by the same dynamic forces present in floating PV panels.

floating↗

Blast from the past: constraining progenitor models of SN 1972E

ABSTRACT We present a novel technique to study Type Ia supernovae (SNe Ia) by constraining surviving companions of historical extragalactic SN by combining archival photographic plates and Hubble Space Telescope(HST) imaging. We demonstrate this technique for Supernova 1972E, the nearest known SN Ia in 125 yr. Some models of SNe Ia describe a white dwarf with a non-degenerate companion that donates enough mass to trigger thermonuclear detonation. Hydrodynamic simulations and stellar evolution models show that these donor stars will survive the explosion, and show increased luminosity for at least a 1000 yr. Thus, late-time observations of the exact location of a supernova can constrain the presence of a surviving donor star and progenitor models. We find the explosion site of SN 1972E by analysing 17 digitized photographic plates taken with the European Southern Observatory 1-m Schmidt and 1 plate taken with the Cerro Tololo Inter-American Observatory 1.5-m telescope. Using the Gaia eDR3 catalogue to determine Supernova 1972E’s location yields: α = 13h39m52${_{.}^{\rm s}}$708 ± 0${_{.}^{\rm s}}$004 and δ = −31°40’9${_{.}^{\prime\prime}}$00 ± 0${_{.}^{\prime\prime}}$04 (ICRS). In 2005, HST/ACS imaged the host galaxy of SN 1972E with the F435W, F555W, and F814W filters covering the explosion site. The nearest detected source is offset by 3.0 times our positional precision, and is inconsistent with the colours expected of a surviving donor star. Thus, the limiting magnitude of the HST observation (F555W > 28 mag) rules out all He star companion models and the most luminous main-sequence companion model currently in the literature. The remaining main-sequence companion models could be tested with a 10 orbit HST exposure in the F606W filter.

Do, Aaron↗

Modeling for pRad Explosive Experiments

The hydrodynamics code CTH was used to simulate high-explosive detonations in various geometries to design a planar, high-pressure drive with Taylor wave release for experiments in the proton radiography containment vessel. The most successful geometry is based on a P25 lens with some of its outer diameter and output end removed. Shock planarity was significantly improved relative to an unmodified P25 lens for this case, where its full diameter was not needed, and the calculated timing was adequate to trigger the experiments successfully.

36 MATERIALS SCIENCE↗

An Improved Method for Coupling Hydrodynamics with Astrophysical Reaction Networks

Reacting astrophysical flows can be challenging to model, because of the difficulty in accurately coupling hydrodynamics and reactions. This can be particularly acute during explosive burning or at high temperatures where nuclear statistical equilibrium is established. We develop a new approach, based on the ideas of spectral deferred corrections (SDC) coupling of explicit hydrodynamics and stiff reaction sources as an alternative to operator splitting, that is simpler than the more comprehensive SDC approach we demonstrated previously. We apply the new method to a double-detonation problem with a moderately sized astrophysical nuclear reaction network and explore the time step size and reaction network tolerances, to show that the simplified-SDC approach provides improved coupling with decreased computational expense compared to traditional Strang operator splitting. This is all done in the framework of the Castro hydrodynamics code, and all algorithm implementations are freely available.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Experiments on the Motion of Drops on a Horizontal Solid Surface due to a Wettability Gradient

Results from experiments performed on the motion of drops of tetraethylene glycol in a wettability gradient present on a silicon surface are reported and compared with predictions from a recently developed theoretical model. The gradient in wettability was formed by exposing strips cut from a silicon wafer to decyltrichlorosiland vapors. Video images of the drops captured during the experiments were subsequently analyzed for drop size and velocity as functions of position along the gradient. In separate experiments on the same strips, the static contact angle formed by small drops was measured and used to obtain the local wettability gradient to which a drop is subjected. The velocity of the drops was found to be a strong function of position along the gradient. A quasi-steady theoretical model that balances the local hydrodynamic resistance with the local driving force generally describes the observations; possible reasons for the remaining discrepancies are discussed. It is shown that a model in which the driving force is reduced to accomodate the hysteresis effect inferred from the data is able to remove most of the discrepancy between the observed and predicted velocities.

Moumen, Nadjoua↗

Quantitative assessment of eddy viscosity rans models for turbulent mixed convection in a differentially heated plane channel

Turbulent mixed convection between two vertical, infinite parallel plates at different temperatures is studied using various two-equation turbulence models. The numerical simulations are performed at a turbulent Reynolds number of Re τ = 150 and a Grashof number of Gr = 9.6 × 10 5 . Comparisons are made against the highly trusted DNS results. Consistent with the DNS approach, the current simulations are performed using constant properties and the Boussinesq approximation to predict the influence of buoyancy. Previous studies have provided assessments of two-equation turbulence models for various scenarios, but often rely on a qualitative “eye” test in order to determine the most appropriate model to predict a given flow. This study aims to provide a new form of quantitative assessment that accounts for both the physics captured by the turbulence model as well as the magnitude of the system response quantities (SRQ) using a modified symmetric mean absolute percent error (SMAPE) method. This method is designed to be approachable to researchers at any level and can be applied to system response quantities from multiple research fields. Uncertainty quantification is also performed to determine the discretization error for each turbulence model. Recommendations are made as to which turbulence models best capture the physics – hydrodynamically and thermally – using both local and global validation metrics. Lastly, a sensitivity analysis is performed on the damping functions used in the most accurate models. This underpins the potential of model developments and adjustments most worth pursuing for buoyant flows. Finally, this framework provides a more physics-based comparative analysis of the selected turbulence models.

42 ENGINEERING↗

Geothermal reservoir engineering research

The Stanford University research program on the study of stimulation and reservoir engineering of geothermal resources commenced as an interdisciplinary program in September, 1972. The broad objectives of this program have been: (1) the development of experimental and computational data to evaluate the optimum performance of fracture-stimulated geothermal reservoirs; (2) the development of a geothermal reservoir model to evaluate important thermophysical, hydrodynamic, and chemical parameters based on fluid-energy-volume balances as part of standard reservoir engineering practice; and (3) the construction of a laboratory model of an explosion-produced chimney to obtain experimental data on the processes of in-place boiling, moving flash fronts, and two-phase flow in porous and fractured hydrothermal reservoirs.

Ramey, H. J., Jr.↗

Modeling the Influence of Moulin Shape on Subglacial Hydrology

Subglacial models represent moulins as cylinders or cones, but field observations suggest that the upper part of moulins in the Greenland Ice Sheet has more complex shapes. These more complex shapes should cause englacial water storage within moulins to vary as a function of depth, a relationship not currently accounted for in models. Here, we use a coupled englacial-subglacial channel model to explore how moulin shape affects depth-dependent moulin water storage and water pressure dynamics within a subglacial channel. We simulate seven different moulin shapes across a range of moulin sizes. We find that the englacial storage capacity at the water level is the main control over the daily water level oscillation range and that depth-varying changes in englacial water storage control the temporal shape of this oscillation. Further, the cross-sectional area of the moulin within the daily oscillation range, but not above or below this range, controls pressures within the connected subglacial channel. Specifically, large cross-sectional areas can dampen daily to weekly oscillations that occur in the surface meltwater supply. Our findings suggest that further knowledge of the shape of moulins around the equilibrium water level would improve englacial storage parameterization in subglacial hydrological models and aid predictions of hydrodynamic coupling.

Moulin↗