SEARCH · Search NASA

Results for “hydrodynamic efficiency”

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 343 records · Page 19

All sources NASA NTRS DOE OSTI

DOE OSTI2010

Total angular momentum and reaction plane in collisions of high-energy particles

The reaction plane is an efficient tool for studying nucleus-nucleus collisions and other particle interactions. In order to determine it in those cases, use is made of model (hydrodynamic) concepts of the relationship between the impact parameter and collective flows of particles. It is shown that there is a direct possibility for reconstructing, in a model-independent way, the orientation of the plane in which particle collisions in reactions occur. The idea of the method is based on the relationship between experimentally measured quantities and the total angular momentum of the system of particles involved in the collision process. The concept of the reaction plane is defined for the case of multiparticle-production processes; it is shown that a normal to it is aligned with the total angular momentum. This circumstance ultimately leads to the interplay of particle azimuthal angles. This effect is interpreted in the literature as a manifestation of a collective production of particle flows in nuclear collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2013

Influence of the chlorine concentration on the radiation efficiency of a XeCl exciplex lamp

The influence of the chlorine concentration on the radiation efficiency of coaxial exciplex lamps (excilamps) excited by a dielectric barrier discharge (DBD) in binary Xe-Cl{sub 2} mixtures at pressures of 240–250 Torr is investigated experimentally and theoretically. The experiments were carried out at Cl{sub 2} concentrations in the range of 0.01–1%. The DBD characteristics were calculated in the framework of a one-dimensional hydrodynamic model at Cl{sub 2} concentrations in the range of 0.1–5%. It is found that the radiation intensities of the emission bands of Xe*{sub 2}(172 nm) and XeCl* (308 nm) are comparable when the chlorine concentration in the mixture is in the range of 0.01–0.1%. In this case, in the mixture, the radiation intensity of the Xe*{sub 2} molecule rapidly decreases with increasing Cl{sub 2} concentration and, at a chlorine concentration of ≥0.2%, the radiation of the B → X band of XeCl* molecules with a peak at 308 nm dominates in the discharge radiation. The radiation efficiency of this band reaches its maximum value at chlorine concentrations in the range of 0.4–0.5%. The calculated efficiencies of DBD radiation exceed those obtained experimentally. This is due to limitations of the one-dimensional model, which assumes the discharge to be uniform in the transverse direction, whereas the actual excilamp discharge is highly inhomogeneous. The influence of the chlorine concentration on the properties of the DBD plasma in binary Xe-Cl{sub 2} mixtures is studied numerically. It is shown that an increase in the Cl{sub 2} concentration in the mixture leads to the attachment of electrons to chlorine atoms and a decrease in the electron density and discharge conductivity. As a result, the electric field and the voltage drop across the discharge gap increase, which, in turn, leads to an increase in the average electron energy and the probability of dissociation of Cl{sub 2} molecules and ionization of Xe atoms and Cl{sub 2} molecules. The total energy deposited in the discharge rises with increasing chlorine concentration due to an increase in the power spent on the heating of positive and negative ions. The power dissipated by electrons decreases with increasing chlorine concentration in the working mixture. Recommendations on the choice of the chlorine content in the mixture for reducing the intensity of VUV radiation of the second continuum of the Xe*{sub 2} excimer without a substantial decrease in the excilamp efficiency are formulated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

The distribution-based remapping of the nodal mass and momentum between arbitrary meshes for staggered arbitrary Lagrangian-Eulerian hydrodynamics

We present a new distribution-based remapping of the nodal mass and momentum between arbitrary source (Lagrangian) and arbitrary target (rezoned) meshes for indirect staggered arbitrary Lagrangian-Eulerian hydrodynamics. Here, the method is based on the following ideas: define cell-centered momentum and mass on the source mesh; conservatively remap those cell-centered quantities from source to target meshes; and use local constrained optimization for each cell of the target mesh to conservatively distribute cell mass and momentum to the nodes of the cell. This new method is efficient, conservative, accurate and bounds preserving.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2017

Chemical hydrodynamics of a downward microbubble flow for intensification of gas‐fed bioreactors

