Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fluids”

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,441 records · Page 80
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Computing the ankle-brachial index with parallel computational fluid dynamics

The ankle-brachial index (ABI), a ratio of arterial blood pressure in the ankles and upper arms, is used to diagnose and monitor circulatory conditions such as coarctation of the aorta and peripheral artery disease. Computational simulations of the ABI can potentially determine the parameters that produce an ABI indicative of ischemia or other abnormalities in blood flow. However, 0- and 1-D computational methods are limited in describing a 3-D patient-derived geometry. Thus, we present a massively parallel framework for computational fluid dynamics (CFD) simulations in the full arterial system. Additionally, using the lattice Boltzmann method to solve the Navier–Stokes equations, we employ highly parallelized and scalable methods to generate the simulation domain and efficiently distribute the computational load among processors. For the first time, we compute an ABI with 3-D CFD. In this proof-of-concept study, we investigate the dependence of ABI on the presence of stenoses, or narrowed regions of the arteries, by directly modifying the arterial geometry. As a result, our framework enables the computation a hemodynamic factor characterizing flow at the scale of the full arterial system, in a manner that is extensible to patient-specific imaging data and holds potential for treatment planning.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2014

Gradient Theory simulations of pure fluid interfaces using a generalized expression for influence parameters and a Helmholtz energy equation of state for fundamentally consistent two-phase calculations

The fidelity of Gradient Theory simulations depends on the accuracy of saturation properties and influence parameters, and require equations of state (EoS) which exhibit a fundamentally consistent behavior in the two-phase regime. Widely applied multi-parameter EoS, however, are generally invalid inside this region. Hence, they may not be fully suitable for application in concert with Gradient Theory despite their ability to accurately predict saturation properties. The commonly assumed temperature-dependence of pure component influence parameters usually restricts their validity to subcritical temperature regimes. This may distort predictions for general multi-component interfaces where temperatures often exceed the critical temperature of vapor phase components. Then, the calculation of influence parameters is not well defined. In this paper, one of the first studies is presented in which Gradient Theory is combined with a next-generation Helmholtz energy EoS which facilitates fundamentally consistent calculations over the entire two-phase regime. Illustrated on pentafluoroethane as an example, reference simulations using this method are performed. They demonstrate the significance of such high-accuracy and fundamentally consistent calculations for the computation of interfacial properties. These reference simulations are compared to corresponding results from cubic PR EoS, widely-applied in combination with Gradient Theory, and mBWR EoS. The analysis reveals that neither of those two methods succeeds to consistently capture the qualitative distribution of obtained key thermodynamic properties in Gradient Theory. Furthermore, a generalized expression of the pure component influence parameter is presented. This development is informed by its fundamental definition based on the direct correlation function of the homogeneous fluid and by presented high-fidelity simulations of interfacial density profiles. As a result, the new model preserves the accuracy of previous temperature-dependent expressions, remains well-defined at supercritical temperatures, and is fully suitable for calculations of general multi-component two-phase interfaces.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2013

Numerical schemes for dynamically orthogonal equations of stochastic fluid and ocean flows

The quantification of uncertainties is critical when systems are nonlinear and have uncertain terms in their governing equations or are constrained by limited knowledge of initial and boundary conditions. Such situations are common in multiscale, intermittent and non-homogeneous fluid and ocean flows. The dynamically orthogonal (DO) field equations provide an adaptive methodology to predict the probability density functions of such flows. The present work derives efficient computational schemes for the DO methodology applied to unsteady stochastic Navier-Stokes and Boussinesq equations, and illustrates and studies the numerical aspects of these schemes. Semi-implicit projection methods are developed for the mean and for the DO modes, and time-marching schemes of first to fourth order are used for the stochastic coefficients. Conservative second-order finite-volumes are employed in physical space with new advection schemes based on total variation diminishing methods. Other results include: (i) the definition of pseudo-stochastic pressures to obtain a number of pressure equations that is linear in the subspace size instead of quadratic; (ii) symmetric advection schemes for the stochastic velocities; (iii) the use of generalized inversion to deal with singular subspace covariances or deterministic modes; and (iv) schemes to maintain orthonormal modes at the numerical level. To verify our implementation and study the properties of our schemes and their variations, a set of stochastic flow benchmarks are defined including asymmetric Dirac and symmetric lock-exchange flows, lid-driven cavity flows, and flows past objects in a confined channel. Different Reynolds number and Grashof number regimes are employed to illustrate robustness. Optimal convergence under both time and space refinements is shown as well as the convergence of the probability density functions with the number of stochastic realizations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Fluid preconditioning for Newton–Krylov-based, fully implicit, electrostatic particle-in-cell simulations

