Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “FLUID MECHANICS”

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,027 records · Page 57

Computational Investigation of Melt Pool Process Dynamics and Pore Formation in Laser Powder Bed Fusion

In the laser powder bed fusion additive manufacturing process, the presence of porosity may result in cracks and significantly affects the part performance. A comprehensive understanding of the melt pool process dynamics and porosity evolution can help to improve build quality. In this study, a novel multi-physics computational fluid dynamics (CFD) model has been applied to investigate the fluid dynamics in melt pools and resultant pore defects. To accurately capture the melting and solidification process, major process physics, such as the surface tension, evaporation as well as laser multi-reflection, have been considered in the model. A discrete element method is utilized to model the generation of powder spreading upon build plate by additional numerical simulations. Multiple single track experiments have been performed to obtain the melt pool shape and cross-sectional dimension information. The predicted melt pool dimensions were found to have a reasonable agreement with experimental measurements, e.g., the errors are in the range of 1.3 to 10.6% for melt pool width, while they are between 1.4 and 15.9% for melt depth. Pores are captured by both CFD simulation and x-ray computed tomography measurement for the case with a laser power of 350 W and laser speed of 100 mm/s. The formation of keyholes maybe related to the melt pool front wall angle, and it is found that the front wall angle increases with the increase in laser line energy density. In addition, a larger laser power or smaller scanning speed can help to generate keyhole-induced pores; they also contribute to produce larger sized pores.

36 MATERIALS SCIENCE↗

Systematic Ion Source Parameter Assessment by Automated Determination of the Distribution of Ion Acceptance (DIA) Using APLI

The determination of the spatially resolved ion signal with atmospheric pressure laser ionization (APLI), which was introduced as distribution of ion acceptance (DIA), serves as a valuable tool for the understanding of complex and highly dynamical conditions in modern atmospheric pressure (AP) ion sources. DIA provides information about fluid dynamics, ion transport, and ion transformation processes in such sources and is an ideal basis for the validation of numerical models of the dynamics in the ion source enclosure. We present a fully automated setup for DIA measurements, which enabled us to acquire a comprehensive dataset of over 700 individual DIA measurements in a commercial AP ion source (Bruker Multi Purpose Ion Source, MPIS). Ion source parameters as voltages, gas heater temperatures and gas flows, were varied, and the effect of those parameters on the DIA of a chemically inert analyte, pyrene, was systematically investigated. It is shown that the response of the DIA is nonlinear and that gas dynamics largely dominates the ion transport in the ion source. Particularly, the position of the heated nebulizer, which is used to introduce one of two gas flows and the analyte into the ion source chamber, had a profound effect on the DIA. This suggests that the gas dynamics in the source switches between different states. The now available comprehensive DIA dataset reveals such critical effects and will guide further numerical modeling efforts to understand the details of the dynamics of ions in the source chamber. .

74 ATOMIC AND MOLECULAR PHYSICS↗

Protorheology in practice: Avoiding misinterpretation

Protorheology is the paradigm that any observed flow or deformation is a chance to infer quantitative rheological properties. While this creates many opportunities for insight, there is significant risk of misunderstanding the physics involved, e.g. misinterpreting a liquid as a solid or mistaking viscous flow time as viscoelastic relaxation time. We describe these and other potential mistakes, use case studies to show how serious the problems can be, and contrast misinterpretations with correct approaches and interpretations. Some issues are especially important with materials involving colloidal particles and flows involving surface tension. Whether the reader is making inference from a tilted vial, time-lapse gravity-driven flow, a bounce test, die swell, or any other protorheology observation, the examples here serve as a guide for avoiding bad data in protorheology.

42 ENGINEERING↗

GAMBL – A dual-cooled fusion blanket using SiC-based structures