Bioreactors are of interest for value‐upgrading of stranded or waste industrial gases. Reactor intensification requires development of low cost bioreactors with fast gas–liquid mass transfer rate. Here we assess published reactor technology in comparison with a novel downward bubble flow created by a micro‐jet array. Compared to known technology, the advanced design achieves higher volumetric gas transfer efficiency (k L a per power density) and can operate at higher k L a. We measure the effect of four reactor heights (height‐to‐diameter ratios of 12, 9, 6, and 3) on the gas transfer coefficient k L , total interfacial area a, liquid residence time distribution, energy consumption, and turbulent hydrodynamics. Leading models for predicting k L and a are appraised with experimental data. The results show k L is governed by “entrance effects” due to Higbie penetration dominate at short distances below the micro‐jet array, while turbulence dominates at intermediate distances, and finally terminal rise velocity dominates at large distances. © 2017 American Institute of Chemical Engineers AIChE J , 64: 1399–1411, 2018

Ansari, Manizheh↗

DOE OSTI2012

Hydrodynamic Modeling Analysis of Tidal Wetland Restoration in Snohomish River, Washington

To re-establish the intertidal wetlands with full tidal interaction and improve salmonid rearing habitat in the Lower Snohomish River estuary, a diked wetland along Union Slough of the Snohomish River was restored by breaching the existing dike and constructing bridges across the breaches. However, post-restoration monitoring indicated that the restored project site could not drain as efficiently as desired. To improve the drainage conditions at the restoration site during low tides, a modeling study was conducted to evaluate additional restoration scenarios and to provide recommendations for finish-grade ground elevations to achieve the desired drainage. To accurately simulate the drainage of the project site, an unstructured-grid hydrodynamic model with fine-grid resolution down to a few meters was used in this study. The model was first validated with observed water level data collected in the project site and then applied to assess the feasibility of different proposed restoration scenarios. A spatial varying bottom roughness option in the model is also implemented to better represent the high roughness due to the presence of dense vegetation in the project site. The methodology, error statistics of model validation and uncertainty of the modeling analysis are presented and discussed.

Yang, Zhaoqing↗

DOE OSTI2017

Numerical investigation of a joint approach to thermal energy storage and compressed air energy storage in aquifers

Different from conventional compressed air energy storage (CAES) systems, the advanced adiabatic compressed air energy storage (AA-CAES) system can store the compression heat which can be used to reheat air during the electricity generation stage. Thus, AA-CAES system can achieve a higher energy storage efficiency. Similar to the AA-CAES system, a compressed air energy storage in aquifers (CAESA) system, which is integrated with an aquifer thermal energy storage (ATES) could possibly achieve the same objective. In order to investigate the impact of ATES on the performance of CAESA, different injection air temperature schemes are designed and analyzed by using numerical simulations. Key parameters relative to energy recovery efficiencies of the different injection schemes, such as pressure distribution and temperature variation within the aquifers as well as energy flow rate in the injection well, are also investigated in this study. The simulations show that, although different injection schemes have a similar overall energy recovery efficiency (~97%) as well as a thermal energy recovery efficiency (~79.2%), the higher injection air temperature has a higher energy storage capability. Our results show the total energy storage for the injection air temperature at 80 °C is about 10% greater than the base model scheme at 40 °C. Sensitivity analysis reveal that permeability of the reservoir boundary could have significant impact on the system performance. However, other hydrodynamic and thermodynamic properties, such as the storage reservoir permeability, thermal conductivity, rock grain specific heat and rock grain density, have little impact on storage capability and the energy flow rate. Our study suggests that the combination of ATES and CAESA can help keep the high efficiency of energy storage so as to make CAESA system more efficiency.

25 ENERGY STORAGE↗

NASA NTRSDate not supplied

Two-Phase Flow in Filters and Orifices: In preparation for Measurements on ISS Packed Bed Reactor Experiment-Water Recovery Series (PBRE-WRS)