A recent proof-of-principle study proposes an energy- and charge-conserving, nonlinearly implicit electrostatic particle-in-cell (PIC) algorithm in one dimension [9]. The algorithm in the reference employs an unpreconditioned Jacobian-free Newton–Krylov method, which ensures nonlinear convergence at every timestep (resolving the dynamical timescale of interest). Kinetic enslavement, which is one key component of the algorithm, not only enables fully implicit PIC as a practical approach, but also allows preconditioning the kinetic solver with a fluid approximation. This study proposes such a preconditioner, in which the linearized moment equations are closed with moments computed from particles. Effective acceleration of the linear GMRES solve is demonstrated, on both uniform and non-uniform meshes. The algorithm performance is largely insensitive to the electron–ion mass ratio. Numerical experiments are performed on a 1D multi-scale ion acoustic wave test problem.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

A hyperbolic-equation system approach for magnetized electron fluids in quasi-neutral plasmas

A new approach using a hyperbolic-equation system (HES) is proposed to solve for the electron fluids in quasi-neutral plasmas. The HES approach avoids treatments of cross-diffusion terms which cause numerical instabilities in conventional approaches using an elliptic equation (EE). A test calculation reveals that the HES approach can robustly solve problems of strong magnetic confinement by using an upwind method. The computation time of the HES approach is compared with that of the EE approach in terms of the size of the problem and the strength of magnetic confinement. The results indicate that the HES approach can be used to solve problems in a simple structured mesh without increasing computational time compared to the EE approach and that it features fast convergence in conditions of strong magnetic confinement.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

A flux-splitting method for hyperbolic-equation system of magnetized electron fluids in quasi-neutral plasmas

A flux-splitting method is proposed for the hyperbolic-equation system (HES) of magnetized electron fluids in quasi-neutral plasmas. The numerical fluxes are split into four categories, which are computed by using an upwind method which incorporates a flux-vector splitting (FVS) and advection upstream splitting method (AUSM). The method is applied to a test calculation condition of uniformly distributed and angled magnetic lines of force. All of the pseudo-time advancement terms converge monotonically and the conservation laws are strictly satisfied in the steady state. The calculation results are compared with those computed by using the elliptic–parabolic-equation system (EPES) approach using a magnetic-field-aligned mesh (MFAM). Both qualitative and quantitative comparisons yield good agreements of results, indicating that the HES approach with the flux-splitting method attains a high computational accuracy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

A stochastic asymptotic-preserving scheme for a kinetic-fluid model for disperse two-phase flows with uncertainty

In this paper we consider a kinetic-fluid model for disperse two-phase flows with uncertainty. We propose a stochastic asymptotic-preserving (s-AP) scheme in the generalized polynomial chaos stochastic Galerkin (gPC-sG) framework, which allows the efficient computation of the problem in both kinetic and hydrodynamic regimes. The s-AP property is proved by deriving the equilibrium of the gPC version of the Fokker–Planck operator. The coefficient matrices that arise in a Helmholtz equation and a Poisson equation, essential ingredients of the algorithms, are proved to be positive definite under reasonable and mild assumptions. The computation of the gPC version of a translation operator that arises in the inversion of the Fokker–Planck operator is accelerated by a spectrally accurate splitting method. Numerical examples illustrate the s-AP property and the efficiency of the gPC-sG method in various asymptotic regimes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

A hybrid, coupled approach for modeling charged fluids from the nano to the mesoscale

Here, we develop and demonstrate a new, hybrid simulation approach for charged fluids, which combines the accuracy of the nonlocal, classical density functional theory (cDFT) with the efficiency of the Poisson–Nernst–Planck (PNP) equations. The approach is motivated by the fact that the more accurate description of the physics in the cDFT model is required only near the charged surfaces, while away from these regions the PNP equations provide an acceptable representation of the ionic system. We formulate the hybrid approach in two stages. The first stage defines a coupled hybrid model in which the PNP and cDFT equations act independently on two overlapping domains, subject to suitable interface coupling conditions. At the second stage we apply the principles of the alternating Schwarz method to the hybrid model by using the interface conditions to define the appropriate boundary conditions and volume constraints exchanged between the PNP and the cDFT subdomains. Numerical examples with two representative examples of ionic systems demonstrate the numerical properties of the method and its potential to reduce the computational cost of a full cDFT calculation, while retaining the accuracy of the latter near the charged surfaces.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Assessment of diffuse-interface methods for compressible multiphase fluid flows and elastic-plastic deformation in solids