We describe a novel fusion reactor blanket concept called GAMBL – GA Modular BLanket. This design concept exploits the advantages of SiC-based structures while minimizing their limitations. Material microstructures are tailored for different functions throughout the first wall and blanket, including the use of graded W/SiC composites to enhance heat transfer capabilities and resist erosion at the plasma-facing surface. The design contains a decoupled first wall and breeding blanket, which we show can provide adequate tritium breeding. By decoupling the breeding part of the blanket, we can operate at very low PbLi breeder pressure and support gravity loads in simple, radiatively-cooled structural beams that do not contain coolant. First wall performance capabilities are enhanced by eliminating constraints imposed by the deep blanket. We have performed initial fluid, thermal, mechanical and neutronics analysis of this concept to show its capabilities to meet the design requirements. R&D needs are mostly related to materials development. Except for unique manufacturing needs, no new facilities are required beyond existing or planned facilities for the US base program on steel-based DCLL blankets.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Thermal loading effects on chalk hydromechanical behavior for nuclear waste disposal

Safe disposal of heat-generating nuclear waste depends on host rock stability under thermal, hydrological, and mechanical stresses. This study investigates the effect of thermal loading on mechanical behavior of the shallowly buried Ghareb formation chalk through triaxial and hydrostatic constant strain rate and creep tests at temperatures up to 100 ˚C and effective pressures up to 20.7 MPa. Experimental results show that thermal loading reduces the elastic moduli of chalk by 50–75%, and a transition occurs above 60 ˚C where creep rates increase rapidly. Water saturation nearly doubles the thermally induced strain compared to dry conditions and strongly decreases material rigidity. Thermal loading also leads to significant pore pressure increases under undrained conditions and leads to reductions in the apparent permeability during drained conditions. Laboratory experimental data were used to parameterize and develop a preliminary constitutive model for predicting future deformation during repository operations in the Ghareb. The strongly coupled effects – mechanical weakening, fluid pressure fluctuations, and permeability modification – demonstrate that elevated repository temperatures will have a pronounced effect on the near field Ghareb behavior during waste disposal operations. The findings indicate that the coupled interactions must be considered in predictive models and repository design to ensure long-term nuclear waste isolation and safety.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Lagrangian particle method for compressible fluid dynamics

A new Lagrangian particle method for solving Euler equations for compressible inviscid fluid or gas flows is proposed. Similar to smoothed particle hydrodynamics (SPH), the method represents fluid cells with Lagrangian particles and is suitable for the simulation of complex free surface / multi-phase flows. The main contributions of our method, which is different from SPH in all other aspects, are (a) significant improvement of approximation of differential operators based on a polynomial fit via weighted least squares approximation and the convergence of prescribed order, (b) a second-order particle-based algorithm that reduces to the first-order upwind method at local extremal points, providing accuracy and long term stability, and (c) more accurate resolution of entropy discontinuities and states at free inter-faces. While the method is consistent and convergent to a prescribed order, the conservation of momentum and energy is not exact and depends on the convergence order . The method is generalizable to coupled hyperbolic-elliptic systems. As a result, numerical verification tests demonstrating the convergence order are presented as well as examples of complex multiphase flows.

97 MATHEMATICS AND COMPUTING↗

Physics-informed neural networks: A deep learning framework for solving forward and inverse problems involving nonlinear partial differential equations

Hejre, we introduce physics-informed neural networks – neural networks that are trained to solve supervised learning tasks while respecting any given laws of physics described by general nonlinear partial differential equations. In this work, we present our developments in the context of solving two main classes of problems: data-driven solution and data-driven discovery of partial differential equations. Depending on the nature and arrangement of the available data, we devise two distinct types of algorithms, namely continuous time and discrete time models. The first type of models forms a new family of data-efficient spatio-temporal function approximators, while the latter type allows the use of arbitrarily accurate implicit Runge–Kutta time stepping schemes with unlimited number of stages. Moreover, the effectiveness of the proposed framework is demonstrated through a collection of classical problems in fluids, quantum mechanics, reaction–diffusion systems, and the propagation of nonlinear shallow-water waves.

97 MATHEMATICS AND COMPUTING↗