Understanding the hydrodynamics of adiabatic two-phase flows in packed bed reactors offers numerous benefits. It enables the optimization of chemical reactions rates and products, crucial for pharmaceuticals and energy production industries. Additionally, this understanding aids in designing more efficient and compact reactors, reducing space and resource requirements. Insights gained from studying such flows in microgravity contribute to advancements of space technologies and the enhancement of our capabilities for undertaking future long duration safe, and sustainable space exploration missions. The Packed Bed Reactor Experiment-Water Recovery Series (PBRE-WRS) is a flight experiment planned for operation on the Microgravity Science Glove box (MSG) facility of the International Space Station (ISS). The objective of the experiment is to validate hydrodynamic models at a range of gas and liquid flow rates. The microgravity environment allows for measurement and observation of aspects of fluid dynamics that are unique when compared to observations made in full or partial gravity environments. The experiment consists of testing two-phase flow hydrodynamics in three different filters, four (4) orifices and one check valve test articles. The experiment utilizes the hardware of the previously flown PBRE and PBRE-2 with packed beds of different packings materials and sizes. The fluid system of the PBRE-WRS consists of a nitrogen gas delivery and water delivery subsystems. The gas delivery subsystem can accommodate two ranges of gas flow rates, just as for the water delivery system. The nitrogen gas loop is open whereas the water loop is closed. Gas is separated from the water using a gas-liquid phase separator which is located in the water delivery module. Preliminary results from ground testing show the pressure drop increasing at different rates with the liquid and gas flow rates in the Brine filter. In this work, a detailed system description is presented along with a summary of results from ground performance testing of different test sections in an end-to-end preflight testing campaign.

Packed Bed Reactor Experiment↗

DOE OSTI2014

Fluctuating hydrodynamics of multispecies nonreactive mixtures

In this study we discuss the formulation of the fluctuating Navier-Stokes equations for multispecies, nonreactive fluids. In particular, we establish a form suitable for numerical solution of the resulting stochastic partial differential equations. An accurate and efficient numerical scheme, based on our previous methods for single species and binary mixtures, is presented and tested at equilibrium as well as for a variety of nonequilibrium problems. These include the study of giant nonequilibrium concentration fluctuations in a ternary mixture in the presence of a diffusion barrier, the triggering of a Rayleigh-Taylor instability by diffusion in a four-species mixture, as well as reverse diffusion in a ternary mixture. Finally, good agreement with theory and experiment demonstrates that the formulation is robust and can serve as a useful tool in the study of thermal fluctuations for multispecies fluids.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2026

Enabling microbial electrolysis cell scale-up via electrochemistry-, hydrodynamic-, and microbial ecology-informed framework

Microbial electrolysis cells (MECs) can produce green hydrogen while removing organic contaminants from liquid waste streams by leveraging the metabolic activity of electroactive microorganisms. Despite their potential in a sustainable, circular economy, large-scale MECs that can treat relevant volumes of wastewater have failed to deliver performance proportional to their lab-scale counterparts. The reason behind this lower performance at scale remains unclear. Here, in this study, we developed a combined electrochemistry-, hydrodynamic-, and microbial ecology-informed framework to analyze and optimize MEC performance during scale-up, enabling accurate quantification of major limitations and the identification of strategies to overcome them, ultimately facilitating equivalent performance at scale. Applying this framework to the scale-up of a zero-gap MEC from 9 cm 2 electrode area to 100 cm 2 electrode area, resulted in similar maximum current densities in a 100 cm 2 MEC (21.7 ± 1.1 A/m 2 ) compared to a 9 cm 2 system (25.1 ± 2.7 A/m 2 ), as well as equivalent hydrogen production rates of 69.3 L/L-d (100 cm 2 ) and 67.7 ± 2.4 L/L-d (9 cm 2 ). COMSOL flow dynamics simulations were used to scale up the reactor configuration without negatively affecting electrolyte velocity and distribution in the cell, minimizing the increase in internal resistances during scale-up (11.7 ± 0.5 mΩm 2 at 9 cm 2 ; 19.7 ± 1.3 mΩm 2 at 100 cm 2 ). Microbial community structures were assessed at both scales using high-throughput sequencing, highlighting the differences of populations across electrode dimensions and operational parameters. The framework presented here accelerates the development of effective strategies toward the scale-up of MECs by furthering the understanding of how electrochemical, hydrodynamic, and microbial ecology parameters change as the reactor dimension is increased. Ultimately, this approach contributes to advancing electrochemical biotechnology toward practical deployment in energy-efficient wastewater treatment systems.

Flow path↗

DOE OSTI2017

Three-dimensional hydrodynamic simulations of OMEGA implosions