This work describes three diffuse-interface methods for the simulation of immiscible, compressible multiphase fluid flows and elastic-plastic deformation in solids. The first method is the localized-artificial-diffusivity approach of Cook, Subramaniam et al., and Adler and Lele, in which artificial diffusion terms are added to the individual phase mass fraction transport equations and are coupled with the other conservation equations. The second method is the gradient-form approach that is based on the quasi-conservative method of Shukla et al., in which the diffusion and sharpening terms (together called regularization terms) are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. The third approach is the divergence-form approach that is based on the fully conservative method of Jain et al., in which the regularization terms are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. In the present study, all three diffuse-interface methods are used in conjunction with a four-equation, multicomponent mixture model, in which pressure and temperature equilibria are assumed among the various phases. The primary objective of this work is to compare these three methods in terms of their ability to: maintain constant interface thickness throughout the simulation; conserve mass, momentum, and energy; and maintain accurate interface shape for long-time integration. The second objective of this work is to consistently extend these methods to model interfaces between solid materials with strength. To assess and compare the methods, they are used to simulate a wide variety of problems, including (1) advection of an air bubble in water, (2) shock interaction with a helium bubble in air, (3) shock interaction and the collapse of an air bubble in water, and (4) Richtmyer–Meshkov instability of a copper–aluminum interface. The current work focuses on comparing these methods in the limit of relatively coarse grid resolution, which illustrates the true performance of these methods. In conclusion, this is because it is rarely practical to use hundreds of grid points to resolve a single bubble or drop in large-scale simulations of engineering interest.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Multi-plane moment-of-fluid interface reconstruction in 3D

Moment-of-fluid (MOF) methods for interface reconstruction approximate the region occupied by material in each mesh element only through reference to its geometric moments. Here, we present a 3D MOF method that represents the material (POM) in each cell as the convex intersection of the cell and multiple half-spaces, each selected to minimize the least-squares error between computed moments of the approximated material and provided reference moments. This optimization problem is highly non-linear and non-convex, making the numerical result very sensitive to the initial guess. To create an effective initial guess in each cell, we construct an ellipsoid from 0th–2nd order reference moments such that its shape corresponds with that of the POM. Within this ellipsoid we inscribe a polyhedron, and initialize the minimization problem with the half-spaces defined by each of its faces. The inscribed polyhedron has minimally 4 faces, and using up to 3rd order moments permits optimization over up to 20 unknown values. We therefore define MOF methods that utilize 4, 5, or 6 half-spaces, correspondingly initialized with the faces of a single inscribed tetrahedron, triangular prism, or hexahedron. Stability of the non-linear optimization is further improved with a prepossessing step that normalizes the reference moments according to the axes of the reference ellipsoid. Using this approach, the non-linear least-squares solver reliably converges to a near-global minimum from a single initial guess. We demonstrate accuracy and robustness using single-cell and multi-cell examples over a wide spectrum of geometry. In particular, we demonstrate our ability to exactly reproduce several important and complex features defined by up to four half-spaces, such as corners, filaments, filament tips, and embedded material in the cell.

3D interface reconstruction↗
DOE OSTI2020

A thin-walled composite beam model for light-weighted structures interacting with fluids

A thin-walled beam model is proposed for structures of variable cross-section, which can be either open or closed and includes multicellular cross-sections with either isotropic or orthotropic materials. The proposed model does not require any priori definition of cross-sectional warping which instead results from the solution of the problem. To achieve that a special deformation pattern is superimposed on the bending deformation described by Euler–Bernoulli beam theory. All sectional properties are automatically incorporated in the analysis as a result of the usual variational formulation of the system of equations. The proposed model is specifically designed to simulate the dynamics of wind/hydrokinetic turbine blade with low computational cost, especially in fluid–structure interaction (FSI) simulation. A number of test cases have been carried out to validate the proposed structural model which show good agreement between the results obtained her e and the solutions available in literature. Finally, FSI simulation of a hydrokinetic blade under field condition is carried out to illustrate the capability of the current thin-walled beam model in practice.

42 ENGINEERING↗
DOE OSTI2017

2D Materials with Nanoconfined Fluids for Electrochemical Energy Storage

In the quest to develop energy storage with both high power and high energy densities, and while maintaining high volumetric capacity, recent results show that a variety of 2D and layered materials exhibit rapid kinetics of ion transport by the incorporation of nanoconfined fluids.

25 ENERGY STORAGE↗
DOE OSTI2019

Performance portability study for massively parallel computational fluid dynamics application on scalable heterogeneous architectures