Meshfree methods on manifolds for hydrodynamic flows on curved surfaces: A Generalized Moving Least-Squares (GMLS) approach

We utilize generalized moving least squares (GMLS) to develop meshfree techniques for discretizing hydrodynamic flow problems on manifolds. We use exterior calculus to formulate incompressible hydrodynamic equations in the Stokesian regime and handle the divergence-free constraints via a generalized vector potential. This provides less coordinate-centric descriptions and enables the development of efficient numerical methods and splitting schemes for the fourth-order governing equations in terms of a system of second-order elliptic operators. Using a Hodge decomposition, we develop methods for manifolds having spherical topology. We show the methods exhibit high-order convergence rates for solving hydrodynamic flows on curved surfaces. The methods also provide general high-order approximations for the metric, curvature, and other geometric quantities of the manifold and associated exterior calculus operators. The approaches also can be utilized to develop high-order solvers for other scalar-valued and vector-valued problems on manifolds.

97 MATHEMATICS AND COMPUTING↗

Shale gas transport model in 3D fractal porous media with variable pore sizes

A model for gas transport in shale is proposed by accounting for three major fluid flow mechanisms in shale stratum, which is modeled as a 3D fractal media. The proposed apparent permeability of shale is an analytical expression that also accounts for heterogeneous pore sizes in shale stratum, and is verified using experimental datasets for methane and helium flow in shale. Results of sensitivity analysis indicate that surface diffusion of adsorbed gas plays an important role, specifically in smaller pores, while surface diffusion would be negligible in larger pores. Further, the proposed model shows that flow due to surface diffusion decreases moderately with the increase of isosteric adsorption heat, while it increases significantly with the increase of the maximum adsorption capacity. In conclusion, one of the key novelties of the proposed permeability model is that it accounts for pore size distribution to reveal novel insights on gas transport in shale that can be used to optimize gas production by operational controls (e.g. controlling reservoir pressure) as flow regimes change with time.

58 GEOSCIENCES↗

Unsteady aerodynamic loads on pitching aerofoils represented by Gaussian body force distributions

The actuator line model (ALM) is an approach commonly used to represent lifting and dragging devices like wings and blades in large-eddy simulations (LES). The crux of the ALM is the projection of the actuator point forces onto the LES grid by means of a Gaussian regularisation kernel. The minimum width of the kernel is constrained by the grid size; however, for most practical applications like LES of wind turbines, this value is an order of magnitude larger than the optimal value that maximises accuracy. This discrepancy motivated the development of corrections for the actuator line, which, however, neglect the effect of unsteady spanwise shed vorticity. In this work we develop a model for the impact of spanwise shed vorticity on the unsteady loading of an aerofoil modelled as a Gaussian body force distribution, where the model is applicable within the regime of unsteady attached flow. The model solution is derived both in the time and frequency domain and features an explicit dependence on the Gaussian kernel width. We verify the model with ALM-LES for both pitch steps and periodic pitching. The model solution is compared with Theodorsen theory and validated with both computational fluid dynamics using body fitted grids and experiment. It is concluded that the optimal kernel width for unsteady aerodynamics is approximately 40 % of the chord. The ALM is able to predict the magnitude of the unsteady loading up to a reduced frequency of 𝑘 ≈ 0.2.

17 WIND ENERGY↗

Episodic Earthquake Swarms in the Mineral Mountains, Utah Driven by the Roosevelt Hydrothermal System