Here, the effects of large-scale (with Legendre modes ≲10) asymmetries in OMEGA direct-drive implosions caused by laser illumination nonuniformities (beam-power imbalance and beam mispointing and mistiming), target offset, and variation in target-layer thickness were investigated using the low-noise, three-dimensional Eulerian hydrodynamic code ASTER. Simulations indicate that these asymmetries can significantly degrade the implosion performance. The most important sources of the asymmetries are the target offsets (~10 to 20 μm), beam-power imbalance (σ rms ~ 10%), and variations (~5%) in target-layer thickness. Large-scale asymmetries distort implosion cores, resulting in a reduced hot-spot confinement and an increased residual kinetic energy of implosion targets. The ion temperature inferred from the width of simulated neutron spectra is influenced by bulk fuel motion in the distorted hot spot and can result in up to an ~1 -keV increase in apparent temperature. Similar temperature variations along different lines of sight are observed. Demonstrating hydrodynamic equivalence to ignition designs on OMEGA requires a reduction in large-scale target and laser-imposed nonuniformities, minimizing target offset, and employing highly efficient mid-adiabat (α = 4) implosion designs, which mitigate cross-beam energy transfer and suppress short-wavelength Rayleigh–Taylor growth.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2013

ON THE INITIAL MASS FUNCTION OF LOW-METALLICITY STARS: THE IMPORTANCE OF DUST COOLING

The first stars to form in the universe are believed to have distribution of masses biased toward massive stars. This contrasts with the present-day initial mass function, which has a predominance of stars with masses lower than 1 M{sub Sun }. Therefore, the mode of star formation must have changed as the universe evolved. Such a transition is attributed to a more efficient cooling provided by increasing metallicity. Especially dust cooling can overcome the compressional heating, which lowers the gas temperature thus increasing its instability to fragmentation. The purpose of this paper is to verify if dust cooling can efficiently cool the gas, and enhance the fragmentation of gas clouds at the early stages of the universe. To confirm that, we calculate a set of hydrodynamic simulations that include sink particles, which represent contracting protostars. The thermal evolution of the gas during the collapse is followed by making use of a primordial chemical network and also a recipe for dust cooling. We model four clouds with different amounts of metals (10{sup -4}, 10{sup -5}, 10-6 Z{sub Sun }, and 0), and analyze how this property affect the fragmentation of star-forming clouds. We find evidence for fragmentation in all four cases, and hence conclude that there is no critical metallicity below which fragmentation is impossible. Nevertheless, there is a clear change in the behavior of the clouds at Z {approx}< 10{sup -5} Z{sub Sun }, caused by the fact that at this metallicity, fragmentation takes longer to occur than accretion, leading to a flat mass function at lower metallicities.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2021

Partially averaged Navier-Stokes closure modeling for variable-density turbulent flow

We report this work extends the framework of the partially averaged Navier-Stokes (PANS) equations to variable-density flow, i.e., multimaterial and/or compressible mixing problems with density variations and production of turbulence kinetic energy by both shear and buoyancy mechanisms. The proposed methodology is utilized to derive the PANS BHR-LEVM closure. This includes a priori testing to analyze and develop guidelines toward the efficient selection of the parameters controlling the physical resolution and, consequently, the range of resolved scales of PANS. Two archetypal test-cases involving transient turbulence, hydrodynamic instabilities, and coherent structures are used to illustrate the accuracy and potential of the method: the Taylor-Green vortex at Reynolds number Re = 3000, and the Rayleigh-Taylor flow at Atwood number 0.5 and (Re) max ≈ 500. These representative problems, for which turbulence is generated by shear and buoyancy processes, constitute the initial validation space of the new model, and their results are comprehensively discussed in two subsequent studies. The computations indicate that PANS can accurately predict the selected flow problems, resolving only a fraction of the scales of large-eddy simulation and direct numerical simulation strategies. The results also reiterate that the physical resolution of the PANS model must guarantee that the key instabilities and coherent structures of the flow are resolved. The remaining scales can be modeled through an adequate turbulence scale-dependent closure.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2015

First high-convergence cryogenic implosion in a near-vacuum hohlraum