Patient-specific hemodynamic simulations have the possibility to greatly improve both the diagnosis and treatment of a variety of vascular diseases. Portability will enable wider adoption of computational fluid dynamics (CFD) applications in the biomedical research community and targeting to platforms ideally suited to different vascular regions. Here, we introduce a case study in performance portability that assesses (1) the ease of porting an MPI application optimized for one specific architecture to new platforms using variants of hybrid MPI+X programming models; (2) performance portability seen when simulating blood flow in three different vascular regions on diverse heterogeneous architectures; (3) model-based performance prediction for future architectures; and (4) performance scaling of the hybrid MPI+X programming on parallel heterogeneous systems. We report on the lessons learned in porting HARVEY, a massively parallel CFD application, from traditional multicore CPUs to diverse heterogeneous architectures ranging from NVIDIA/AMD GPUs to Intel MICs and Altera FPGAs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Boron doped diamond synthesized from detonation nanodiamond in a C-O-H fluid at high pressure and high temperature

Boron doped diamond (BDD) was synthesized under high pressure and high temperature (HPHT) of 7 GPa, 1230 °C in a short time of 10 s from a powder mixtures of detonation nanodiamond (DND), pentaerythritol C{sub 5}H{sub 8}(OH){sub 4} and amorphous boron. SEM, TEM, XRD, XPS, FTIR and Raman spectroscopy indicated that BDD nano- and micro-crystals have formed by consolidation of DND particles (4 nm in size). XRD showed the enlargement of crystallites size to 6–80 nm and the increase in diamond lattice parameter by 0.02–0.07% without appearance of any microstrains. Raman spectroscopy was used to estimate the content of boron atoms embedded in the diamond lattice. It was found that the Raman diamond peak shifts significantly from 1332 cm{sup −1} to 1290 cm{sup −1} without appearance of any non-diamond carbon. The correlation between Raman peak position, its width, and boron content in diamond is proposed. Hydrogenated diamond carbon in significant amount was detected by IR spectroscopy and XPS. Due to the doping with boron content of about 0.1 at%, the electrical conductivity of the diamond achieved approximately 0.2 Ω{sup −1} cm{sup −1}. Reaction mechanism of diamond growth (models of recrystallization and oriented attachment) is discussed, including the initial stages of pentaerythritol pyrolysis and thermal desorption of functional groups from the surface of DND particles with the generation of supercritical fluid of low-molecular substances (H{sub 2}O, CH{sub 4}, CO, CO{sub 2}, etc.), as well as byproducts formation (B{sub 2}O{sub 3}, B{sub 4}C). - Highlights: • Boron doped diamond BDD synthesized from detonation nanodiamond, C{sub 5}H{sub 8}(OH){sub 4} and boron. • BDD nano/microcrystals are studied by XRD, SEM, TEM, XPS, FTIR, Raman spectroscopy. • Diamond growth and transition to electro-conductive state (0.2 Ω{sup −1} cm{sup −1}) is observed. • Correlations of Raman bands position and width versus boron content are considered. • High-pressure/temperature reaction mechanism (stages, routes, products) is described.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2013

Degradation behavior of hydroxyapatite/poly(lactic-co-glycolic) acid nanocomposite in simulated body fluid

Graphical abstract: In this manuscript, we initiated a systematic study to investigate the effect of HA on thermal properties, inner structure, reduction of mechanical strength, surface morphology and the surface deposit of n-HA/PLGA composite with respect to the soaking time. The results showed that n-HA played an important role in improving the degradation behavior of n-HA/PLGA composite, which can accelerate the degradation of n-HA/PLGA composite and endow it with bioactivity, after n-HA was detached from PLGA during the degradation, so that n-HA/PLGA composite may have a more promising prospect of the clinical application than pure PLGA as bone fracture internal fixation materials, and the results would be of reference significance to predict the in vivo degradation and biological properties. - Highlights: • Effect of n-HA on degradation behavior of n-HA/PLGA composite was investigated. • Degradation behaviors of n-HA/PLGA and PLGA were carried out in SBF for 6 months. • Viscosity, thermal properties, inner structure and bending strength were tested. • n-HA can accelerate the degradation and endows it with bioactivity. - Abstract: To investigate the effect of hydroxyapatite(HA) on the degradation behavior of hydroxyapatite/poly(lactic-co-glycolic) acid (HA/PLGA) nanocomposite, the degradation experiment of n-HA/PLGA composite and pure PLGA were carried out by soaking in simulated body fluid(SBF) at 37 °C for 1, 2, 4 and 6 months. The change of intrinsic viscosity, thermal properties, inner structure, bending strength reduction, surface morphology and the surface deposit of n-HA/PLGA composite and pure PLGA with respect to the soaking time were investigated by means of UbbeloHde Viscometer, differential scanning calorimeter (DSC), scanning electron microscope(SEM), electromechanical universal tester, a conventional camera and X-ray diffraction (XRD). The results showed that n-HA played an important role in improving the degradation behavior of n-HA/PLGA composite, which can accelerate the degradation PLGA and endow it with bioactivity, after n-HA was detached from PLGA during the degradation, so that n-HA/PLGA composite may have a more promising prospect of the clinical application than pure PLGA as bone fracture internal fixation materials.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Carbonated hydroxyapatite starting from calcite and different orthophosphates under moderate hydrothermal conditions: Synthesis and surface reactivity in simulated body fluid