Over 1,000 earthquakes (-2.0 < M < 2.0), identified using a matched-filter method, occurred in the Mineral Mountains, Utah between 2016 and 2019. The enhanced catalog is complete down to M -0.9 and contains roughly 15 times more events than originally cataloged. Earthquake relocation of ~800 earthquakes shows that activity is concentrated in a <2 km long E-W striking narrow zone, ~4 km east of the Roosevelt hydrothermal system. Two fault orientations, both N-S and E-W parallel to the Opal Mound and Mag Lee faults, respectively, are observed after computing composite focal mechanisms of highly similar earthquakes. Looking solely at the temporal distribution of the seismicity, we identify 15 periods of swarm-like activity, with two major clusters occurring in December 2016, recorded by three stations, and in October 2019 recorded by eight stations. The October 2019 swarm, the best recorded sequence in the area, provides evidence for the underlying triggering mechanism. We show that a complex mechanism of fluid diffusion and aseismic slip is responsible for the swarm evolution with migration velocities reaching 10 km/day. We hypothesize that these episodic swarms in the Mineral Mountains are primarily driven by migrating fluids that originate within the Roosevelt hydrothermal system.

58 GEOSCIENCES↗

Constraining maximum event magnitude during injection-triggered seismicity

Abstract Understanding mechanisms controlling fluid injection-triggered seismicity is key in defining strategies to ameliorate it. Recent triggered events (e.g. Pohang, Mw 5.5) have exceeded predictions of average energy release by a factor of >1000x, necessitating robust methodologies to both define critical antecedent conditions and to thereby constrain anticipated event size. We define maximum event magnitudes resulting from triggering as a function of pre-existing critical stresses and fluid injection volume. Fluid injection experiments on prestressed laboratory faults confirm these estimates of triggered moment magnitudes for varied boundary conditions and injection rates. In addition, observed ratios of shear slip to dilation rates on individual faults signal triggering and may serve as a measurable proxy for impending rupture. This new framework provides a robust method of constraining maximum event size for preloaded faults and unifies prior laboratory and field observations that span sixteen decades in injection volume and four decades in length scale.

58 GEOSCIENCES↗

Rethinking materials simulations: Blending direct numerical simulations with neural operators

Abstract Materials simulations based on direct numerical solvers are accurate but computationally expensive for predicting materials evolution across length- and time-scales, due to the complexity of the underlying evolution equations, the nature of multiscale spatiotemporal interactions, and the need to reach long-time integration. We develop a method that blends direct numerical solvers with neural operators to accelerate such simulations. This methodology is based on the integration of a community numerical solver with a U-Net neural operator, enhanced by a temporal-conditioning mechanism to enable accurate extrapolation and efficient time-to-solution predictions of the dynamics. We demonstrate the effectiveness of this hybrid framework on simulations of microstructure evolution via the phase-field method. Such simulations exhibit high spatial gradients and the co-evolution of different material phases with simultaneous slow and fast materials dynamics. We establish accurate extrapolation of the coupled solver with large speed-up compared to DNS depending on the hybrid strategy utilized. This methodology is generalizable to a broad range of materials simulations, from solid mechanics to fluid dynamics, geophysics, climate, and more.

36 MATERIALS SCIENCE↗

Effect of work of adhesion on deep bed filtration process

Collection of aerosol particles in the particular steps of the technology of their production, and purification of the air at the workplace and atmospheric environment, requires the efficient method of separation of particulate matter from the carrier gas. There are many papers published in last few years in which the deposition of particles on fibrous collectors is considered, Most of them assume that collisions between particle and collector surface is 100% effective. In this work we study the influence of particles and fiber properties on the deposition efficiency. For the purpose of this work the lattice-Boltzmann model describes fluid dynamics, while the solid particle motion is modeled by the Brownian dynamics. The interactions between particles and surface are modelled using energy balanced oscillatory model. The work of adhesion was estimated using Atomic Force Microscopy.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Coupling of electrostatic ion cyclotron and ion acoustic waves in the solar wind