Recent experiments on the National Ignition Facility [M. J. Edwards et al., Phys. Plasmas 20, 070501 (2013)] demonstrate that utilizing a near-vacuum hohlraum (low pressure gas-filled) is a viable option for high convergence cryogenic deuterium-tritium (DT) layered capsule implosions. This is made possible by using a dense ablator (high-density carbon), which shortens the drive duration needed to achieve high convergence: a measured 40% higher hohlraum efficiency than typical gas-filled hohlraums, which requires less laser energy going into the hohlraum, and an observed better symmetry control than anticipated by standard hydrodynamics simulations. The first series of near-vacuum hohlraum experiments culminated in a 6.8 ns, 1.2 MJ laser pulse driving a 2-shock, high adiabat (α ~ 3.5) cryogenic DT layered high density carbon capsule. This resulted in one of the best performances so far on the NIF relative to laser energy, with a measured primary neutron yield of 1.8 X 10¹⁵ neutrons, with 20% calculated alpha heating at convergence ~27X.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2023

Recent Hydrodynamic Modeling Enhancements in OpenFAST

To improve the hydrodynamic modeling capabilities of the offshore wind design and modeling tool OpenFAST, a restructuring of the HydroDyn hydrodynamics module was undertaken with several new features implemented. The generation of the wave field is now separated from HydroDyn into a new module called SeaState. Unlike previous versions of HydroDyn, which precomputed the wave kinematics at the undisplaced positions of the hydrodynamic nodes, the new SeaState module computes the complete time history of the wave field at the vertices of a user-defined wave grid. During the simulation, the wave kinematics at any point within the grid can be efficiently interpolated, allowing the strip-theory wave loads to be evaluated based on the wave kinematics at the displaced structure position. A phase correction was also implemented for the potential-flow wave-exciting loads to account for large structure displacements. Several wave-stretching methods were implemented for the strip-theory solution along with a wave-load redistribution method that ensures smooth variation of the nodal loads as the nodes of a discrete hydrodynamic mesh enter and exit the water. The load smoothing is included to avoid the excitation of unphysical high-frequency structural vibration during hydroelastic simulations. The MacCamy-Fuchs diffraction correction for the strip-theory solution and the capability to insert a constrained NewWave into extreme stochastic sea states with directional spreading were also added. The present article documents the formulations of the new improvements to HydroDyn with example applications and numerical results.

constrained wave↗

DOE OSTI2019

Collaboration with NSTX in Calculations of Radiofrequency and Neutral Beam Heating and Current Drive

The overall objective of this project was to advance the understanding of heating and current drive necessary for toroidal magnetic confinement fusion energy reactors, and to provide accurate models for these sources to aid in understanding magneto-hydrodynamic, transport and boundary physics issues. Increased understanding is provided by computer simulation using the underlying physics laws, and will lead to greater reliability and increased efficiency of a power producing fusion reactor. The proposed work entailed applications of existing CompX Fokker-Planck, ray tracing, full-wave microwave, and particle diffusion calculator codes to the National Spherical Torus Experiment at Princeton Plasma Physics Laboratory. Special features of NSTX include that large ion orbits are emphasized in its geometry and thus offer a major opportunity within fusion research for focused exploration of their ubiquitous effects, including comparison with simulation codes. The CompX CQL3D Fokker-Planck code was be further adapted to the finite-orbit-width physics. Comparison with experiment ensured that the computer codes have sufficient accuracy to accurately model the experiments. Nonthermal distributions of energetic ions are produced in experiments-- and modeled with the codes-- exhibiting the signatures of auxiliary heating, transport, and radiation processes within the plasma. Using physics-based modeling, in conjunction with experimental data, nonthermal distributions can provide keys to validation and quantification of auxiliary heating physics. Nonthermal distributions also hold particular information on radial transport processes, not obscured by the smoothing Maxwellization effects of Coulomb collisions. Generally, radial transport is not well understood but is a key to efficient fusion energy production. By comparison with Fokker-Planck modeling, the velocity dependence of radial transport coefficients becomes experimentally accessible, particularly in cases of well-modeled sources such as neutral beams, fast wave power, or the toroidal electric field. It is clear from modeling of fast-wave interactions with ions that finite-orbit-width effects play a significant role. For example, a different set of ions is computed to intersect a given cyclotron resonance depending on whether or not radial orbit drifts are taken into account. Also, radial transport, associated bootstrap current, and alpha channeling depend on radial drifts. Ion cyclotron heating above the first harmonic interacts preferentially with higher velocity ions, giving large radial drifts. The CQL3D Fokker-Planck code, maintained by the proposers, is a physics based model used extensively within the US fusion community to compare with and interpret experimental observations. At the beginning of this grant, bounce-averages of the rf quasilinear and other operators have been performed neglecting the radial drifts which give finite orbit widths to the particle orbits. The proposed work included a major upgrade to CQL3D to include finite-orbit-width effects be accounted for by performing the bounce-averages of rf and neutral beam fast ions along guiding center orbits. This CQL3D enhancement provided microwave and neutral beam contributions to neoclassical radial transport coefficients for nonthermal distributions, and substantially increased the comprehensiveness of the ion model. Additional work validated Fokker-Planck and ray tracing modeling of neutral beam and high harmonic fast wave absorption and current drive against the NSTX experiment, as had been very successful for electron cyclotron current drive in the DIII-D tokamak at General Atomics and lower hybrid waves on C-Mod at MIT. The validation work was in conjunction with staff and graduate students at Princeton Plasma Physics Laboratory. This DOE grant supported fusion energy research, a potential long-term solution to the world's energy needs. The research grant under consideration, DE-SC0006614, provided research support at the 0.3 Full-Time-Equivalent level for application of computer models to aid in understanding and projecting efficacy of heating and current drive sources in the National Spherical Torus Experiment, a tokamak variant, at the Princeton Plasma Physics Laboratory. The NSTX experiment explores the physics of very tight aspect ratio, almost spherical tokamaks, aiming at producing steady-state fusion plasmas. The current drive is an integral part of the steady-state concept, maintaining the magnetic geometry in the steady-state tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2016