Highlights: • Carbonated apatite (CAP) could be easily obtained from CaCO{sub 3} and orthophosphates. • Highest CaCO{sub 3} dissolution and apatitic carbonate content were obtained with H{sub 3}PO{sub 4}. • A-B-type CAP was formed. • The synthesized CAP was thermally stable up to 1000 °C. • This CAP showed high biomineralization activity before and after thermal treatment. - Abstract: The one-step synthesis of carbonated hydroxyapatite (CAP) using calcite and different orthophosphates was investigated in a closed batch reactor. Only orthophosphoric acid could lead to the complete decomposition of calcite particles, when the reaction temperature was set at 80 °C. On the other hand, the reaction time and the dilution of the initial calcite suspension had no significant influence on the formation of the solid products. CAP was formed as the main crystalline calcium phosphate with the carbonate content in the range of 4.2–4.6 wt.%. The thermal decarbonation of the synthesized CAP started at 750 °C but it was only significant at 1000 °C under air atmosphere. This thermal decarbonation was total at 1200 °C or above. All CAP samples and products following thermal treatments were found bioactive in the test using simulated body fluid (SBF) solution.

36 MATERIALS SCIENCE↗
DOE OSTI2016

Application of Jacobian-free Newton–Krylov method in implicitly solving two-fluid six-equation two-phase flow problems: Implementation, validation and benchmark

This work represents a first-of-its-kind successful application to employ advanced numerical methods in solving realistic two-phase flow problems with two-fluid six-equation two-phase flow model. These advanced numerical methods include high-resolution spatial discretization scheme with staggered grids (high-order) fully implicit time integration schemes, and Jacobian-free Newton–Krylov (JFNK) method as the nonlinear solver. The computer code developed in this work has been extensively validated with existing experimental flow boiling data in vertical pipes and rod bundles, which cover wide ranges of experimental conditions, such as pressure, inlet mass flux, wall heat flux and exit void fraction. Additional code-to-code benchmark with the RELAP5-3D code further verifies the correct code implementation. The combined methods employed in this work exhibit strong robustness in solving two-phase flow problems even when phase appearance (boiling) and realistic discrete flow regimes are considered. Transitional flow regimes used in existing system analysis codes, normally introduced to overcome numerical difficulty, were completely removed in this work. As a result, this in turn provides the possibility to utilize more sophisticated flow regime maps in the future to further improve simulation accuracy.

42 ENGINEERING↗
DOE OSTI2021

Alteration in micro-mechanical characteristics of wellbore cement fracture surfaces due to fluid exposure

The stiffness of wellbore cement fracture surfaces was measured after exposing to the advective flow of nitrogen, silicone oil, and medium sweet dead crude oil for different exposure periods. The test specimens were extracted from fractured cement cylinders, where the cement fracture surfaces were exposed to the different fluids up to 15 weeks. A nanoindenter with a Berkovich indenter tip was used to measure load-indentation depth data, which was used to extract the elastic modulus (E) and nano-hardness (H) of the cement fracture surfaces. A reduction in the elastic modulus compared with an unexposed specimen were observed in all the specimens. Both elastic modulus and nano-hardness for the specimens exposed to silicone oil were lower than specimens exposed to nitrogen gas and varied with the period of exposure. The elastic modulus and nano-hardness of the specimens exposed to crude oil were the lowest with a significant decrement with the exposure period. The frequency distribution of the nanoindentation measurements shows that the volume-fraction ratio of the two types of cement hydrated nanocomposites for both the unexposed and test specimens is about 70:30%. Phase transformation beneath the indenter is observed for all of the specimens, with more obvious plastic deformation in specimens exposed to crude oil. Analytical measurements (SEM, EDS, FT-IR, and XRD) on exposed cement fracture surfaces reveal different levels of physical and chemical alteration that are consistent with the reduction in stiffness measured by nanoindentation. Furthermore, the study suggests that cement stiffness will decrease due to crude oil exposure, and the fracture will be sensitive to stress and pore pressure with time.

42 ENGINEERING↗