The coupling of electrostatic ion cyclotron and ion acoustic waves is examined in three component magnetized plasma consisting of electrons, protons, and alpha particles. In the theoretical model relevant to solar wind plasma, electrons are assumed to be superthermal with kappa distribution and protons as well as alpha particles follow the fluid dynamical equations. A general linear dispersion relation is derived for such a plasma system which is analyzed both analytically and numerically. For parallel propagation, electrostatic ion cyclotron (proton and helium cyclotron) and ion acoustic (slow and fast) modes are decoupled. For oblique propagation, coupling between the cyclotron and acoustic modes occurs. Furthermore, when the angle of propagation is increased, the separation between acoustic and cyclotron modes increases which is an indication of weaker coupling at large angle of propagation. For perpendicular propagation, only cyclotron modes are observed. The effect of various parameters such as number density and temperature of alpha particles and superthermality on dispersion characteristics is examined in details. The coupling between various modes occurs for small values of wavenumber.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fluid dynamic design and experimental study of an aspirated temperature measurement platform used in climate observation

Due to the solar radiation effect, current air temperature sensors inside a thermometer screen or radiation shield may produce measurement errors that are 0.8 °C or higher. To improve the observation accuracy, an aspirated temperature measurement platform is designed. A computational fluid dynamics (CFD) method is implemented to analyze and calculate the radiation error of the aspirated temperature measurement platform under various environmental conditions. Then, a radiation error correction equation is obtained by fitting the CFD results using a genetic algorithm (GA) method. In order to verify the performance of the temperature sensor, the aspirated temperature measurement platform, temperature sensors with a naturally ventilated radiation shield, and a thermometer screen are characterized in the same environment to conduct the intercomparison. The average radiation errors of the sensors in the naturally ventilated radiation shield and the thermometer screen are 0.44 °C and 0.25 °C, respectively. In contrast, the radiation error of the aspirated temperature measurement platform is as low as 0.05 °C. This aspirated temperature sensor allows the radiation error to be reduced by approximately 88.6% compared to the naturally ventilated radiation shield, and allows the error to be reduced by a percentage of approximately 80% compared to the thermometer screen. The mean absolute error and root mean square error between the correction equation and experimental results are 0.032 °C and 0.036 °C, respectively, which demonstrates the accuracy of the CFD and GA methods proposed in this research.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Capsule physics comparison of National Ignition Facility implosion designs using plastic, high density carbon, and beryllium ablators

Indirect drive implosion experiments on the National Ignition Facility (NIF) [E. I. Moses et al., Phys. Plasmas 16, 041006 (2009)] have now tested three different ablator materials: glow discharge polymer plastic, high density carbon, and beryllium. How do these different ablators compare in current and proposed implosion experiments on NIF? What are the relative advantages and disadvantages of each? This paper compares these different ablator options in capsule-only simulations of current NIF experiments and potential future designs. The simulations compare the impact of the capsule fill tube, support tent, and interface surface roughness for each case, as well as all perturbations in combination. According to the simulations, each ablator is impacted by the various perturbation sources differently, and each material poses unique challenges in the pursuit of ignition on NIF.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Shock-driven hydrodynamic instability of a sinusoidally perturbed, high-Atwood number, oblique interface

A shock incident on an interface will cause any initial perturbations on that interface to grow. When the shock front is parallel to the interface, the perturbations grow due to the Richtmyer-Meshkov (RM) process. When there is some tilt between the shock front and the interface, shear flow will result across the postshock interface. Recent experiments on the OMEGA EP laser have studied the hydrodynamic instability growth which results from a supported shock interacting with a sinusoidally perturbed, oblique interface. The observed instability growth was dominated by Richtmyer-Meshkov at early times but became Kelvin-Helmholtz (KH)-like at late times. Previously, this instability growth was described using an analytic model for the deposition of baroclinic vorticity on the interface by a shock combined with a discrete vortex model. Here, we utilize the same baroclinic vorticity deposition model in conjunction with a desingularized, periodic Birkhoff-Rott equation to model instability evolution. The Birkhoff-Rott equation takes into account the vorticity distribution along the interface, whereas the discrete vortex model assumed that all vorticity over each wavelength of the perturbation is confined to a point. We compare the new model to xRAGE simulations and experiments. In conclusion, the model is found to overpredict both the instability growth and shear across the interface by about a factor of two, but correctly predicts that the growth is RM-like at early times and KH-like at late times.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