Angular momentum loss in the envelope-disk transition region of the HH 111 protostellar system: evidence for magnetic braking?

HH 111 is a Class I protostellar system at a distance of ∼400 pc, with the central source VLA 1 associated with a rotating disk deeply embedded in a flattened envelope. Here we present the observations of this system at ∼0.″6 (240 au) resolution in C{sup 18}O (J = 2 — 1) and a 230 GHz continuum obtained with the Atacama Large Millimeter/Submillimeter Array, and in SO (N{sub J} = 5{sub 6}−4{sub 5}) obtained with the Submillimeter Array. The observations show for the first time how a Keplerian rotating disk can be formed inside a flattened envelope. The flattened envelope is detected in C{sup 18}O, extending out to ≳2400 au from the VLA 1 source. It has a differential rotation, with the outer part (≳2000 au) better described by a rotation that has constant specific angular momentum, and the innermost part (≲160 au) by a Keplerian rotation. The rotationally supported disk is therefore relatively compact in this system, which is consistent with the dust continuum observations. Most interestingly, if the flow is in steady state, there is a substantial drop in specific angular momentum in the envelope–disk transition region from 2000 to 160 au, by a factor of ∼3. Such a decrease is not expected outside a disk formed from simple hydrodynamic core collapse, but can happen naturally if the core is significantly magnetized, because magnetic fields can be trapped in the transition region outside the disk by the ram pressure of the protostellar accretion flow, which can lead to efficient magnetic braking. In addition, SO shock emission is detected around the outer radius of the disk and could trace an accretion shock around the disk.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Recent Hydrodynamic Modeling Enhancements in OpenFAST: Preprint

constrained wave↗

DOE OSTI2010

A Novel Method for Modeling Neumann and Robin Boundary Conditions in Smoothed Particle Hydrodynamics

In this paper we present an improved method for handling Neumann or Robin boundary conditions in smoothed particle hydrodynamics. The Neumann and Robin boundary conditions are common to many physical problems (such as heat/mass transfer), and can prove challenging to model in volumetric modeling techniques such as smoothed particle hydrodynamics (SPH). A new SPH method for diffusion type equations subject to Neumann or Robin boundary conditions is proposed. The new method is based on the continuum surface force model [1] and allows an efficient implementation of the Neumann and Robin boundary conditions in the SPH method for geometrically complex boundaries. The paper discusses the details of the method and the criteria needed to apply the model. The model is used to simulate diffusion and surface reactions and its accuracy is demonstrated through test cases for boundary conditions describing different surface reactions.

97 MATHEMATICS AND COMPUTING↗