Computational Models for Nanoscale Fluid Dynamics and Transport Inspired by Nonequilibrium Thermodynamics1

PubMed Central

Radhakrishnan, Ravi; Yu, Hsiu-Yu; Eckmann, David M.; Ayyaswamy, Portonovo S.

2017-01-01

Traditionally, the numerical computation of particle motion in a fluid is resolved through computational fluid dynamics (CFD). However, resolving the motion of nanoparticles poses additional challenges due to the coupling between the Brownian and hydrodynamic forces. Here, we focus on the Brownian motion of a nanoparticle coupled to adhesive interactions and confining-wall-mediated hydrodynamic interactions. We discuss several techniques that are founded on the basis of combining CFD methods with the theory of nonequilibrium statistical mechanics in order to simultaneously conserve thermal equipartition and to show correct hydrodynamic correlations. These include the fluctuating hydrodynamics (FHD) method, the generalized Langevin method, the hybrid method, and the deterministic method. Through the examples discussed, we also show a top-down multiscale progression of temporal dynamics from the colloidal scales to the molecular scales, and the associated fluctuations, hydrodynamic correlations. While the motivation and the examples discussed here pertain to nanoscale fluid dynamics and mass transport, the methodologies presented are rather general and can be easily adopted to applications in convective heat transfer. PMID:28035168

Simultaneous Multiple-Location Separation Control

NASA Technical Reports Server (NTRS)

Greenblatt, David (Inventor)

2009-01-01

A method of controlling a shear layer for a fluid dynamic body introduces first periodic disturbances into the fluid medium at a first flow separation location. Simultaneously, second periodic disturbances are introduced into the fluid medium at a second flow separation location. A phase difference between the first and second periodic disturbances is adjusted to control flow separation of the shear layer as the fluid medium moves over the fluid dynamic body.

Computational fluid mechanics utilizing the variational principle of modeling damping seals

NASA Technical Reports Server (NTRS)

Abernathy, J. M.

1986-01-01

A computational fluid dynamics code for application to traditional incompressible flow problems has been developed. The method is actually a slight compressibility approach which takes advantage of the bulk modulus and finite sound speed of all real fluids. The finite element numerical analog uses a dynamic differencing scheme based, in part, on a variational principle for computational fluid dynamics. The code was developed in order to study the feasibility of damping seals for high speed turbomachinery. Preliminary seal analyses have been performed.

Spectral methods for CFD

NASA Technical Reports Server (NTRS)

Zang, Thomas A.; Streett, Craig L.; Hussaini, M. Yousuff

1989-01-01

One of the objectives of these notes is to provide a basic introduction to spectral methods with a particular emphasis on applications to computational fluid dynamics. Another objective is to summarize some of the most important developments in spectral methods in the last two years. The fundamentals of spectral methods for simple problems will be covered in depth, and the essential elements of several fluid dynamical applications will be sketched.

Multipole Vortex Blobs (MVB): Symplectic Geometry and Dynamics.

PubMed

Holm, Darryl D; Jacobs, Henry O

2017-01-01

Vortex blob methods are typically characterized by a regularization length scale, below which the dynamics are trivial for isolated blobs. In this article, we observe that the dynamics need not be trivial if one is willing to consider distributional derivatives of Dirac delta functionals as valid vorticity distributions. More specifically, a new singular vortex theory is presented for regularized Euler fluid equations of ideal incompressible flow in the plane. We determine the conditions under which such regularized Euler fluid equations may admit vorticity singularities which are stronger than delta functions, e.g., derivatives of delta functions. We also describe the symplectic geometry associated with these augmented vortex structures, and we characterize the dynamics as Hamiltonian. Applications to the design of numerical methods similar to vortex blob methods are also discussed. Such findings illuminate the rich dynamics which occur below the regularization length scale and enlighten our perspective on the potential for regularized fluid models to capture multiscale phenomena.

Smoothed dissipative particle dynamics with angular momentum conservation

NASA Astrophysics Data System (ADS)

Müller, Kathrin; Fedosov, Dmitry A.; Gompper, Gerhard

2015-01-01

Smoothed dissipative particle dynamics (SDPD) combines two popular mesoscopic techniques, the smoothed particle hydrodynamics and dissipative particle dynamics (DPD) methods, and can be considered as an improved dissipative particle dynamics approach. Despite several advantages of the SDPD method over the conventional DPD model, the original formulation of SDPD by Español and Revenga (2003) [9], lacks angular momentum conservation, leading to unphysical results for problems where the conservation of angular momentum is essential. To overcome this limitation, we extend the SDPD method by introducing a particle spin variable such that local and global angular momentum conservation is restored. The new SDPD formulation (SDPD+a) is directly derived from the Navier-Stokes equation for fluids with spin, while thermal fluctuations are incorporated similarly to the DPD method. We test the new SDPD method and demonstrate that it properly reproduces fluid transport coefficients. Also, SDPD with angular momentum conservation is validated using two problems: (i) the Taylor-Couette flow with two immiscible fluids and (ii) a tank-treading vesicle in shear flow with a viscosity contrast between inner and outer fluids. For both problems, the new SDPD method leads to simulation predictions in agreement with the corresponding analytical theories, while the original SDPD method fails to capture properly physical characteristics of the systems due to violation of angular momentum conservation. In conclusion, the extended SDPD method with angular momentum conservation provides a new approach to tackle fluid problems such as multiphase flows and vesicle/cell suspensions, where the conservation of angular momentum is essential.

Hybrid Method for Power Control Simulation of a Single Fluid Plasma Thruster

NASA Astrophysics Data System (ADS)

Jaisankar, S.; Sheshadri, T. S.

2018-05-01

Propulsive plasma flow through a cylindrical-conical diverging thruster is simulated by a power controlled hybrid method to obtain the basic flow, thermodynamic and electromagnetic variables. Simulation is based on a single fluid model with electromagnetics being described by the equations of potential Poisson, Maxwell and the Ohm's law while the compressible fluid dynamics by the Navier Stokes in cylindrical form. The proposed method solved the electromagnetics and fluid dynamics separately, both to segregate the two prominent scales for an efficient computation and for the delivery of voltage controlled rated power. The magnetic transport is solved for steady state while fluid dynamics is allowed to evolve in time along with an electromagnetic source using schemes based on generalized finite difference discretization. The multistep methodology with power control is employed for simulating fully ionized propulsive flow of argon plasma through the thruster. Numerical solution shows convergence of every part of the solver including grid stability causing the multistep hybrid method to converge for a rated power delivery. Simulation results are reasonably in agreement with the reported physics of plasma flow in the thruster thus indicating the potential utility of this hybrid computational framework, especially when single fluid approximation of plasma is relevant.

Fluid-structure interaction dynamic simulation of spring-loaded pressure relief valves under seismic wave

NASA Astrophysics Data System (ADS)

Lv, Dongwei; Zhang, Jian; Yu, Xinhai

2018-05-01

In this paper, a fluid-structure interaction dynamic simulation method of spring-loaded pressure relief valve was established. The dynamic performances of the fluid regions and the stress and strain of the structure regions were calculated at the same time by accurately setting up the contact pairs between the solid parts and the coupling surfaces between the fluid regions and the structure regions. A two way fluid-structure interaction dynamic simulation of a simplified pressure relief valve model was carried out. The influence of vertical sinusoidal seismic waves on the performance of the pressure relief valve was preliminarily investigated by loading sine waves. Under vertical seismic waves, the pressure relief valve will flutter, and the reseating pressure was affected by the amplitude and frequency of the seismic waves. This simulation method of the pressure relief valve under vertical seismic waves can provide effective means for investigating the seismic performances of the valves, and make up for the shortcomings of the experiment.

Finite element analysis in fluids; Proceedings of the Seventh International Conference on Finite Element Methods in Flow Problems, University of Alabama, Huntsville, Apr. 3-7, 1989

NASA Technical Reports Server (NTRS)

Chung, T. J. (Editor); Karr, Gerald R. (Editor)

1989-01-01

Recent advances in computational fluid dynamics are examined in reviews and reports, with an emphasis on finite-element methods. Sections are devoted to adaptive meshes, atmospheric dynamics, combustion, compressible flows, control-volume finite elements, crystal growth, domain decomposition, EM-field problems, FDM/FEM, and fluid-structure interactions. Consideration is given to free-boundary problems with heat transfer, free surface flow, geophysical flow problems, heat and mass transfer, high-speed flow, incompressible flow, inverse design methods, MHD problems, the mathematics of finite elements, and mesh generation. Also discussed are mixed finite elements, multigrid methods, non-Newtonian fluids, numerical dissipation, parallel vector processing, reservoir simulation, seepage, shallow-water problems, spectral methods, supercomputer architectures, three-dimensional problems, and turbulent flows.

Freezing point and solid-liquid interfacial free energy of Stockmayer dipolar fluids: a molecular dynamics simulation study.

PubMed

Wang, Jun; Apte, Pankaj A; Morris, James R; Zeng, Xiao Cheng

2013-09-21

Stockmayer fluids are a prototype model system for dipolar fluids. We have computed the freezing temperatures of Stockmayer fluids at zero pressure using three different molecular-dynamics simulation methods, namely, the superheating-undercooling method, the constant-pressure and constant-temperature two-phase coexistence method, and the constant-pressure and constant-enthalpy two-phase coexistence method. The best estimate of the freezing temperature (in reduced unit) for the Stockmayer (SM) fluid with the dimensionless dipole moment μ*=1, √2, √3 is 0.656 ± 0.001, 0.726 ± 0.002, and 0.835 ± 0.005, respectively. The freezing temperature increases with the dipolar strength. Moreover, for the first time, the solid-liquid interfacial free energies γ of the fcc (111), (110), and (100) interfaces are computed using two independent methods, namely, the cleaving-wall method and the interfacial fluctuation method. Both methods predict that the interfacial free energy increases with the dipole moment. Although the interfacial fluctuation method suggests a weaker interfacial anisotropy, particularly for strongly dipolar SM fluids, both methods predicted the same trend of interfacial anisotropy, i.e., γ100 > γ110 > γ111.

A mesoscopic bridging scale method for fluids and coupling dissipative particle dynamics with continuum finite element method

PubMed Central

Kojic, Milos; Filipovic, Nenad; Tsuda, Akira

2012-01-01

A multiscale procedure to couple a mesoscale discrete particle model and a macroscale continuum model of incompressible fluid flow is proposed in this study. We call this procedure the mesoscopic bridging scale (MBS) method since it is developed on the basis of the bridging scale method for coupling molecular dynamics and finite element models [G.J. Wagner, W.K. Liu, Coupling of atomistic and continuum simulations using a bridging scale decomposition, J. Comput. Phys. 190 (2003) 249–274]. We derive the governing equations of the MBS method and show that the differential equations of motion of the mesoscale discrete particle model and finite element (FE) model are only coupled through the force terms. Based on this coupling, we express the finite element equations which rely on the Navier–Stokes and continuity equations, in a way that the internal nodal FE forces are evaluated using viscous stresses from the mesoscale model. The dissipative particle dynamics (DPD) method for the discrete particle mesoscale model is employed. The entire fluid domain is divided into a local domain and a global domain. Fluid flow in the local domain is modeled with both DPD and FE method, while fluid flow in the global domain is modeled by the FE method only. The MBS method is suitable for modeling complex (colloidal) fluid flows, where continuum methods are sufficiently accurate only in the large fluid domain, while small, local regions of particular interest require detailed modeling by mesoscopic discrete particles. Solved examples – simple Poiseuille and driven cavity flows illustrate the applicability of the proposed MBS method. PMID:23814322

The nonlinear dynamics of a spacecraft coupled to the vibration of a contained fluid

NASA Technical Reports Server (NTRS)

Peterson, Lee D.; Crawley, Edward F.; Hansman, R. John

1988-01-01

The dynamics of a linear spacecraft mode coupled to a nonlinear low gravity slosh of a fluid in a cylindrical tank is investigated. Coupled, nonlinear equations of motion for the fluid-spacecraft dynamics are derived through an assumed mode Lagrangian method. Unlike linear fluid slosh models, this nonlinear slosh model retains two fundamental slosh modes and three secondary modes. An approximate perturbation solution of the equations of motion indicates that the nonlinear coupled system response involves fluid-spacecraft modal resonances not predicted by either a linear, or a nonlinear, uncoupled slosh analysis. Experimental results substantiate the analytical predictions.

System and method for reducing combustion dynamics in a combustor

DOEpatents

Uhm, Jong Ho; Ziminsky, Willy Steve; Johnson, Thomas Edward; Srinivasan, Shiva; York, William David

2016-11-29

A system for reducing combustion dynamics in a combustor includes an end cap that extends radially across the combustor and includes an upstream surface axially separated from a downstream surface. A combustion chamber is downstream of the end cap, and tubes extend from the upstream surface through the downstream surface. Each tube provides fluid communication through the end cap to the combustion chamber. The system further includes means for reducing combustion dynamics in the combustor. A method for reducing combustion dynamics in a combustor includes flowing a working fluid through tubes that extend axially through an end cap that extends radially across the combustor and obstructing at least a portion of the working fluid flowing through a first set of the tubes.

Modeling and control of magnetorheological fluid dampers using neural networks

NASA Astrophysics Data System (ADS)

Wang, D. H.; Liao, W. H.

2005-02-01

Due to the inherent nonlinear nature of magnetorheological (MR) fluid dampers, one of the challenging aspects for utilizing these devices to achieve high system performance is the development of accurate models and control algorithms that can take advantage of their unique characteristics. In this paper, the direct identification and inverse dynamic modeling for MR fluid dampers using feedforward and recurrent neural networks are studied. The trained direct identification neural network model can be used to predict the damping force of the MR fluid damper on line, on the basis of the dynamic responses across the MR fluid damper and the command voltage, and the inverse dynamic neural network model can be used to generate the command voltage according to the desired damping force through supervised learning. The architectures and the learning methods of the dynamic neural network models and inverse neural network models for MR fluid dampers are presented, and some simulation results are discussed. Finally, the trained neural network models are applied to predict and control the damping force of the MR fluid damper. Moreover, validation methods for the neural network models developed are proposed and used to evaluate their performance. Validation results with different data sets indicate that the proposed direct identification dynamic model using the recurrent neural network can be used to predict the damping force accurately and the inverse identification dynamic model using the recurrent neural network can act as a damper controller to generate the command voltage when the MR fluid damper is used in a semi-active mode.

Advances in modelling of biomimetic fluid flow at different scales

PubMed Central

2011-01-01

The biomimetic flow at different scales has been discussed at length. The need of looking into the biological surfaces and morphologies and both geometrical and physical similarities to imitate the technological products and processes has been emphasized. The complex fluid flow and heat transfer problems, the fluid-interface and the physics involved at multiscale and macro-, meso-, micro- and nano-scales have been discussed. The flow and heat transfer simulation is done by various CFD solvers including Navier-Stokes and energy equations, lattice Boltzmann method and molecular dynamics method. Combined continuum-molecular dynamics method is also reviewed. PMID:21711847

An immersed-shell method for modelling fluid–structure interactions

PubMed Central

Viré, A.; Xiang, J.; Pain, C. C.

2015-01-01

The paper presents a novel method for numerically modelling fluid–structure interactions. The method consists of solving the fluid-dynamics equations on an extended domain, where the computational mesh covers both fluid and solid structures. The fluid and solid velocities are relaxed to one another through a penalty force. The latter acts on a thin shell surrounding the solid structures. Additionally, the shell is represented on the extended domain by a non-zero shell-concentration field, which is obtained by conservatively mapping the shell mesh onto the extended mesh. The paper outlines the theory underpinning this novel method, referred to as the immersed-shell approach. It also shows how the coupling between a fluid- and a structural-dynamics solver is achieved. At this stage, results are shown for cases of fundamental interest. PMID:25583857

Simulating coupled dynamics of a rigid-flexible multibody system and compressible fluid

NASA Astrophysics Data System (ADS)

Hu, Wei; Tian, Qiang; Hu, HaiYan

2018-04-01

As a subsequent work of previous studies of authors, a new parallel computation approach is proposed to simulate the coupled dynamics of a rigid-flexible multibody system and compressible fluid. In this approach, the smoothed particle hydrodynamics (SPH) method is used to model the compressible fluid, the natural coordinate formulation (NCF) and absolute nodal coordinate formulation (ANCF) are used to model the rigid and flexible bodies, respectively. In order to model the compressible fluid properly and efficiently via SPH method, three measures are taken as follows. The first is to use the Riemann solver to cope with the fluid compressibility, the second is to define virtual particles of SPH to model the dynamic interaction between the fluid and the multibody system, and the third is to impose the boundary conditions of periodical inflow and outflow to reduce the number of SPH particles involved in the computation process. Afterwards, a parallel computation strategy is proposed based on the graphics processing unit (GPU) to detect the neighboring SPH particles and to solve the dynamic equations of SPH particles in order to improve the computation efficiency. Meanwhile, the generalized-alpha algorithm is used to solve the dynamic equations of the multibody system. Finally, four case studies are given to validate the proposed parallel computation approach.

Models for the dynamics of dust-like matter in the self-gravity field: The method of hydrodynamic substitutions

NASA Astrophysics Data System (ADS)

Zhuravlev, V. M.

2017-09-01

Models for the dynamics of a dust-like medium in the self-gravity field are investigated. Solutions of the corresponding problems are constructed by the method of hydrodynamic substitutions generalizing the Cole-Hopf substitutions. The method is extended to multidimensional ideal and viscous fluid flows with cylindrical and spherical symmetries for which exact solutions are constructed. Solutions for the dynamics of self-gravitating dust with arbitrary initial distributions of both fluid density and velocity are constructed using special coordinate transformations. In particular, the problem of cosmological expansion is considered in terms of Newton's gravity theory. Models of a one-dimensional viscous dust fluid flow and some problems of gas hydrodynamics are considered. Examples of exact solutions and their brief analysis are provided.

Computational Fluid Dynamics Symposium on Aeropropulsion

NASA Technical Reports Server (NTRS)

1991-01-01

Recognizing the considerable advances that have been made in computational fluid dynamics, the Internal Fluid Mechanics Division of NASA Lewis Research Center sponsored this symposium with the objective of providing a forum for exchanging information regarding recent developments in numerical methods, physical and chemical modeling, and applications. This conference publication is a compilation of 4 invited and 34 contributed papers presented in six sessions: algorithms one and two, turbomachinery, turbulence, components application, and combustors. Topics include numerical methods, grid generation, chemically reacting flows, turbulence modeling, inlets, nozzles, and unsteady flows.

Fluid technology (selected components, devices, and systems): A compilation

NASA Technical Reports Server (NTRS)

1974-01-01

Developments in fluid technology and hydraulic equipment are presented. The subjects considered are: (1) the use of fluids in the operation of switches, amplifiers, and servo devices, (2) devices and data for laboratory use in the study of fluid dynamics, and (3) the use of fluids as controls and certain methods of controlling fluids.

Scaling effects in spiral capsule robots.

PubMed

Liang, Liang; Hu, Rong; Chen, Bai; Tang, Yong; Xu, Yan

2017-04-01

Spiral capsule robots can be applied to human gastrointestinal tracts and blood vessels. Because of significant variations in the sizes of the inner diameters of the intestines as well as blood vessels, this research has been unable to meet the requirements for medical applications. By applying the fluid dynamic equations, using the computational fluid dynamics method, to a robot axial length ranging from 10 -5 to 10 -2  m, the operational performance indicators (axial driving force, load torque, and maximum fluid pressure on the pipe wall) of the spiral capsule robot and the fluid turbulent intensity around the robot spiral surfaces was numerically calculated in a straight rigid pipe filled with fluid. The reasonableness and validity of the calculation method adopted in this study were verified by the consistency of the calculated values by the computational fluid dynamics method and the experimental values from a relevant literature. The results show that the greater the fluid turbulent intensity, the greater the impact of the fluid turbulence on the driving performance of the spiral capsule robot and the higher the energy consumption of the robot. For the same level of size of the robot, the axial driving force, the load torque, and the maximum fluid pressure on the pipe wall of the outer spiral robot were larger than those of the inner spiral robot. For different requirements of the operating environment, we can choose a certain kind of spiral capsule robot. This study provides a theoretical foundation for spiral capsule robots.

Improving students’ conceptions on fluid dynamics through peer teaching model with PDEODE (PTM-PDEODE)

NASA Astrophysics Data System (ADS)

Samsudin, A.; Fratiwi, N.; Amin, N.; Wiendartun; Supriyatman; Wibowo, F.; Faizin, M.; Costu, B.

2018-05-01

This study based on an importance of improving students’ conceptions and reduces students’ misconceptions on fluid dynamics concepts. Consequently, should be done the study through combining Peer Teaching Model (PTM) and PDEODE (Prediction, Discuss, Explain, Observe, Discuss and Explain) learning strategy (PTM-PDEODE). For the research methods, we used the 4D model (Defining, Designing, Developing, and Disseminating). The samples are 38 students (their ages were an average of 17 years-old) at one of the senior high schools in Bandung. The improvement of students’ conceptions was diagnosed through a four-tier test of fluid dynamics. At the disseminating phase, students’ conceptions of fluid dynamics concepts are increase after the use of PTM-PDEODE. In conclusion, the development of PTM-PDEODE is respectable enough to improve students’ conceptions on dinamics fluid.

An immersed boundary method for fluid-structure interaction with compressible multiphase flows

NASA Astrophysics Data System (ADS)

Wang, Li; Currao, Gaetano M. D.; Han, Feng; Neely, Andrew J.; Young, John; Tian, Fang-Bao

2017-10-01

This paper presents a two-dimensional immersed boundary method for fluid-structure interaction with compressible multiphase flows involving large structure deformations. This method involves three important parts: flow solver, structure solver and fluid-structure interaction coupling. In the flow solver, the compressible multiphase Navier-Stokes equations for ideal gases are solved by a finite difference method based on a staggered Cartesian mesh, where a fifth-order accuracy Weighted Essentially Non-Oscillation (WENO) scheme is used to handle spatial discretization of the convective term, a fourth-order central difference scheme is employed to discretize the viscous term, the third-order TVD Runge-Kutta scheme is used to discretize the temporal term, and the level-set method is adopted to capture the multi-material interface. In this work, the structure considered is a geometrically non-linear beam which is solved by using a finite element method based on the absolute nodal coordinate formulation (ANCF). The fluid dynamics and the structure motion are coupled in a partitioned iterative manner with a feedback penalty immersed boundary method where the flow dynamics is defined on a fixed Lagrangian grid and the structure dynamics is described on a global coordinate. We perform several validation cases (including fluid over a cylinder, structure dynamics, flow induced vibration of a flexible plate, deformation of a flexible panel induced by shock waves in a shock tube, an inclined flexible plate in a hypersonic flow, and shock-induced collapse of a cylindrical helium cavity in the air), and compare the results with experimental and other numerical data. The present results agree well with the published data and the current experiment. Finally, we further demonstrate the versatility of the present method by applying it to a flexible plate interacting with multiphase flows.

SPAR improved structure/fluid dynamic analysis capability

NASA Technical Reports Server (NTRS)

Oden, J. T.; Pearson, M. L.

1983-01-01

The capability of analyzing a coupled dynamic system of flowing fluid and elastic structure was added to the SPAR computer code. A method, developed and adopted for use in SPAR utilizes the existing assumed stress hybrid plan element in SPAR. An operational mode was incorporated in SPAR which provides the capability for analyzing the flaw of a two dimensional, incompressible, viscous fluid within rigid boundaries. Equations were developed to provide for the eventual analysis of the interaction of such fluids with an elastic solid.

Technical Competencies Applied in Experimental Fluid Dynamics

NASA Astrophysics Data System (ADS)

Tagg, Randall

2017-11-01

The practical design, construction, and operation of fluid dynamics experiments require a broad range of competencies. Three types are instrumental, procedural, and design. Respective examples would be operation of a spectrum analyzer, soft-soldering or brazing flow plumbing, and design of a small wind tunnel. Some competencies, such as the selection and installation of pumping systems, are unique to fluid dynamics and fluids engineering. Others, such as the design and construction of electronic amplifiers or optical imaging systems, overlap with other fields. Thus the identification and development of learning materials and methods for instruction are part of a larger effort to identify competencies needed in active research and technical innovation.

SPAR improved structure-fluid dynamic analysis capability, phase 2

NASA Technical Reports Server (NTRS)

Pearson, M. L.

1984-01-01

An efficient and general method of analyzing a coupled dynamic system of fluid flow and elastic structures is investigated. The improvement of Structural Performance Analysis and Redesign (SPAR) code is summarized. All error codes are documented and the SPAR processor/subroutine cross reference is included.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Müller, Kathrin, E-mail: k.mueller@fz-juelich.de; Fedosov, Dmitry A., E-mail: d.fedosov@fz-juelich.de; Gompper, Gerhard, E-mail: g.gompper@fz-juelich.de

Smoothed dissipative particle dynamics (SDPD) combines two popular mesoscopic techniques, the smoothed particle hydrodynamics and dissipative particle dynamics (DPD) methods, and can be considered as an improved dissipative particle dynamics approach. Despite several advantages of the SDPD method over the conventional DPD model, the original formulation of SDPD by Español and Revenga (2003) [9], lacks angular momentum conservation, leading to unphysical results for problems where the conservation of angular momentum is essential. To overcome this limitation, we extend the SDPD method by introducing a particle spin variable such that local and global angular momentum conservation is restored. The new SDPDmore » formulation (SDPD+a) is directly derived from the Navier–Stokes equation for fluids with spin, while thermal fluctuations are incorporated similarly to the DPD method. We test the new SDPD method and demonstrate that it properly reproduces fluid transport coefficients. Also, SDPD with angular momentum conservation is validated using two problems: (i) the Taylor–Couette flow with two immiscible fluids and (ii) a tank-treading vesicle in shear flow with a viscosity contrast between inner and outer fluids. For both problems, the new SDPD method leads to simulation predictions in agreement with the corresponding analytical theories, while the original SDPD method fails to capture properly physical characteristics of the systems due to violation of angular momentum conservation. In conclusion, the extended SDPD method with angular momentum conservation provides a new approach to tackle fluid problems such as multiphase flows and vesicle/cell suspensions, where the conservation of angular momentum is essential.« less

Dynamic modeling method for infrared smoke based on enhanced discrete phase model

NASA Astrophysics Data System (ADS)

Zhang, Zhendong; Yang, Chunling; Zhang, Yan; Zhu, Hongbo

2018-03-01

The dynamic modeling of infrared (IR) smoke plays an important role in IR scene simulation systems and its accuracy directly influences the system veracity. However, current IR smoke models cannot provide high veracity, because certain physical characteristics are frequently ignored in fluid simulation; simplifying the discrete phase as a continuous phase and ignoring the IR decoy missile-body spinning. To address this defect, this paper proposes a dynamic modeling method for IR smoke, based on an enhanced discrete phase model (DPM). A mathematical simulation model based on an enhanced DPM is built and a dynamic computing fluid mesh is generated. The dynamic model of IR smoke is then established using an extended equivalent-blackbody-molecule model. Experiments demonstrate that this model realizes a dynamic method for modeling IR smoke with higher veracity.

Smart algorithms and adaptive methods in computational fluid dynamics

NASA Astrophysics Data System (ADS)

Tinsley Oden, J.

1989-05-01

A review is presented of the use of smart algorithms which employ adaptive methods in processing large amounts of data in computational fluid dynamics (CFD). Smart algorithms use a rationally based set of criteria for automatic decision making in an attempt to produce optimal simulations of complex fluid dynamics problems. The information needed to make these decisions is not known beforehand and evolves in structure and form during the numerical solution of flow problems. Once the code makes a decision based on the available data, the structure of the data may change, and criteria may be reapplied in order to direct the analysis toward an acceptable end. Intelligent decisions are made by processing vast amounts of data that evolve unpredictably during the calculation. The basic components of adaptive methods and their application to complex problems of fluid dynamics are reviewed. The basic components of adaptive methods are: (1) data structures, that is what approaches are available for modifying data structures of an approximation so as to reduce errors; (2) error estimation, that is what techniques exist for estimating error evolution in a CFD calculation; and (3) solvers, what algorithms are available which can function in changing meshes. Numerical examples which demonstrate the viability of these approaches are presented.

Algorithmic trends in computational fluid dynamics; The Institute for Computer Applications in Science and Engineering (ICASE)/LaRC Workshop, NASA Langley Research Center, Hampton, VA, US, Sep. 15-17, 1991

NASA Technical Reports Server (NTRS)

Hussaini, M. Y. (Editor); Kumar, A. (Editor); Salas, M. D. (Editor)

1993-01-01

The purpose here is to assess the state of the art in the areas of numerical analysis that are particularly relevant to computational fluid dynamics (CFD), to identify promising new developments in various areas of numerical analysis that will impact CFD, and to establish a long-term perspective focusing on opportunities and needs. Overviews are given of discretization schemes, computational fluid dynamics, algorithmic trends in CFD for aerospace flow field calculations, simulation of compressible viscous flow, and massively parallel computation. Also discussed are accerelation methods, spectral and high-order methods, multi-resolution and subcell resolution schemes, and inherently multidimensional schemes.

Strongly Coupled Fluid-Body Dynamics in the Immersed Boundary Projection Method

NASA Astrophysics Data System (ADS)

Wang, Chengjie; Eldredge, Jeff D.

2014-11-01

A computational algorithm is developed to simulate dynamically coupled interaction between fluid and rigid bodies. The basic computational framework is built upon a multi-domain immersed boundary method library, whirl, developed in previous work. In this library, the Navier-Stokes equations for incompressible flow are solved on a uniform Cartesian grid by the vorticity-based immersed boundary projection method of Colonius and Taira. A solver for the dynamics of rigid-body systems is also included. The fluid and rigid-body solvers are strongly coupled with an iterative approach based on the block Gauss-Seidel method. Interfacial force, with its intimate connection with the Lagrange multipliers used in the fluid solver, is used as the primary iteration variable. Relaxation, developed from a stability analysis of the iterative scheme, is used to achieve convergence in only 2-4 iterations per time step. Several two- and three-dimensional numerical tests are conducted to validate and demonstrate the method, including flapping of flexible wings, self-excited oscillations of a system of linked plates and three-dimensional propulsion of flexible fluked tail. This work has been supported by AFOSR, under Award FA9550-11-1-0098.

Interpolation Method Needed for Numerical Uncertainty Analysis of Computational Fluid Dynamics

NASA Technical Reports Server (NTRS)

Groves, Curtis; Ilie, Marcel; Schallhorn, Paul

2014-01-01

Using Computational Fluid Dynamics (CFD) to predict a flow field is an approximation to the exact problem and uncertainties exist. There is a method to approximate the errors in CFD via Richardson's Extrapolation. This method is based off of progressive grid refinement. To estimate the errors in an unstructured grid, the analyst must interpolate between at least three grids. This paper describes a study to find an appropriate interpolation scheme that can be used in Richardson's extrapolation or other uncertainty method to approximate errors. Nomenclature

Least-squares finite element method for fluid dynamics

NASA Technical Reports Server (NTRS)

Jiang, Bo-Nan; Povinelli, Louis A.

1989-01-01

An overview is given of new developments of the least squares finite element method (LSFEM) in fluid dynamics. Special emphasis is placed on the universality of LSFEM; the symmetry and positiveness of the algebraic systems obtained from LSFEM; the accommodation of LSFEM to equal order interpolations for incompressible viscous flows; and the natural numerical dissipation of LSFEM for convective transport problems and high speed compressible flows. The performance of LSFEM is illustrated by numerical examples.

Entropic Lattice Boltzmann Methods

DTIC Science & Technology

2001-12-10

model of fluid dynamics in one dimension, first considered by Renda et al. in 1997 [14]. Here the geometric picture involves a four dimensional polytope...convention of including constant terms in an extra column of the matrix, using the device of appending 1 to the column vector of unknowns. In general, there...we apply the entropic lattice Boltzmann method to a simple five-velocity model of fluid dynamics in one dimension, first considered by Renda et al

Finite elements and fluid dynamics. [instability effects on solution of nonlinear equations

NASA Technical Reports Server (NTRS)

Fix, G.

1975-01-01

Difficulties concerning a use of the finite element method in the solution of the nonlinear equations of fluid dynamics are partly related to various 'hidden' instabilities which often arise in fluid calculations. The instabilities are typically due to boundary effects or nonlinearities. It is shown that in certain cases these instabilities can be avoided if certain conservation laws are satisfied, and that the latter are often intimately related to finite elements.

Blending and nudging in fluid dynamics: some simple observations

NASA Astrophysics Data System (ADS)

Germano, M.

2017-10-01

Blending and nudging methods have been recently applied in fluid dynamics, particularly regarding the assimilation of experimental data into the computations. In the paper we formally derive the differential equation associated to blending and compare it to the standard nudging equation. Some simple considerations related to these techniques and their mutual relations are exposed.

High Order Approximations for Compressible Fluid Dynamics on Unstructured and Cartesian Meshes

NASA Technical Reports Server (NTRS)

Barth, Timothy (Editor); Deconinck, Herman (Editor)

1999-01-01

The development of high-order accurate numerical discretization techniques for irregular domains and meshes is often cited as one of the remaining challenges facing the field of computational fluid dynamics. In structural mechanics, the advantages of high-order finite element approximation are widely recognized. This is especially true when high-order element approximation is combined with element refinement (h-p refinement). In computational fluid dynamics, high-order discretization methods are infrequently used in the computation of compressible fluid flow. The hyperbolic nature of the governing equations and the presence of solution discontinuities makes high-order accuracy difficult to achieve. Consequently, second-order accurate methods are still predominately used in industrial applications even though evidence suggests that high-order methods may offer a way to significantly improve the resolution and accuracy for these calculations. To address this important topic, a special course was jointly organized by the Applied Vehicle Technology Panel of NATO's Research and Technology Organization (RTO), the von Karman Institute for Fluid Dynamics, and the Numerical Aerospace Simulation Division at the NASA Ames Research Center. The NATO RTO sponsored course entitled "Higher Order Discretization Methods in Computational Fluid Dynamics" was held September 14-18, 1998 at the von Karman Institute for Fluid Dynamics in Belgium and September 21-25, 1998 at the NASA Ames Research Center in the United States. During this special course, lecturers from Europe and the United States gave a series of comprehensive lectures on advanced topics related to the high-order numerical discretization of partial differential equations with primary emphasis given to computational fluid dynamics (CFD). Additional consideration was given to topics in computational physics such as the high-order discretization of the Hamilton-Jacobi, Helmholtz, and elasticity equations. This volume consists of five articles prepared by the special course lecturers. These articles should be of particular relevance to those readers with an interest in numerical discretization techniques which generalize to very high-order accuracy. The articles of Professors Abgrall and Shu consider the mathematical formulation of high-order accurate finite volume schemes utilizing essentially non-oscillatory (ENO) and weighted essentially non-oscillatory (WENO) reconstruction together with upwind flux evaluation. These formulations are particularly effective in computing numerical solutions of conservation laws containing solution discontinuities. Careful attention is given by the authors to implementational issues and techniques for improving the overall efficiency of these methods. The article of Professor Cockburn discusses the discontinuous Galerkin finite element method. This method naturally extends to high-order accuracy and has an interpretation as a finite volume method. Cockburn addresses two important issues associated with the discontinuous Galerkin method: controlling spurious extrema near solution discontinuities via "limiting" and the extension to second order advective-diffusive equations (joint work with Shu). The articles of Dr. Henderson and Professor Schwab consider the mathematical formulation and implementation of the h-p finite element methods using hierarchical basis functions and adaptive mesh refinement. These methods are particularly useful in computing high-order accurate solutions containing perturbative layers and corner singularities. Additional flexibility is obtained using a mortar FEM technique whereby nonconforming elements are interfaced together. Numerous examples are given by Henderson applying the h-p FEM method to the simulation of turbulence and turbulence transition.

The Influence of Dynamic Contact Angle on Wetting Dynamics

NASA Technical Reports Server (NTRS)

Rame, Enrique; Garoff, Steven

2005-01-01

When surface tension forces dominate, and regardless of whether the situation is static or dynamic, the contact angle (the angle the interface between two immiscible fluids makes when it contacts a solid) is the key parameter that determines the shape of a fluid-fluid interface. The static contact angle is easy to measure and implement in models predicting static capillary surface shapes and such associated quantities as pressure drops. By contrast, when the interface moves relative to the solid (as in dynamic wetting processes) the dynamic contact angle is not identified unambiguously because it depends on the geometry of the system Consequently, its determination becomes problematic and measurements in one geometry cannot be applied in another for prediction purposes. However, knowing how to measure and use the dynamic contact angle is crucial to determine such dynamics as a microsystem throughput reliably. In this talk we will present experimental and analytical efforts aimed at resolving modeling issues present in dynamic wetting. We will review experiments that show the inadequacy of the usual hydrodynamic model when a fluid-fluid meniscus moves over a solid surface such as the wall of a small tube or duct. We will then present analytical results that show how to parametrize these problems in a predictive manner. We will illustrate these ideas by showing how to implement the method in numerical fluid mechanical calculations.

ISCFD Nagoya 1989 - International Symposium on Computational Fluid Dynamics, 3rd, Nagoya, Japan, Aug. 28-31, 1989, Technical Papers

NASA Astrophysics Data System (ADS)

Recent advances in computational fluid dynamics are discussed in reviews and reports. Topics addressed include large-scale LESs for turbulent pipe and channel flows, numerical solutions of the Euler and Navier-Stokes equations on parallel computers, multigrid methods for steady high-Reynolds-number flow past sudden expansions, finite-volume methods on unstructured grids, supersonic wake flow on a blunt body, a grid-characteristic method for multidimensional gas dynamics, and CIC numerical simulation of a wave boundary layer. Consideration is given to vortex simulations of confined two-dimensional jets, supersonic viscous shear layers, spectral methods for compressible flows, shock-wave refraction at air/water interfaces, oscillatory flow in a two-dimensional collapsible channel, the growth of randomness in a spatially developing wake, and an efficient simplex algorithm for the finite-difference and dynamic linear-programming method in optimal potential control.

Axisymmetric Lattice Boltzmann Model of Droplet Impact on Solid Surfaces

NASA Astrophysics Data System (ADS)

Dalgamoni, Hussein; Yong, Xin

2017-11-01

Droplet impact is a ubiquitous fluid phenomena encountered in scientific and engineering applications such as ink-jet printing, coating, electronics manufacturing, and many others. It is of great technological importance to understand the detailed dynamics of drop impact on various surfaces. The lattice Boltzmann method (LBM) emerges as an efficient method for modeling complex fluid systems involving rapidly evolving fluid-fluid and fluid-solid interfaces with complex geometries. In this work, we model droplet impact on flat solid substrates with well-defined wetting behavior using a two-phase axisymmetric LBM with high density and viscosity contrasts. We extend the two-dimensional Lee and Liu model to capture axisymmetric effect in the normal impact. First we compare the 2D axisymmetric results with the 2D and 3D results reported by Lee and Liu to probe the effect of axisymmetric terms. Then, we explore the effects of Weber number, Ohnesorge number, and droplet-surface equilibrium contact angle on the impact. The dynamic contact angle and spreading factor of the droplet during impact are investigated to qualitatively characterize the impact dynamics.

Dynamic Deployment Simulations of Inflatable Space Structures

NASA Technical Reports Server (NTRS)

Wang, John T.

2005-01-01

The feasibility of using Control Volume (CV) method and the Arbitrary Lagrangian Eulerian (ALE) method in LSDYNA to simulate the dynamic deployment of inflatable space structures is investigated. The CV and ALE methods were used to predict the inflation deployments of three folded tube configurations. The CV method was found to be a simple and computationally efficient method that may be adequate for modeling slow inflation deployment sine the inertia of the inflation gas can be neglected. The ALE method was found to be very computationally intensive since it involves the solving of three conservative equations of fluid as well as dealing with complex fluid structure interactions.

Computational fluid dynamics - An introduction for engineers

NASA Astrophysics Data System (ADS)

Abbott, Michael Barry; Basco, David R.

An introduction to the fundamentals of CFD for engineers and physical scientists is presented. The principal definitions, basic ideas, and most common methods used in CFD are presented, and the application of these methods to the description of free surface, unsteady, and turbulent flow is shown. Emphasis is on the numerical treatment of incompressible unsteady fluid flow with primary applications to water problems using the finite difference method. While traditional areas of application like hydrology, hydraulic and coastal engineering and oceanography get the main emphasis, newer areas of application such as medical fluid dynamics, bioengineering, and soil physics and chemistry are also addressed. The possibilities and limitations of CFD are pointed out along with the relations of CFD to other branches of science.

Implicit mesh discontinuous Galerkin methods and interfacial gauge methods for high-order accurate interface dynamics, with applications to surface tension dynamics, rigid body fluid-structure interaction, and free surface flow: Part I

NASA Astrophysics Data System (ADS)

Saye, Robert

2017-09-01

In this two-part paper, a high-order accurate implicit mesh discontinuous Galerkin (dG) framework is developed for fluid interface dynamics, facilitating precise computation of interfacial fluid flow in evolving geometries. The framework uses implicitly defined meshes-wherein a reference quadtree or octree grid is combined with an implicit representation of evolving interfaces and moving domain boundaries-and allows physically prescribed interfacial jump conditions to be imposed or captured with high-order accuracy. Part one discusses the design of the framework, including: (i) high-order quadrature for implicitly defined elements and faces; (ii) high-order accurate discretisation of scalar and vector-valued elliptic partial differential equations with interfacial jumps in ellipticity coefficient, leading to optimal-order accuracy in the maximum norm and discrete linear systems that are symmetric positive (semi)definite; (iii) the design of incompressible fluid flow projection operators, which except for the influence of small penalty parameters, are discretely idempotent; and (iv) the design of geometric multigrid methods for elliptic interface problems on implicitly defined meshes and their use as preconditioners for the conjugate gradient method. Also discussed is a variety of aspects relating to moving interfaces, including: (v) dG discretisations of the level set method on implicitly defined meshes; (vi) transferring state between evolving implicit meshes; (vii) preserving mesh topology to accurately compute temporal derivatives; (viii) high-order accurate reinitialisation of level set functions; and (ix) the integration of adaptive mesh refinement. In part two, several applications of the implicit mesh dG framework in two and three dimensions are presented, including examples of single phase flow in nontrivial geometry, surface tension-driven two phase flow with phase-dependent fluid density and viscosity, rigid body fluid-structure interaction, and free surface flow. A class of techniques known as interfacial gauge methods is adopted to solve the corresponding incompressible Navier-Stokes equations, which, compared to archetypical projection methods, have a weaker coupling between fluid velocity, pressure, and interface position, and allow high-order accurate numerical methods to be developed more easily. Convergence analyses conducted throughout the work demonstrate high-order accuracy in the maximum norm for all of the applications considered; for example, fourth-order spatial accuracy in fluid velocity, pressure, and interface location is demonstrated for surface tension-driven two phase flow in 2D and 3D. Specific application examples include: vortex shedding in nontrivial geometry, capillary wave dynamics revealing fine-scale flow features, falling rigid bodies tumbling in unsteady flow, and free surface flow over a submersed obstacle, as well as high Reynolds number soap bubble oscillation dynamics and vortex shedding induced by a type of Plateau-Rayleigh instability in water ripple free surface flow. These last two examples compare numerical results with experimental data and serve as an additional means of validation; they also reveal physical phenomena not visible in the experiments, highlight how small-scale interfacial features develop and affect macroscopic dynamics, and demonstrate the wide range of spatial scales often at play in interfacial fluid flow.

Implicit mesh discontinuous Galerkin methods and interfacial gauge methods for high-order accurate interface dynamics, with applications to surface tension dynamics, rigid body fluid-structure interaction, and free surface flow: Part II

NASA Astrophysics Data System (ADS)

Saye, Robert

2017-09-01

In this two-part paper, a high-order accurate implicit mesh discontinuous Galerkin (dG) framework is developed for fluid interface dynamics, facilitating precise computation of interfacial fluid flow in evolving geometries. The framework uses implicitly defined meshes-wherein a reference quadtree or octree grid is combined with an implicit representation of evolving interfaces and moving domain boundaries-and allows physically prescribed interfacial jump conditions to be imposed or captured with high-order accuracy. Part one discusses the design of the framework, including: (i) high-order quadrature for implicitly defined elements and faces; (ii) high-order accurate discretisation of scalar and vector-valued elliptic partial differential equations with interfacial jumps in ellipticity coefficient, leading to optimal-order accuracy in the maximum norm and discrete linear systems that are symmetric positive (semi)definite; (iii) the design of incompressible fluid flow projection operators, which except for the influence of small penalty parameters, are discretely idempotent; and (iv) the design of geometric multigrid methods for elliptic interface problems on implicitly defined meshes and their use as preconditioners for the conjugate gradient method. Also discussed is a variety of aspects relating to moving interfaces, including: (v) dG discretisations of the level set method on implicitly defined meshes; (vi) transferring state between evolving implicit meshes; (vii) preserving mesh topology to accurately compute temporal derivatives; (viii) high-order accurate reinitialisation of level set functions; and (ix) the integration of adaptive mesh refinement. In part two, several applications of the implicit mesh dG framework in two and three dimensions are presented, including examples of single phase flow in nontrivial geometry, surface tension-driven two phase flow with phase-dependent fluid density and viscosity, rigid body fluid-structure interaction, and free surface flow. A class of techniques known as interfacial gauge methods is adopted to solve the corresponding incompressible Navier-Stokes equations, which, compared to archetypical projection methods, have a weaker coupling between fluid velocity, pressure, and interface position, and allow high-order accurate numerical methods to be developed more easily. Convergence analyses conducted throughout the work demonstrate high-order accuracy in the maximum norm for all of the applications considered; for example, fourth-order spatial accuracy in fluid velocity, pressure, and interface location is demonstrated for surface tension-driven two phase flow in 2D and 3D. Specific application examples include: vortex shedding in nontrivial geometry, capillary wave dynamics revealing fine-scale flow features, falling rigid bodies tumbling in unsteady flow, and free surface flow over a submersed obstacle, as well as high Reynolds number soap bubble oscillation dynamics and vortex shedding induced by a type of Plateau-Rayleigh instability in water ripple free surface flow. These last two examples compare numerical results with experimental data and serve as an additional means of validation; they also reveal physical phenomena not visible in the experiments, highlight how small-scale interfacial features develop and affect macroscopic dynamics, and demonstrate the wide range of spatial scales often at play in interfacial fluid flow.

Effect of ultrasound on dynamics characteristic of the cavitation bubble in grinding fluids during honing process.

PubMed

Guo, Ce; Zhu, Xijing

2018-03-01

The effect of ultrasound on generating and controlling the cavitation bubble of the grinding fluid during ultrasonic vibration honing was investigated. The grinding fluid on the surface of the honing stone was measured by utilizing the digital microscope VHX-600ESO. Based on analyzing the cavitation mechanism of the grinding fluid, the bubble dynamics model under conventional honing (CH) and ultrasonic vibration honing (UVH) was established respectively. Difference of dynamic behaviors of the bubble between the cases in UVH and CH was compared respectively, and the effects of acoustic amplitude and ultrasonic frequency on the bubble dynamics were simulated numerically using the Runge-Kutta fourth order method with variable step size adaptive control. Finally, the cavitation intensity of grinding fluids under ultrasound was measured quantitatively using acoustimeter. The results showed that the grinding fluid subjected to ultrasound can generate many bubbles and further forms numerous groups of araneose cavitation bubbles on the surface of the honing stone. The oscillation of the bubble under UVH is more intense than the case under CH, and the maximum velocity of the bubble wall under UVH is higher two magnitudes than the case under CH. For lower acoustic amplitude, the dynamic behaviors of the bubble under UVH are similar to that case under CH. As increasing acoustic amplitude, the cavitation intensity of the bubble is growing increased. Honing pressure has an inhabitation effect on cavitation effect of the grinding fluid. The perfect performance of cavitation of the grinding fluid can be obtained when the device of UVH is in the resonance. However, the cavitation intensity of the grinding fluid can be growing weakened with increasing ultrasonic frequency, when the device of UVH is in the off-resonance. The experimental results agree with the theoretical and numerical analysis, which provides a method for exploring applications of the cavitation effect in ultrasonic assisted machining. Copyright © 2017 Elsevier B.V. All rights reserved.

2D modeling of direct laser metal deposition process using a finite particle method

NASA Astrophysics Data System (ADS)

Anedaf, T.; Abbès, B.; Abbès, F.; Li, Y. M.

2018-05-01

Direct laser metal deposition is one of the material additive manufacturing processes used to produce complex metallic parts. A thorough understanding of the underlying physical phenomena is required to obtain a high-quality parts. In this work, a mathematical model is presented to simulate the coaxial laser direct deposition process tacking into account of mass addition, heat transfer, and fluid flow with free surface and melting. The fluid flow in the melt pool together with mass and energy balances are solved using the Computational Fluid Dynamics (CFD) software NOGRID-points, based on the meshless Finite Pointset Method (FPM). The basis of the computations is a point cloud, which represents the continuum fluid domain. Each finite point carries all fluid information (density, velocity, pressure and temperature). The dynamic shape of the molten zone is explicitly described by the point cloud. The proposed model is used to simulate a single layer cladding.

The Repeated Replacement Method: A Pure Lagrangian Meshfree Method for Computational Fluid Dynamics

PubMed Central

Walker, Wade A.

2012-01-01

In this paper we describe the repeated replacement method (RRM), a new meshfree method for computational fluid dynamics (CFD). RRM simulates fluid flow by modeling compressible fluids’ tendency to evolve towards a state of constant density, velocity, and pressure. To evolve a fluid flow simulation forward in time, RRM repeatedly “chops out” fluid from active areas and replaces it with new “flattened” fluid cells with the same mass, momentum, and energy. We call the new cells “flattened” because we give them constant density, velocity, and pressure, even though the chopped-out fluid may have had gradients in these primitive variables. RRM adaptively chooses the sizes and locations of the areas it chops out and replaces. It creates more and smaller new cells in areas of high gradient, and fewer and larger new cells in areas of lower gradient. This naturally leads to an adaptive level of accuracy, where more computational effort is spent on active areas of the fluid, and less effort is spent on inactive areas. We show that for common test problems, RRM produces results similar to other high-resolution CFD methods, while using a very different mathematical framework. RRM does not use Riemann solvers, flux or slope limiters, a mesh, or a stencil, and it operates in a purely Lagrangian mode. RRM also does not evaluate numerical derivatives, does not integrate equations of motion, and does not solve systems of equations. PMID:22866175

Active learning of constitutive relation from mesoscopic dynamics for macroscopic modeling of non-Newtonian flows

NASA Astrophysics Data System (ADS)

Zhao, Lifei; Li, Zhen; Caswell, Bruce; Ouyang, Jie; Karniadakis, George Em

2018-06-01

We simulate complex fluids by means of an on-the-fly coupling of the bulk rheology to the underlying microstructure dynamics. In particular, a continuum model of polymeric fluids is constructed without a pre-specified constitutive relation, but instead it is actively learned from mesoscopic simulations where the dynamics of polymer chains is explicitly computed. To couple the bulk rheology of polymeric fluids and the microscale dynamics of polymer chains, the continuum approach (based on the finite volume method) provides the transient flow field as inputs for the (mesoscopic) dissipative particle dynamics (DPD), and in turn DPD returns an effective constitutive relation to close the continuum equations. In this multiscale modeling procedure, we employ an active learning strategy based on Gaussian process regression (GPR) to minimize the number of expensive DPD simulations, where adaptively selected DPD simulations are performed only as necessary. Numerical experiments are carried out for flow past a circular cylinder of a non-Newtonian fluid, modeled at the mesoscopic level by bead-spring chains. The results show that only five DPD simulations are required to achieve an effective closure of the continuum equations at Reynolds number Re = 10. Furthermore, when Re is increased to 100, only one additional DPD simulation is required for constructing an extended GPR-informed model closure. Compared to traditional message-passing multiscale approaches, applying an active learning scheme to multiscale modeling of non-Newtonian fluids can significantly increase the computational efficiency. Although the method demonstrated here obtains only a local viscosity from the polymer dynamics, it can be extended to other multiscale models of complex fluids whose macro-rheology is unknown.

Three-dimensional Computational Fluid Dynamics Investigation of a Spinning Helicopter Slung Load

NASA Technical Reports Server (NTRS)

Theorn, J. N.; Duque, E. P. N.; Cicolani, L.; Halsey, R.

2005-01-01

After performing steady-state Computational Fluid Dynamics (CFD) calculations using OVERFLOW to validate the CFD method against static wind-tunnel data of a box-shaped cargo container, the same setup was used to investigate unsteady flow with a moving body. Results were compared to flight test data previously collected in which the container is spinning.

Dissolution Dynamic Nuclear Polarization capability study with fluid path

NASA Astrophysics Data System (ADS)

Malinowski, Ronja M.; Lipsø, Kasper W.; Lerche, Mathilde H.; Ardenkjær-Larsen, Jan H.

2016-11-01

Signal enhancement by hyperpolarization is a way of overcoming the low sensitivity in magnetic resonance; MRI in particular. One of the most well-known methods, dissolution Dynamic Nuclear Polarization, has been used clinically in cancer patients. One way of ensuring a low bioburden of the hyperpolarized product is by use of a closed fluid path that constitutes a barrier to contamination. The fluid path can be filled with the pharmaceuticals, i.e. imaging agent and solvents, in a clean room, and then stored or immediately used at the polarizer. In this study, we present a method of filling the fluid path that allows it to be reused. The filling method has been investigated in terms of reproducibility at two extrema, high dose for patient use and low dose for rodent studies, using [1-13C]pyruvate as example. We demonstrate that the filling method allows high reproducibility of six quality control parameters with standard deviations 3-10 times smaller than the acceptance criteria intervals in clinical studies.

Prospects and expectations for unstructured methods

NASA Technical Reports Server (NTRS)

Baker, Timothy J.

1995-01-01

The last decade has witnessed a vigorous and sustained research effort on unstructured methods for computational fluid dynamics. Unstructured mesh generators and flow solvers have evolved to the point where they are now in use for design purposes throughout the aerospace industry. In this paper we survey the various mesh types, structured as well as unstructured, and examine their relative strengths and weaknesses. We argue that unstructured methodology does offer the best prospect for the next generation of computational fluid dynamics algorithms.

High-Order Methods for Computational Fluid Dynamics: A Brief Review of Compact Differential Formulations on Unstructured Grids

NASA Technical Reports Server (NTRS)

Huynh, H. T.; Wang, Z. J.; Vincent, P. E.

2013-01-01

Popular high-order schemes with compact stencils for Computational Fluid Dynamics (CFD) include Discontinuous Galerkin (DG), Spectral Difference (SD), and Spectral Volume (SV) methods. The recently proposed Flux Reconstruction (FR) approach or Correction Procedure using Reconstruction (CPR) is based on a differential formulation and provides a unifying framework for these high-order schemes. Here we present a brief review of recent developments for the FR/CPR schemes as well as some pacing items.

On the micromechanics of slip events in sheared, fluid-saturated fault gouge

NASA Astrophysics Data System (ADS)

Dorostkar, Omid; Guyer, Robert A.; Johnson, Paul A.; Marone, Chris; Carmeliet, Jan

2017-06-01

We used a three-dimensional discrete element method coupled with computational fluid dynamics to study the poromechanical properties of dry and fluid-saturated granular fault gouge. The granular layer was sheared under dry conditions to establish a steady state condition of stick-slip dynamic failure, and then fluid was introduced to study its effect on subsequent failure events. The fluid-saturated case showed increased stick-slip recurrence time and larger slip events compared to the dry case. Particle motion induces fluid flow with local pressure variation, which in turn leads to high particle kinetic energy during slip due to increased drag forces from fluid on particles. The presence of fluid during the stick phase of loading promotes a more stable configuration evidenced by higher particle coordination number. Our coupled fluid-particle simulations provide grain-scale information that improves understanding of slip instabilities and illuminates details of phenomenological, macroscale observations.

Comments on Frequency Swept Rotating Input Perturbation Techniques and Identification of the Fluid Force Models in Rotor/bearing/seal Systems and Fluid Handling Machines

NASA Technical Reports Server (NTRS)

Muszynska, Agnes; Bently, Donald E.

1991-01-01

Perturbation techniques used for identification of rotating system dynamic characteristics are described. A comparison between two periodic frequency-swept perturbation methods applied in identification of fluid forces of rotating machines is presented. The description of the fluid force model identified by inputting circular periodic frequency-swept force is given. This model is based on the existence and strength of the circumferential flow, most often generated by the shaft rotation. The application of the fluid force model in rotor dynamic analysis is presented. It is shown that the rotor stability is an entire rotating system property. Some areas for further research are discussed.

Dynamic response of fluid inside a penny shaped crack

DOE Office of Scientific and Technical Information (OSTI.GOV)

Hayashi, Kazuo; Seki, Hitoshi

1997-12-31

In order to discuss the method for estimating the geometric characteristics of geothermal reservoir cracks, a theoretical study is performed on the dynamic response of the fluid inside a reservoir crack in a rock mass subjected to a dynamic excitation due to propagation of an elastic wave. As representative models of reservoir cracks, a penny shaped crack and a two-dimensional crack which are connected to a borehole are considered. It is found that the resonance frequency of the fluid motion is dependent on the crack size, the fluid`s viscosity and the permeability of the formation. The intensity of the resonancemore » is dependent on the fluid`s viscosity when the size, the aperture and the permeability are fixed. It is also found that, at a value of the fluid`s viscosity, the resonance of fluid pressure becomes strongest. The optimum value of the fluid`s viscosity is found to be almost perfectly determined by the permeability of the formation. Furthermore, it is revealed that, if the fluid`s viscosity is fixed to be the optimum value, the resonance frequency is almost independent of the permeability and aperture, but is dependent on the size of crack. Inversely speaking, this implies that the size of the reservoir crack can be estimated from the resonance frequency, if the fluid with the above mentioned optimum value of viscosity is employed for hydraulic fracturing.« less

A parallel direct-forcing fictitious domain method for simulating microswimmers

NASA Astrophysics Data System (ADS)

Gao, Tong; Lin, Zhaowu

2017-11-01

We present a 3D parallel direct-forcing fictitious domain method for simulating swimming micro-organisms at small Reynolds numbers. We treat the motile micro-swimmers as spherical rigid particles using the ``Squirmer'' model. The particle dynamics are solved on the moving Larangian meshes that overlay upon a fixed Eulerian mesh for solving the fluid motion, and the momentum exchange between the two phases is resolved by distributing pseudo body-forces over the particle interior regions which constrain the background fictitious fluids to follow the particle movement. While the solid and fluid subproblems are solved separately, no inner-iterations are required to enforce numerical convergence. We demonstrate the accuracy and robustness of the method by comparing our results with the existing analytical and numerical studies for various cases of single particle dynamics and particle-particle interactions. We also perform a series of numerical explorations to obtain statistical and rheological measurements to characterize the dynamics and structures of Squirmer suspensions. NSF DMS 1619960.

Potential applications of computational fluid dynamics to biofluid analysis

NASA Technical Reports Server (NTRS)

Kwak, D.; Chang, J. L. C.; Rogers, S. E.; Rosenfeld, M.; Kwak, D.

1988-01-01

Computational fluid dynamics was developed to the stage where it has become an indispensable part of aerospace research and design. In view of advances made in aerospace applications, the computational approach can be used for biofluid mechanics research. Several flow simulation methods developed for aerospace problems are briefly discussed for potential applications to biofluids, especially to blood flow analysis.

A phase-field method to analyze the dynamics of immiscible fluids in porous media

NASA Astrophysics Data System (ADS)

de Paoli, Marco; Roccon, Alessio; Zonta, Francesco; Soldati, Alfredo

2017-11-01

Liquid carbon dioxide (CO2) injected into geological formations (filled with brine) is not completely soluble in the surrounding fluid. For this reason, complex transport phenomena may occur across the interface that separates the two phases (CO2+brine and brine). Inspired by this geophysical instance, we used a Phase-Field Method (PFM) to describe the dynamics of two immiscible fluids in satured porous media. The basic idea of the PFM is to introduce an order parameter (ϕ) that varies continuously across the interfacial layer between the phases and is uniform in the bulk. The equation that describes the distribution of ϕ is the Cahn-Hilliard (CH) equation, which is coupled with the Darcy equation (to evaluate fluid velocity) through the buoyancy and Korteweg stress terms. The governing equations are solved through a pseudo-spectral technique (Fourier-Chebyshev). Our results show that the value of the surface tension between the two phases strongly influences the initial and the long term dynamics of the system. We believe that the proposed numerical approach, which grants an accurate evaluation of the interfacial fluxes of momentum/energy/species, is attractive to describe the transfer mechanism and the overall dynamics of immiscible and partially miscible phases.

Computational and analytical methods in nonlinear fluid dynamics

NASA Astrophysics Data System (ADS)

Walker, James

1993-09-01

The central focus of the program was on the application and development of modern analytical and computational methods to the solution of nonlinear problems in fluid dynamics and reactive gas dynamics. The research was carried out within the Division of Engineering Mathematics in the Department of Mechanical Engineering and Mechanics and principally involved Professors P.A. Blythe, E. Varley and J.D.A. Walker. In addition. the program involved various international collaborations. Professor Blythe completed work on reactive gas dynamics with Professor D. Crighton FRS of Cambridge University in the United Kingdom. Professor Walker and his students carried out joint work with Professor F.T. Smith, of University College London, on various problems in unsteady flow and turbulent boundary layers.

Comparison of two methods to determine fan performance curves using computational fluid dynamics

NASA Astrophysics Data System (ADS)

Onma, Patinya; Chantrasmi, Tonkid

2018-01-01

This work investigates a systematic numerical approach that employs Computational Fluid Dynamics (CFD) to obtain performance curves of a backward-curved centrifugal fan. Generating the performance curves requires a number of three-dimensional simulations with varying system loads at a fixed rotational speed. Two methods were used and their results compared to experimental data. The first method incrementally changes the mass flow late through the inlet boundary condition while the second method utilizes a series of meshes representing the physical damper blade at various angles. The generated performance curves from both methods are compared with an experiment setup in accordance with the AMCA fan performance testing standard.

Overview of Sensitivity Analysis and Shape Optimization for Complex Aerodynamic Configurations

NASA Technical Reports Server (NTRS)

Newman, Perry A.; Newman, James C., III; Barnwell, Richard W.; Taylor, Arthur C., III; Hou, Gene J.-W.

1998-01-01

This paper presents a brief overview of some of the more recent advances in steady aerodynamic shape-design sensitivity analysis and optimization, based on advanced computational fluid dynamics. The focus here is on those methods particularly well- suited to the study of geometrically complex configurations and their potentially complex associated flow physics. When nonlinear state equations are considered in the optimization process, difficulties are found in the application of sensitivity analysis. Some techniques for circumventing such difficulties are currently being explored and are included here. Attention is directed to methods that utilize automatic differentiation to obtain aerodynamic sensitivity derivatives for both complex configurations and complex flow physics. Various examples of shape-design sensitivity analysis for unstructured-grid computational fluid dynamics algorithms are demonstrated for different formulations of the sensitivity equations. Finally, the use of advanced, unstructured-grid computational fluid dynamics in multidisciplinary analyses and multidisciplinary sensitivity analyses within future optimization processes is recommended and encouraged.

A generalized transport-velocity formulation for smoothed particle hydrodynamics

DOE Office of Scientific and Technical Information (OSTI.GOV)

Zhang, Chi; Hu, Xiangyu Y., E-mail: xiangyu.hu@tum.de; Adams, Nikolaus A.

The standard smoothed particle hydrodynamics (SPH) method suffers from tensile instability. In fluid-dynamics simulations this instability leads to particle clumping and void regions when negative pressure occurs. In solid-dynamics simulations, it results in unphysical structure fragmentation. In this work the transport-velocity formulation of Adami et al. (2013) is generalized for providing a solution of this long-standing problem. Other than imposing a global background pressure, a variable background pressure is used to modify the particle transport velocity and eliminate the tensile instability completely. Furthermore, such a modification is localized by defining a shortened smoothing length. The generalized formulation is suitable formore » fluid and solid materials with and without free surfaces. The results of extensive numerical tests on both fluid and solid dynamics problems indicate that the new method provides a unified approach for multi-physics SPH simulations.« less

Experimental investigation of the flow dynamics and rheology of complex fluids in pipe flow by hybrid multi-scale velocimetry

NASA Astrophysics Data System (ADS)

Haavisto, Sanna; Cardona, Maria J.; Salmela, Juha; Powell, Robert L.; McCarthy, Michael J.; Kataja, Markku; Koponen, Antti I.

2017-11-01

A hybrid multi-scale velocimetry method utilizing Doppler optical coherence tomography in combination with either magnetic resonance imaging or ultrasound velocity profiling is used to investigate pipe flow of four rheologically different working fluids under varying flow regimes. These fluids include water, an aqueous xanthan gum solution, a softwood fiber suspension, and a microfibrillated cellulose suspension. The measurement setup enables not only the analysis of the rheological (bulk) behavior of a studied fluid but gives simultaneously information on their wall layer dynamics, both of which are needed for analyzing and solving practical fluid flow-related problems. Preliminary novel results on rheological and boundary layer flow properties of the working fluids are reported and the potential of the hybrid measurement setup is demonstrated.

Simulation of plume dynamics by the Lattice Boltzmann Method

NASA Astrophysics Data System (ADS)

Mora, Peter; Yuen, David A.

2017-09-01

The Lattice Boltzmann Method (LBM) is a semi-microscopic method to simulate fluid mechanics by modelling distributions of particles moving and colliding on a lattice. We present 2-D simulations using the LBM of a fluid in a rectangular box being heated from below, and cooled from above, with a Rayleigh of Ra = 108, similar to current estimates of the Earth's mantle, and a Prandtl number of 5000. At this Prandtl number, the flow is found to be in the non-inertial regime where the inertial terms denoted I ≪ 1. Hence, the simulations presented lie within the regime of relevance for geodynamical problems. We obtain narrow upwelling plumes with mushroom heads and chutes of downwelling fluid as expected of a flow in the non-inertial regime. The method developed demonstrates that the LBM has great potential for simulating thermal convection and plume dynamics relevant to geodynamics, albeit with some limitations.

A matrix-free implicit unstructured multigrid finite volume method for simulating structural dynamics and fluid structure interaction

NASA Astrophysics Data System (ADS)

Lv, X.; Zhao, Y.; Huang, X. Y.; Xia, G. H.; Su, X. H.

2007-07-01

A new three-dimensional (3D) matrix-free implicit unstructured multigrid finite volume (FV) solver for structural dynamics is presented in this paper. The solver is first validated using classical 2D and 3D cantilever problems. It is shown that very accurate predictions of the fundamental natural frequencies of the problems can be obtained by the solver with fast convergence rates. This method has been integrated into our existing FV compressible solver [X. Lv, Y. Zhao, et al., An efficient parallel/unstructured-multigrid preconditioned implicit method for simulating 3d unsteady compressible flows with moving objects, Journal of Computational Physics 215(2) (2006) 661-690] based on the immersed membrane method (IMM) [X. Lv, Y. Zhao, et al., as mentioned above]. Results for the interaction between the fluid and an immersed fixed-free cantilever are also presented to demonstrate the potential of this integrated fluid-structure interaction approach.

Data-Driven Model Reduction and Transfer Operator Approximation

NASA Astrophysics Data System (ADS)

Klus, Stefan; Nüske, Feliks; Koltai, Péter; Wu, Hao; Kevrekidis, Ioannis; Schütte, Christof; Noé, Frank

2018-06-01

In this review paper, we will present different data-driven dimension reduction techniques for dynamical systems that are based on transfer operator theory as well as methods to approximate transfer operators and their eigenvalues, eigenfunctions, and eigenmodes. The goal is to point out similarities and differences between methods developed independently by the dynamical systems, fluid dynamics, and molecular dynamics communities such as time-lagged independent component analysis, dynamic mode decomposition, and their respective generalizations. As a result, extensions and best practices developed for one particular method can be carried over to other related methods.

Three-dimensional simulation of pseudopod-driven swimming of amoeboid cells

NASA Astrophysics Data System (ADS)

Campbell, Eric; Bagchi, Prosenjit

2016-11-01

Pseudopod-driven locomotion is common in eukaryotic cells, such as amoeba, neutrophils, and cancer cells. Pseudopods are protrusions of the cell body that grow, bifurcate, and retract. Due to the dynamic nature of pseudopods, the shape of a motile cell constantly changes. The actin-myosin protein dynamics is a likely mechanism for pseudopod growth. Existing theoretical models often focus on the acto-myosin dynamics, and not the whole cell shape dynamics. Here we present a full 3D simulation of pseudopod-driven motility by coupling a surface-bound reaction-diffusion (RD) model for the acto-myosin dynamics, a continuum model for the cell membrane deformation, and flow of the cytoplasmic and extracellular fluids. The whole cell is represented as a viscous fluid surrounded by a membrane. A finite-element method is used to solve the membrane deformation, and the RD model on the deforming membrane, while a finite-difference/spectral method is used to solve the flow fields inside and outside the cell. The fluid flow and cell deformation are coupled by the immersed-boundary method. The model predicts pseudopod growth, bifurcation, and retraction as observed for a swimming amoeba. The work provides insights on the role of membrane stiffness and cytoplasmic viscosity on amoeboid swimming. Funded by NSF CBET 1438255.

Reference Computational Meshing Strategy for Computational Fluid Dynamics Simulation of Departure from Nucleate BoilingReference Computational Meshing Strategy for Computational Fluid Dynamics Simulation of Departure from Nucleate Boiling

DOE Office of Scientific and Technical Information (OSTI.GOV)

Pointer, William David

The objective of this effort is to establish a strategy and process for generation of suitable computational mesh for computational fluid dynamics simulations of departure from nucleate boiling in a 5 by 5 fuel rod assembly held in place by PWR mixing vane spacer grids. This mesh generation process will support ongoing efforts to develop, demonstrate and validate advanced multi-phase computational fluid dynamics methods that enable more robust identification of dryout conditions and DNB occurrence.Building upon prior efforts and experience, multiple computational meshes were developed using the native mesh generation capabilities of the commercial CFD code STAR-CCM+. These meshes weremore » used to simulate two test cases from the Westinghouse 5 by 5 rod bundle facility. The sensitivity of predicted quantities of interest to the mesh resolution was then established using two evaluation methods, the Grid Convergence Index method and the Least Squares method. This evaluation suggests that the Least Squares method can reliably establish the uncertainty associated with local parameters such as vector velocity components at a point in the domain or surface averaged quantities such as outlet velocity magnitude. However, neither method is suitable for characterization of uncertainty in global extrema such as peak fuel surface temperature, primarily because such parameters are not necessarily associated with a fixed point in space. This shortcoming is significant because the current generation algorithm for identification of DNB event conditions relies on identification of such global extrema. Ongoing efforts to identify DNB based on local surface conditions will address this challenge« less

CFD Analysis of the SBXC Glider Airframe

DTIC Science & Technology

2016-06-01

mathematically on finite element methods. To validate and verify the methodology developed, a mathematical comparison was made with the previous research data...greater than 15 m/s. 14. SUBJECT TERMS finite element method, computational fluid dynamics, Y Plus, mesh element quality, aerodynamic data, fluid...based mathematically on finite element methods. To validate and verify the methodology developed, a mathematical comparison was made with the

Computational fluid dynamics for propulsion technology: Geometric grid visualization in CFD-based propulsion technology research

NASA Technical Reports Server (NTRS)

Ziebarth, John P.; Meyer, Doug

1992-01-01

The coordination is examined of necessary resources, facilities, and special personnel to provide technical integration activities in the area of computational fluid dynamics applied to propulsion technology. Involved is the coordination of CFD activities between government, industry, and universities. Current geometry modeling, grid generation, and graphical methods are established to use in the analysis of CFD design methodologies.

Parallel Three-Dimensional Computation of Fluid Dynamics and Fluid-Structure Interactions of Ram-Air Parachutes

NASA Technical Reports Server (NTRS)

Tezduyar, Tayfun E.

1998-01-01

This is a final report as far as our work at University of Minnesota is concerned. The report describes our research progress and accomplishments in development of high performance computing methods and tools for 3D finite element computation of aerodynamic characteristics and fluid-structure interactions (FSI) arising in airdrop systems, namely ram-air parachutes and round parachutes. This class of simulations involves complex geometries, flexible structural components, deforming fluid domains, and unsteady flow patterns. The key components of our simulation toolkit are a stabilized finite element flow solver, a nonlinear structural dynamics solver, an automatic mesh moving scheme, and an interface between the fluid and structural solvers; all of these have been developed within a parallel message-passing paradigm.

Model structure identification for wastewater treatment simulation based on computational fluid dynamics.

PubMed

Alex, J; Kolisch, G; Krause, K

2002-01-01

The objective of this presented project is to use the results of an CFD simulation to automatically, systematically and reliably generate an appropriate model structure for simulation of the biological processes using CSTR activated sludge compartments. Models and dynamic simulation have become important tools for research but also increasingly for the design and optimisation of wastewater treatment plants. Besides the biological models several cases are reported about the application of computational fluid dynamics ICFD) to wastewater treatment plants. One aim of the presented method to derive model structures from CFD results is to exclude the influence of empirical structure selection to the result of dynamic simulations studies of WWTPs. The second application of the approach developed is the analysis of badly performing treatment plants where the suspicion arises that bad flow behaviour such as short cut flows is part of the problem. The method suggested requires as the first step the calculation of fluid dynamics of the biological treatment step at different loading situations by use of 3-dimensional CFD simulation. The result of this information is used to generate a suitable model structure for conventional dynamic simulation of the treatment plant by use of a number of CSTR modules with a pattern of exchange flows between the tanks automatically. The method is explained in detail and the application to the WWTP Wuppertal Buchenhofen is presented.

Analysis of nanoscale two-phase flow of argon using molecular dynamics

NASA Astrophysics Data System (ADS)

Verma, Abhishek Kumar; Kumar, Rakesh

2014-12-01

Two phase flows through micro and nanochannels have attracted a lot of attention because of their immense applicability to many advanced fields such as MEMS/NEMS, electronic cooling, bioengineering etc. In this work, a molecular dynamics simulation method is employed to study the condensation process of superheated argon vapor force driven flow through a nanochannel combining fluid flow and heat transfer. A simple and effective particle insertion method is proposed to model phase change of argon based on non-periodic boundary conditions in the simulation domain. Starting from a crystalline solid wall of channel, the condensation process evolves from a transient unsteady state where we study the influence of different wall temperatures and fluid wall interactions on interfacial and heat transport properties of two phase flows. Subsequently, we analyzed transient temperature, density and velocity fields across the channel and their dependency on varying wall temperature and fluid wall interaction, after a dynamic equilibrium is achieved in phase transition. Quasi-steady nonequilibrium temperature profile, heat flux and interfacial thermal resistance were analyzed. The results demonstrate that the molecular dynamics method, with the proposed particle insertion method, effectively solves unsteady nonequilibrium two phase flows at nanoscale resolutions whose interphase between liquid and vapor phase is typically of the order of a few molecular diameters.

A Versatile Nonlinear Method for Predictive Modeling

NASA Technical Reports Server (NTRS)

Liou, Meng-Sing; Yao, Weigang

2015-01-01

As computational fluid dynamics techniques and tools become widely accepted for realworld practice today, it is intriguing to ask: what areas can it be utilized to its potential in the future. Some promising areas include design optimization and exploration of fluid dynamics phenomena (the concept of numerical wind tunnel), in which both have the common feature where some parameters are varied repeatedly and the computation can be costly. We are especially interested in the need for an accurate and efficient approach for handling these applications: (1) capturing complex nonlinear dynamics inherent in a system under consideration and (2) versatility (robustness) to encompass a range of parametric variations. In our previous paper, we proposed to use first-order Taylor expansion collected at numerous sampling points along a trajectory and assembled together via nonlinear weighting functions. The validity and performance of this approach was demonstrated for a number of problems with a vastly different input functions. In this study, we are especially interested in enhancing the method's accuracy; we extend it to include the second-orer Taylor expansion, which however requires a complicated evaluation of Hessian matrices for a system of equations, like in fluid dynamics. We propose a method to avoid these Hessian matrices, while maintaining the accuracy. Results based on the method are presented to confirm its validity.

Deterministic and stochastic algorithms for resolving the flow fields in ducts and networks using energy minimization

NASA Astrophysics Data System (ADS)

Sochi, Taha

2016-09-01

Several deterministic and stochastic multi-variable global optimization algorithms (Conjugate Gradient, Nelder-Mead, Quasi-Newton and global) are investigated in conjunction with energy minimization principle to resolve the pressure and volumetric flow rate fields in single ducts and networks of interconnected ducts. The algorithms are tested with seven types of fluid: Newtonian, power law, Bingham, Herschel-Bulkley, Ellis, Ree-Eyring and Casson. The results obtained from all those algorithms for all these types of fluid agree very well with the analytically derived solutions as obtained from the traditional methods which are based on the conservation principles and fluid constitutive relations. The results confirm and generalize the findings of our previous investigations that the energy minimization principle is at the heart of the flow dynamics systems. The investigation also enriches the methods of computational fluid dynamics for solving the flow fields in tubes and networks for various types of Newtonian and non-Newtonian fluids.

Dissolution Dynamic Nuclear Polarization capability study with fluid path.

PubMed

Malinowski, Ronja M; Lipsø, Kasper W; Lerche, Mathilde H; Ardenkjær-Larsen, Jan H

2016-11-01

Signal enhancement by hyperpolarization is a way of overcoming the low sensitivity in magnetic resonance; MRI in particular. One of the most well-known methods, dissolution Dynamic Nuclear Polarization, has been used clinically in cancer patients. One way of ensuring a low bioburden of the hyperpolarized product is by use of a closed fluid path that constitutes a barrier to contamination. The fluid path can be filled with the pharmaceuticals, i.e. imaging agent and solvents, in a clean room, and then stored or immediately used at the polarizer. In this study, we present a method of filling the fluid path that allows it to be reused. The filling method has been investigated in terms of reproducibility at two extrema, high dose for patient use and low dose for rodent studies, using [1-13C]pyruvate as example. We demonstrate that the filling method allows high reproducibility of six quality control parameters with standard deviations 3-10 times smaller than the acceptance criteria intervals in clinical studies. Copyright © 2016 Elsevier Inc. All rights reserved.

Nanoscale hydrodynamics near solids

NASA Astrophysics Data System (ADS)

Camargo, Diego; de la Torre, J. A.; Duque-Zumajo, D.; Español, Pep; Delgado-Buscalioni, Rafael; Chejne, Farid

2018-02-01

Density Functional Theory (DFT) is a successful and well-established theory for the study of the structure of simple and complex fluids at equilibrium. The theory has been generalized to dynamical situations when the underlying dynamics is diffusive as in, for example, colloidal systems. However, there is no such a clear foundation for Dynamic DFT (DDFT) for the case of simple fluids in contact with solid walls. In this work, we derive DDFT for simple fluids by including not only the mass density field but also the momentum density field of the fluid. The standard projection operator method based on the Kawasaki-Gunton operator is used for deriving the equations for the average value of these fields. The solid is described as featureless under the assumption that all the internal degrees of freedom of the solid relax much faster than those of the fluid (solid elasticity is irrelevant). The fluid moves according to a set of non-local hydrodynamic equations that include explicitly the forces due to the solid. These forces are of two types, reversible forces emerging from the free energy density functional, and accounting for impenetrability of the solid, and irreversible forces that involve the velocity of both the fluid and the solid. These forces are localized in the vicinity of the solid surface. The resulting hydrodynamic equations should allow one to study dynamical regimes of simple fluids in contact with solid objects in isothermal situations.

Physically-Based Modelling and Real-Time Simulation of Fluids.

NASA Astrophysics Data System (ADS)

Chen, Jim Xiong

1995-01-01

Simulating physically realistic complex fluid behaviors presents an extremely challenging problem for computer graphics researchers. Such behaviors include the effects of driving boats through water, blending differently colored fluids, rain falling and flowing on a terrain, fluids interacting in a Distributed Interactive Simulation (DIS), etc. Such capabilities are useful in computer art, advertising, education, entertainment, and training. We present a new method for physically-based modeling and real-time simulation of fluids in computer graphics and dynamic virtual environments. By solving the 2D Navier -Stokes equations using a CFD method, we map the surface into 3D using the corresponding pressures in the fluid flow field. This achieves realistic real-time fluid surface behaviors by employing the physical governing laws of fluids but avoiding extensive 3D fluid dynamics computations. To complement the surface behaviors, we calculate fluid volume and external boundary changes separately to achieve full 3D general fluid flow. To simulate physical activities in a DIS, we introduce a mechanism which uses a uniform time scale proportional to the clock-time and variable time-slicing to synchronize physical models such as fluids in the networked environment. Our approach can simulate many different fluid behaviors by changing the internal or external boundary conditions. It can model different kinds of fluids by varying the Reynolds number. It can simulate objects moving or floating in fluids. It can also produce synchronized general fluid flows in a DIS. Our model can serve as a testbed to simulate many other fluid phenomena which have never been successfully modeled previously.

Data-driven sensor placement from coherent fluid structures

NASA Astrophysics Data System (ADS)

Manohar, Krithika; Kaiser, Eurika; Brunton, Bingni W.; Kutz, J. Nathan; Brunton, Steven L.

2017-11-01

Optimal sensor placement is a central challenge in the prediction, estimation and control of fluid flows. We reinterpret sensor placement as optimizing discrete samples of coherent fluid structures for full state reconstruction. This permits a drastic reduction in the number of sensors required for faithful reconstruction, since complex fluid interactions can often be described by a small number of coherent structures. Our work optimizes point sensors using the pivoted matrix QR factorization to sample coherent structures directly computed from flow data. We apply this sampling technique in conjunction with various data-driven modal identification methods, including the proper orthogonal decomposition (POD) and dynamic mode decomposition (DMD). In contrast to POD-based sensors, DMD demonstrably enables the optimization of sensors for prediction in systems exhibiting multiple scales of dynamics. Finally, reconstruction accuracy from pivot sensors is shown to be competitive with sensors obtained using traditional computationally prohibitive optimization methods.

Multiphase Fluid Dynamics for Spacecraft Applications

NASA Astrophysics Data System (ADS)

Shyy, W.; Sim, J.

2011-09-01

Multiphase flows involving moving interfaces between different fluids/phases are observed in nature as well as in a wide range of engineering applications. With the recent development of high fidelity computational techniques, a number of challenging multiphase flow problems can now be computed. We introduce the basic notion of the main categories of multiphase flow computation; Lagrangian, Eulerian, and Eulerian-Lagrangian techniques to represent and follow interface, and sharp and continuous interface methods to model interfacial dynamics. The marker-based adaptive Eulerian-Lagrangian method, which is one of the most popular methods, is highlighted with microgravity and space applications including droplet collision and spacecraft liquid fuel tank surface stability.

Study of the long-time dynamics of a viscous vortex sheet with a fully adaptive nonstiff method

NASA Astrophysics Data System (ADS)

Ceniceros, Hector D.; Roma, Alexandre M.

2004-12-01

A numerical investigation of the long-time dynamics of two immiscible two-dimensional fluids shearing past one another is presented. The fluids are incompressible and the interface between the bulk phases is subjected to surface tension. The simple case of density and viscosity matched fluids is considered. The two-dimensional Navier-Stokes equations are solved numerically with a fully adaptive nonstiff strategy based on the immersed boundary method. Dynamically adaptive mesh refinements are used to cover at all times the separately tracked fluid interface at the finest grid level. In addition, by combining adaptive front tracking, in the form of continuous interface marker equidistribution, with a predictor-corrector discretization an efficient method is introduced to successfully treat the well-known numerical difficulties associated with surface tension. The resulting numerical method can be used to compute stably and with high resolution the flow for wide-ranging Weber numbers but this study focuses on the computationally challenging cases for which elongated fingering and interface roll-up are observed. To assess the importance of the viscous and vortical effects in the interfacial dynamics the full viscous flow simulations are compared with inviscid counterparts computed with a state-of-the-art boundary integral method. In the examined cases of roll-up, it is found that in contrast to the inviscid flow in which the interface undergoes a topological reconfiguration, the viscous interface remarkably escapes self-intersection and rich long-time dynamics due to separation, transport, and diffusion of vorticity is observed. An even more striking motion occurs at an intermediate Weber number for which elongated interpenetrating fingers of fluid develop. In this case, it is found that the Kelvin-Helmholtz instability weakens due to shedding of vorticity and unlike the inviscid counterpart in which there is indefinite finger growth the viscous interface is pulled back by surface tension. As the interface recedes, thin necks connecting pockets of fluid with the rest of the fingers form. Narrow jets are observed at the necking regions but the vorticity there ultimately appears to be insufficient to drain all the fluid and cause reconnection. However, at another point, two disparate portions of the interface come in close proximity as the interface continues to contract. Large curvature points and an intense concentration of vorticity are observed in this region and then the motion is abruptly terminated by the collapse of the interface.

Immersed Boundary Simulations of Active Fluid Droplets

PubMed Central

Hawkins, Rhoda J.

2016-01-01

We present numerical simulations of active fluid droplets immersed in an external fluid in 2-dimensions using an Immersed Boundary method to simulate the fluid droplet interface as a Lagrangian mesh. We present results from two example systems, firstly an active isotropic fluid boundary consisting of particles that can bind and unbind from the interface and generate surface tension gradients through active contractility. Secondly, a droplet filled with an active polar fluid with homeotropic anchoring at the droplet interface. These two systems demonstrate spontaneous symmetry breaking and steady state dynamics resembling cell motility and division and show complex feedback mechanisms with minimal degrees of freedom. The simulations outlined here will be useful for quantifying the wide range of dynamics observable in these active systems and modelling the effects of confinement in a consistent and adaptable way. PMID:27606609

On the role of fluids in stick-slip dynamics of saturated granular fault gouge using a coupled computational fluid dynamics-discrete element approach

NASA Astrophysics Data System (ADS)

Dorostkar, Omid; Guyer, Robert A.; Johnson, Paul A.; Marone, Chris; Carmeliet, Jan

2017-05-01

The presence of fault gouge has considerable influence on slip properties of tectonic faults and the physics of earthquake rupture. The presence of fluids within faults also plays a significant role in faulting and earthquake processes. In this paper, we present 3-D discrete element simulations of dry and fluid-saturated granular fault gouge and analyze the effect of fluids on stick-slip behavior. Fluid flow is modeled using computational fluid dynamics based on the Navier-Stokes equations for an incompressible fluid and modified to take into account the presence of particles. Analysis of a long time train of slip events shows that the (1) drop in shear stress, (2) compaction of granular layer, and (3) the kinetic energy release during slip all increase in magnitude in the presence of an incompressible fluid, compared to dry conditions. We also observe that on average, the recurrence interval between slip events is longer for fluid-saturated granular fault gouge compared to the dry case. This observation is consistent with the occurrence of larger events in the presence of fluid. It is found that the increase in kinetic energy during slip events for saturated conditions can be attributed to the increased fluid flow during slip. Our observations emphasize the important role that fluid flow and fluid-particle interactions play in tectonic fault zones and show in particular how discrete element method (DEM) models can help understand the hydromechanical processes that dictate fault slip.

Grain scale observations of stick-slip dynamics in fluid saturated granular fault gouge

NASA Astrophysics Data System (ADS)

Johnson, P. A.; Dorostkar, O.; Guyer, R. A.; Marone, C.; Carmeliet, J.

2017-12-01

We are studying granular mechanics during slip. In the present work, we conduct coupled computational fluid dynamics (CFD) and discrete element method (DEM) simulations to study grain scale characteristics of slip instabilities in fluid saturated granular fault gouge. The granular sample is confined with constant normal load (10 MPa), and sheared with constant velocity (0.6 mm/s). This loading configuration is chosen to promote stick-slip dynamics, based on a phase-space study. Fluid is introduced in the beginning of stick phase and characteristics of slip events i.e. macroscopic friction coefficient, kinetic energy and layer thickness are monitored. At the grain scale, we monitor particle coordination number, fluid-particle interaction forces as well as particle and fluid kinetic energy. Our observations show that presence of fluids in a drained granular fault gouge stabilizes the layer in the stick phase and increases the recurrence time. In saturated model, we observe that average particle coordination number reaches higher values compared to dry granular gouge. Upon slip, we observe that a larger portion of the granular sample is mobilized in saturated gouge compared to dry system. We also observe that regions with high particle kinetic energy are correlated with zones of high fluid motion. Our observations highlight that spatiotemporal profile of fluid dynamic pressure affects the characteristics of slip instabilities, increasing macroscopic friction coefficient drop, kinetic energy release and granular layer compaction. We show that numerical simulations help characterize the micromechanics of fault mechanics.

Oxazine-based sensor for contaminant detection, fabrication method therefor, and uses thereof

DOEpatents

Nnanna, Agbai Agwu; Jalal, Ahmed Hasnian

2014-05-27

A sensor, a method for its fabrication, and a method for its use to detect contaminants, for example, ammonia, in stagnant and dynamic fluid media, especially liquid media. The sensor is an opto-chemical sensor that includes a polymer optical fiber, a sensing layer comprising oxazine 170 perchlorate on the polymer optical fiber, and a membrane layer on the sensing layer. The membrane layer is gas permeable and not permeable to the fluid in the fluid system, and moisture is entrapped by and between the sensing and membrane layers.

Evolution of Advection Upstream Splitting Method Schemes

NASA Technical Reports Server (NTRS)

Liou, Meng-Sing

2010-01-01

This paper focuses on the evolution of advection upstream splitting method(AUSM) schemes. The main ingredients that have led to the development of modern computational fluid dynamics (CFD) methods have been reviewed, thus the ideas behind AUSM. First and foremost is the concept of upwinding. Second, the use of Riemann problem in constructing the numerical flux in the finite-volume setting. Third, the necessity of including all physical processes, as characterised by the linear (convection) and nonlinear (acoustic) fields. Fourth, the realisation of separating the flux into convection and pressure fluxes. The rest of this review briefly outlines the technical evolution of AUSM and more details can be found in the cited references. Keywords: Computational fluid dynamics methods, hyperbolic systems, advection upstream splitting method, conservation laws, upwinding, CFD

Method for inverting reflection trace data from 3-D and 4-D seismic surveys and identifying subsurface fluid and pathways in and among hydrocarbon reservoirs based on impedance models

DOEpatents

He, W.; Anderson, R.N.

1998-08-25

A method is disclosed for inverting 3-D seismic reflection data obtained from seismic surveys to derive impedance models for a subsurface region, and for inversion of multiple 3-D seismic surveys (i.e., 4-D seismic surveys) of the same subsurface volume, separated in time to allow for dynamic fluid migration, such that small scale structure and regions of fluid and dynamic fluid flow within the subsurface volume being studied can be identified. The method allows for the mapping and quantification of available hydrocarbons within a reservoir and is thus useful for hydrocarbon prospecting and reservoir management. An iterative seismic inversion scheme constrained by actual well log data which uses a time/depth dependent seismic source function is employed to derive impedance models from 3-D and 4-D seismic datasets. The impedance values can be region grown to better isolate the low impedance hydrocarbon bearing regions. Impedance data derived from multiple 3-D seismic surveys of the same volume can be compared to identify regions of dynamic evolution and bypassed pay. Effective Oil Saturation or net oil thickness can also be derived from the impedance data and used for quantitative assessment of prospective drilling targets and reservoir management. 20 figs.

Method for inverting reflection trace data from 3-D and 4-D seismic surveys and identifying subsurface fluid and pathways in and among hydrocarbon reservoirs based on impedance models

DOEpatents

He, Wei; Anderson, Roger N.

1998-01-01

A method is disclosed for inverting 3-D seismic reflection data obtained from seismic surveys to derive impedance models for a subsurface region, and for inversion of multiple 3-D seismic surveys (i.e., 4-D seismic surveys) of the same subsurface volume, separated in time to allow for dynamic fluid migration, such that small scale structure and regions of fluid and dynamic fluid flow within the subsurface volume being studied can be identified. The method allows for the mapping and quantification of available hydrocarbons within a reservoir and is thus useful for hydrocarbon prospecting and reservoir management. An iterative seismic inversion scheme constrained by actual well log data which uses a time/depth dependent seismic source function is employed to derive impedance models from 3-D and 4-D seismic datasets. The impedance values can be region grown to better isolate the low impedance hydrocarbon bearing regions. Impedance data derived from multiple 3-D seismic surveys of the same volume can be compared to identify regions of dynamic evolution and bypassed pay. Effective Oil Saturation or net oil thickness can also be derived from the impedance data and used for quantitative assessment of prospective drilling targets and reservoir management.

Dynamic analysis of a 5-MW tripod offshore wind turbine by considering fluid-structure interaction

NASA Astrophysics Data System (ADS)

Zhang, Li-wei; Li, Xin

2017-10-01

Fixed offshore wind turbines usually have large underwater supporting structures. The fluid influences the dynamic characteristics of the structure system. The dynamic model of a 5-MW tripod offshore wind turbine considering the pile-soil system and fluid structure interaction (FSI) is established, and the structural modes in air and in water are obtained by use of ANSYS. By comparing low-order natural frequencies and mode shapes, the influence of sea water on the free vibration characteristics of offshore wind turbine is analyzed. On basis of the above work, seismic responses under excitation by El-Centro waves are calculated by the time-history analysis method. The results reveal that the dynamic responses such as the lateral displacement of the foundation and the section bending moment of the tubular piles increase substantially under the influence of the added-mass and hydrodynamic pressure of sea water. The method and conclusions presented in this paper can provide a theoretical reference for structure design and analysis of offshore wind turbines fixed in deep seawater.

Dynamics of a passive micro-vibration isolator based on a pretensioned plane cable net structure and fluid damper

NASA Astrophysics Data System (ADS)

Chen, Yanhao; Lu, Qi; Jing, Bo; Zhang, Zhiyi

2016-09-01

This paper addresses dynamic modelling and experiments on a passive vibration isolator for application in the space environment. The isolator is composed of a pretensioned plane cable net structure and a fluid damper in parallel. Firstly, the frequency response function (FRF) of a single cable is analysed according to the string theory, and the FRF synthesis method is adopted to establish a dynamic model of the plane cable net structure. Secondly, the equivalent damping coefficient of the fluid damper is analysed. Thirdly, experiments are carried out to compare the plane cable net structure, the fluid damper and the vibration isolator formed by the net and the damper, respectively. It is shown that the plane cable net structure can achieve substantial vibration attenuation but has a great amplification at its resonance frequency due to the light damping of cables. The damping effect of fluid damper is acceptable without taking the poor carrying capacity into consideration. Compared to the plane cable net structure and the fluid damper, the isolator has an acceptable resonance amplification as well as vibration attenuation.

Correlating contact line capillarity and dynamic contact angle hysteresis in surfactant-nanoparticle based complex fluids

NASA Astrophysics Data System (ADS)

Harikrishnan, A. R.; Dhar, Purbarun; Agnihotri, Prabhat K.; Gedupudi, Sateesh; Das, Sarit K.

2018-04-01

Dynamic wettability and contact angle hysteresis can be correlated to shed insight onto any solid-liquid interaction. Complex fluids are capable of altering the expected hysteresis and dynamic wetting behavior due to interfacial interactions. We report the effect of capillary number on the dynamic advancing and receding contact angles of surfactant-based nanocolloidal solutions on hydrophilic, near hydrophobic, and superhydrophobic surfaces by performing forced wetting and de-wetting experiments by employing the embedded needle method. A segregated study is performed to infer the contributing effects of the constituents and effects of particle morphology. The static contact angle hysteresis is found to be a function of particle and surfactant concentrations and greatly depends on the nature of the morphology of the particles. An order of estimate of line energy and a dynamic flow parameter called spreading factor and the transient variations of these parameters are explored which sheds light on the dynamics of contact line movement and response to perturbation of three-phase contact. The Cox-Voinov-Tanner law was found to hold for hydrophilic and a weak dependency on superhydrophobic surfaces with capillary number, and even for the complex fluids, with a varying degree of dependency for different fluids.

Communication: A method to compute the transport coefficient of pure fluids diffusing through planar interfaces from equilibrium molecular dynamics simulations.

PubMed

Vermorel, Romain; Oulebsir, Fouad; Galliero, Guillaume

2017-09-14

The computation of diffusion coefficients in molecular systems ranks among the most useful applications of equilibrium molecular dynamics simulations. However, when dealing with the problem of fluid diffusion through vanishingly thin interfaces, classical techniques are not applicable. This is because the volume of space in which molecules diffuse is ill-defined. In such conditions, non-equilibrium techniques allow for the computation of transport coefficients per unit interface width, but their weak point lies in their inability to isolate the contribution of the different physical mechanisms prone to impact the flux of permeating molecules. In this work, we propose a simple and accurate method to compute the diffusional transport coefficient of a pure fluid through a planar interface from equilibrium molecular dynamics simulations, in the form of a diffusion coefficient per unit interface width. In order to demonstrate its validity and accuracy, we apply our method to the case study of a dilute gas diffusing through a smoothly repulsive single-layer porous solid. We believe this complementary technique can benefit to the interpretation of the results obtained on single-layer membranes by means of complex non-equilibrium methods.

Methods for simulation-based analysis of fluid-structure interaction.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Barone, Matthew Franklin; Payne, Jeffrey L.

2005-10-01

Methods for analysis of fluid-structure interaction using high fidelity simulations are critically reviewed. First, a literature review of modern numerical techniques for simulation of aeroelastic phenomena is presented. The review focuses on methods contained within the arbitrary Lagrangian-Eulerian (ALE) framework for coupling computational fluid dynamics codes to computational structural mechanics codes. The review treats mesh movement algorithms, the role of the geometric conservation law, time advancement schemes, wetted surface interface strategies, and some representative applications. The complexity and computational expense of coupled Navier-Stokes/structural dynamics simulations points to the need for reduced order modeling to facilitate parametric analysis. The proper orthogonalmore » decomposition (POD)/Galerkin projection approach for building a reduced order model (ROM) is presented, along with ideas for extension of the methodology to allow construction of ROMs based on data generated from ALE simulations.« less

Computational fluid dynamics - The coming revolution

NASA Technical Reports Server (NTRS)

Graves, R. A., Jr.

1982-01-01

The development of aerodynamic theory is traced from the days of Aristotle to the present, with the next stage in computational fluid dynamics dependent on superspeed computers for flow calculations. Additional attention is given to the history of numerical methods inherent in writing computer codes applicable to viscous and inviscid analyses for complex configurations. The advent of the superconducting Josephson junction is noted to place configurational demands on computer design to avoid limitations imposed by the speed of light, and a Japanese projection of a computer capable of several hundred billion operations/sec is mentioned. The NASA Numerical Aerodynamic Simulator is described, showing capabilities of a billion operations/sec with a memory of 240 million words using existing technology. Near-term advances in fluid dynamics are discussed.

Interferometric characterization of tear film dynamics

NASA Astrophysics Data System (ADS)

Primeau, Brian Christopher

The anterior refracting surface of the eye is the thin tear film that forms on the surface of the cornea. When a contact lens is on worn, the tear film covers the contact lens as it would a bare cornea, and is affected by the contact lens material properties. Tear film irregularity can cause both discomfort and vision quality degradation. Under normal conditions, the tear film is less than 10 microns thick and the thickness and topography change in the time between blinks. In order to both better understand the tear film, and to characterize how contact lenses affect tear film behavior, two interferometers were designed and built to separately measure tear film behavior in vitro and in vivo. An in vitro method of characterizing dynamic fluid layers applied to contact lenses mounted on mechanical substrates has been developed using a phase-shifting Twyman-Green interferometer. This interferometer continuously measures light reflected from the surface of the fluid layer, allowing precision analysis of the dynamic fluid layer. Movies showing this fluid layer behavior can be generated. The fluid behavior on the contact lens surface is measured, allowing quantitative analysis beyond what typical contact angle or visual inspection methods provide. The in vivo interferometer is a similar system, with additional modules included to provide capability for human testing. This tear film measurement allows analysis beyond capabilities of typical fluorescein visual inspection or videokeratometry and provides better sensitivity and resolution than shearing interferometry methods.

Deformation of a Capsule in a Power-Law Shear Flow

PubMed Central

2016-01-01

An immersed boundary-lattice Boltzmann method is developed for fluid-structure interactions involving non-Newtonian fluids (e.g., power-law fluid). In this method, the flexible structure (e.g., capsule) dynamics and the fluid dynamics are coupled by using the immersed boundary method. The incompressible viscous power-law fluid motion is obtained by solving the lattice Boltzmann equation. The non-Newtonian rheology is achieved by using a shear rate-dependant relaxation time in the lattice Boltzmann method. The non-Newtonian flow solver is then validated by considering a power-law flow in a straight channel which is one of the benchmark problems to validate an in-house solver. The numerical results present a good agreement with the analytical solutions for various values of power-law index. Finally, we apply this method to study the deformation of a capsule in a power-law shear flow by varying the Reynolds number from 0.025 to 0.1, dimensionless shear rate from 0.004 to 0.1, and power-law index from 0.2 to 1.8. It is found that the deformation of the capsule increases with the power-law index for different Reynolds numbers and nondimensional shear rates. In addition, the Reynolds number does not have significant effect on the capsule deformation in the flow regime considered. Moreover, the power-law index effect is stronger for larger dimensionless shear rate compared to smaller values. PMID:27840656

Cellular automatons applied to gas dynamic problems

NASA Technical Reports Server (NTRS)

Long, Lyle N.; Coopersmith, Robert M.; Mclachlan, B. G.

1987-01-01

This paper compares the results of a relatively new computational fluid dynamics method, cellular automatons, with experimental data and analytical results. This technique has been shown to qualitatively predict fluidlike behavior; however, there have been few published comparisons with experiment or other theories. Comparisons are made for a one-dimensional supersonic piston problem, Stokes first problem, and the flow past a normal flat plate. These comparisons are used to assess the ability of the method to accurately model fluid dynamic behavior and to point out its limitations. Reasonable results were obtained for all three test cases, but the fundamental limitations of cellular automatons are numerous. It may be misleading, at this time, to say that cellular automatons are a computationally efficient technique. Other methods, based on continuum or kinetic theory, would also be very efficient if as little of the physics were included.

Computational fluid dynamics modelling of hydraulics and sedimentation in process reactors during aeration tank settling.

PubMed

Jensen, M D; Ingildsen, P; Rasmussen, M R; Laursen, J

2006-01-01

Aeration tank settling is a control method allowing settling in the process tank during high hydraulic load. The control method is patented. Aeration tank settling has been applied in several waste water treatment plants using the present design of the process tanks. Some process tank designs have shown to be more effective than others. To improve the design of less effective plants, computational fluid dynamics (CFD) modelling of hydraulics and sedimentation has been applied. This paper discusses the results at one particular plant experiencing problems with partly short-circuiting of the inlet and outlet causing a disruption of the sludge blanket at the outlet and thereby reducing the retention of sludge in the process tank. The model has allowed us to establish a clear picture of the problems arising at the plant during aeration tank settling. Secondly, several process tank design changes have been suggested and tested by means of computational fluid dynamics modelling. The most promising design changes have been found and reported.

Predicting Flow Reversals in a Computational Fluid Dynamics Simulated Thermosyphon Using Data Assimilation.

PubMed

Reagan, Andrew J; Dubief, Yves; Dodds, Peter Sheridan; Danforth, Christopher M

2016-01-01

A thermal convection loop is a annular chamber filled with water, heated on the bottom half and cooled on the top half. With sufficiently large forcing of heat, the direction of fluid flow in the loop oscillates chaotically, dynamics analogous to the Earth's weather. As is the case for state-of-the-art weather models, we only observe the statistics over a small region of state space, making prediction difficult. To overcome this challenge, data assimilation (DA) methods, and specifically ensemble methods, use the computational model itself to estimate the uncertainty of the model to optimally combine these observations into an initial condition for predicting the future state. Here, we build and verify four distinct DA methods, and then, we perform a twin model experiment with the computational fluid dynamics simulation of the loop using the Ensemble Transform Kalman Filter (ETKF) to assimilate observations and predict flow reversals. We show that using adaptively shaped localized covariance outperforms static localized covariance with the ETKF, and allows for the use of less observations in predicting flow reversals. We also show that a Dynamic Mode Decomposition (DMD) of the temperature and velocity fields recovers the low dimensional system underlying reversals, finding specific modes which together are predictive of reversal direction.

Predicting Flow Reversals in a Computational Fluid Dynamics Simulated Thermosyphon Using Data Assimilation

PubMed Central

Reagan, Andrew J.; Dubief, Yves; Dodds, Peter Sheridan; Danforth, Christopher M.

2016-01-01

A thermal convection loop is a annular chamber filled with water, heated on the bottom half and cooled on the top half. With sufficiently large forcing of heat, the direction of fluid flow in the loop oscillates chaotically, dynamics analogous to the Earth’s weather. As is the case for state-of-the-art weather models, we only observe the statistics over a small region of state space, making prediction difficult. To overcome this challenge, data assimilation (DA) methods, and specifically ensemble methods, use the computational model itself to estimate the uncertainty of the model to optimally combine these observations into an initial condition for predicting the future state. Here, we build and verify four distinct DA methods, and then, we perform a twin model experiment with the computational fluid dynamics simulation of the loop using the Ensemble Transform Kalman Filter (ETKF) to assimilate observations and predict flow reversals. We show that using adaptively shaped localized covariance outperforms static localized covariance with the ETKF, and allows for the use of less observations in predicting flow reversals. We also show that a Dynamic Mode Decomposition (DMD) of the temperature and velocity fields recovers the low dimensional system underlying reversals, finding specific modes which together are predictive of reversal direction. PMID:26849061

On the role of fluids in stick-slip dynamics of saturated granular fault gouge using a coupled computational fluid dynamics-discrete element approach: STICK-SLIP IN SATURATED FAULT GOUGE

DOE Office of Scientific and Technical Information (OSTI.GOV)

Dorostkar, Omid; Guyer, Robert A.; Johnson, Paul A.

The presence of fault gouge has considerable influence on slip properties of tectonic faults and the physics of earthquake rupture. The presence of fluids within faults also plays a significant role in faulting and earthquake processes. In this study, we present 3-D discrete element simulations of dry and fluid-saturated granular fault gouge and analyze the effect of fluids on stick-slip behavior. Fluid flow is modeled using computational fluid dynamics based on the Navier-Stokes equations for an incompressible fluid and modified to take into account the presence of particles. Analysis of a long time train of slip events shows that themore » (1) drop in shear stress, (2) compaction of granular layer, and (3) the kinetic energy release during slip all increase in magnitude in the presence of an incompressible fluid, compared to dry conditions. We also observe that on average, the recurrence interval between slip events is longer for fluid-saturated granular fault gouge compared to the dry case. This observation is consistent with the occurrence of larger events in the presence of fluid. It is found that the increase in kinetic energy during slip events for saturated conditions can be attributed to the increased fluid flow during slip. Finally, our observations emphasize the important role that fluid flow and fluid-particle interactions play in tectonic fault zones and show in particular how discrete element method (DEM) models can help understand the hydromechanical processes that dictate fault slip.« less

On the role of fluids in stick-slip dynamics of saturated granular fault gouge using a coupled computational fluid dynamics-discrete element approach: STICK-SLIP IN SATURATED FAULT GOUGE

DOE PAGES

Dorostkar, Omid; Guyer, Robert A.; Johnson, Paul A.; ...

2017-05-01

The presence of fault gouge has considerable influence on slip properties of tectonic faults and the physics of earthquake rupture. The presence of fluids within faults also plays a significant role in faulting and earthquake processes. In this study, we present 3-D discrete element simulations of dry and fluid-saturated granular fault gouge and analyze the effect of fluids on stick-slip behavior. Fluid flow is modeled using computational fluid dynamics based on the Navier-Stokes equations for an incompressible fluid and modified to take into account the presence of particles. Analysis of a long time train of slip events shows that themore » (1) drop in shear stress, (2) compaction of granular layer, and (3) the kinetic energy release during slip all increase in magnitude in the presence of an incompressible fluid, compared to dry conditions. We also observe that on average, the recurrence interval between slip events is longer for fluid-saturated granular fault gouge compared to the dry case. This observation is consistent with the occurrence of larger events in the presence of fluid. It is found that the increase in kinetic energy during slip events for saturated conditions can be attributed to the increased fluid flow during slip. Finally, our observations emphasize the important role that fluid flow and fluid-particle interactions play in tectonic fault zones and show in particular how discrete element method (DEM) models can help understand the hydromechanical processes that dictate fault slip.« less

Fluid mechanics of dynamic stall. II - Prediction of full scale characteristics

NASA Technical Reports Server (NTRS)

Ericsson, L. E.; Reding, J. P.

1988-01-01

Analytical extrapolations are made from experimental subscale dynamics to predict full scale characteristics of dynamic stall. The method proceeds by establishing analytic relationships between dynamic and static aerodynamic characteristics induced by viscous flow effects. The method is then validated by predicting dynamic test results on the basis of corresponding static test data obtained at the same subscale flow conditions, and the effect of Reynolds number on the static aerodynamic characteristics are determined from subscale to full scale flow conditions.

Finite Element Modeling of Non-linear Coupled Interacting Fault System

NASA Astrophysics Data System (ADS)

Xing, H. L.; Zhang, J.; Wyborn, D.

2009-04-01

PANDAS - Parallel Adaptive static/dynamic Nonlinear Deformation Analysis System - a novel supercomputer simulation tool is developed for simulating the highly non-linear coupled geomechanical-fluid flow-thermal systems involving heterogeneously fractured geomaterials. PANDAS includes the following key components: Pandas/Pre, ESyS_Crustal, Pandas/Thermo, Pandas/Fluid and Pandas/Post as detailed in the following: • Pandas/Pre is developed to visualise the microseismicity events recorded during the hydraulic stimulation process to further evaluate the fracture location and evolution and geological setting of a certain reservoir, and then generate the mesh by it and/or other commercial graphics software (such as Patran) for the further finite element analysis of various cases; The Delaunay algorithm is applied as a suitable method for mesh generation using such a point set; • ESyS_Crustal is a finite element code developed for the interacting fault system simulation, which employs the adaptive static/dynamic algorithm to simulate the dynamics and evolution of interacting fault systems and processes that are relevant on short to mediate time scales in which several dynamic phenomena related with stick-slip instability along the faults need to be taken into account, i.e. (a). slow quasi-static stress accumulation, (b) rapid dynamic rupture, (c) wave propagation and (d) corresponding stress redistribution due to the energy release along the multiple fault boundaries; those are needed to better describe ruputure/microseimicity/earthquake related phenomena with applications in earthquake forecasting, hazard quantification, exploration, and environmental problems. It has been verified with various available experimental results[1-3]; • Pandas/Thermo is a finite element method based module for the thermal analysis of the fractured porous media; the temperature distribution is calculated from the heat transfer induced by the thermal boundary conditions without/with the coupled fluid effects and the geomechanical energy conversion for the pure/coupled thermal analysis. • Pandas/Fluid is a finite element method based module for simulating the fluid flow in the fractured porous media; the fluid flow velocity and pressure are calculated from energy equilibrium equations without/together with the coupling effects of the thermal and solid rock deformation for an independent/coupled fluid flow analysis; • Pandas/Post is to visualise the simulation results through the integration of VTK and/or Patran. All the above modules can be used independently/together to simulate individual/coupled phenomena (such as interacting fault system dynamics, heat flow and fluid flow) without/with coupling effects. PANDAS has been applied to the following issues: • visualisation of the microseismic events to monitor and determine where/how the underground rupture proceeds during a hydraulic stimulation, to generate the mesh using the recorded data for determining the domain of the ruptured zone and to evaluate the material parameters (i.e. the permeability) for the further numerical analysis; • interacting fault system simulation to determine the relevant complicated dynamic rupture process. • geomechanical-fluid flow coupling analysis to investigate the interactions between fluid flow and deformation in the fractured porous media under different loading conditions. • thermo-fluid flow coupling analysis of a fractured geothermal reservoir system. PANDAS will be further developed for a multiscale simulation of multiphase dynamic behaviour for a certain fractured geothermal reservoir. More details and additional application examples will be given during the presentation. References [1] Xing, H. L., Makinouchi, A. and Mora, P. (2007). Finite element modeling of interacting fault system, Physics of the Earth and Planetary Interiors, 163, 106-121.doi:10.1016/j.pepi.2007.05.006 [2] Xing, H. L., Mora, P., Makinouchi, A. (2006). An unified friction description and its application to simulation of frictional instability using finite element method. Philosophy Magazine, 86, 3453-3475 [3] Xing, H. L., Mora, P.(2006). Construction of an intraplate fault system model of South Australia, and simulation tool for the iSERVO institute seed project.. Pure and Applied Geophysics. 163, 2297-2316. DOI 10.1007/s00024-006-0127-x

Marc Henry de Frahan | NREL

Science.gov Websites

Computing Project, Marc develops high-fidelity turbulence models to enhance simulation accuracy and efficient numerical algorithms for future high performance computing hardware architectures. Research Interests High performance computing High order numerical methods for computational fluid dynamics Fluid

Attractors of equations of non-Newtonian fluid dynamics

NASA Astrophysics Data System (ADS)

Zvyagin, V. G.; Kondrat'ev, S. K.

2014-10-01

This survey describes a version of the trajectory-attractor method, which is applied to study the limit asymptotic behaviour of solutions of equations of non-Newtonian fluid dynamics. The trajectory-attractor method emerged in papers of the Russian mathematicians Vishik and Chepyzhov and the American mathematician Sell under the condition that the corresponding trajectory spaces be invariant under the translation semigroup. The need for such an approach was caused by the fact that for many equations of mathematical physics for which the Cauchy initial-value problem has a global (weak) solution with respect to the time, the uniqueness of such a solution has either not been established or does not hold. In particular, this is the case for equations of fluid dynamics. At the same time, trajectory spaces invariant under the translation semigroup could not be constructed for many equations of non-Newtonian fluid dynamics. In this connection, a different approach to the construction of trajectory attractors for dissipative systems was proposed in papers of Zvyagin and Vorotnikov without using invariance of trajectory spaces under the translation semigroup and is based on the topological lemma of Shura-Bura. This paper presents examples of equations of non-Newtonian fluid dynamics (the Jeffreys system describing movement of the Earth's crust, the model of motion of weak aqueous solutions of polymers, a system with memory) for which the aforementioned construction is used to prove the existence of attractors in both the autonomous and the non-autonomous cases. At the beginning of the paper there is also a brief exposition of the results of Ladyzhenskaya on the existence of attractors of the two-dimensional Navier-Stokes system and the result of Vishik and Chepyzhov for the case of attractors of the three-dimensional Navier-Stokes system. Bibliography: 34 titles.

Computer animation challenges for computational fluid dynamics

NASA Astrophysics Data System (ADS)

Vines, Mauricio; Lee, Won-Sook; Mavriplis, Catherine

2012-07-01

Computer animation requirements differ from those of traditional computational fluid dynamics (CFD) investigations in that visual plausibility and rapid frame update rates trump physical accuracy. We present an overview of the main techniques for fluid simulation in computer animation, starting with Eulerian grid approaches, the Lattice Boltzmann method, Fourier transform techniques and Lagrangian particle introduction. Adaptive grid methods, precomputation of results for model reduction, parallelisation and computation on graphical processing units (GPUs) are reviewed in the context of accelerating simulation computations for animation. A survey of current specific approaches for the application of these techniques to the simulation of smoke, fire, water, bubbles, mixing, phase change and solid-fluid coupling is also included. Adding plausibility to results through particle introduction, turbulence detail and concentration on regions of interest by level set techniques has elevated the degree of accuracy and realism of recent animations. Basic approaches are described here. Techniques to control the simulation to produce a desired visual effect are also discussed. Finally, some references to rendering techniques and haptic applications are mentioned to provide the reader with a complete picture of the challenges of simulating fluids in computer animation.

Synchrotron phase-contrast X-ray imaging reveals fluid dosing dynamics for gene transfer into mouse airways.

PubMed

Donnelley, M; Siu, K K W; Jamison, R A; Parsons, D W

2012-01-01

Although airway gene transfer research in mouse models relies on bolus fluid dosing into the nose or trachea, the dynamics and immediate fate of delivered gene transfer agents are poorly understood. In particular, this is because there are no in vivo methods able to accurately visualize the movement of fluid in small airways of intact animals. Using synchrotron phase-contrast X-ray imaging, we show that the fate of surrogate fluid doses delivered into live mouse airways can now be accurately and non-invasively monitored with high spatial and temporal resolution. This new imaging approach can help explain the non-homogenous distributions of gene expression observed in nasal airway gene transfer studies, suggests that substantial dose losses may occur at deliver into mouse trachea via immediate retrograde fluid motion and shows the influence of the speed of bolus delivery on the relative targeting of conducting and deeper lung airways. These findings provide insight into some of the factors that can influence gene expression in vivo, and this method provides a new approach to documenting and analyzing dose delivery in small-animal models.

Exploring a multi-scale method for molecular simulation in continuum solvent model: Explicit simulation of continuum solvent as an incompressible fluid.

PubMed

Xiao, Li; Luo, Ray

2017-12-07

We explored a multi-scale algorithm for the Poisson-Boltzmann continuum solvent model for more robust simulations of biomolecules. In this method, the continuum solvent/solute interface is explicitly simulated with a numerical fluid dynamics procedure, which is tightly coupled to the solute molecular dynamics simulation. There are multiple benefits to adopt such a strategy as presented below. At this stage of the development, only nonelectrostatic interactions, i.e., van der Waals and hydrophobic interactions, are included in the algorithm to assess the quality of the solvent-solute interface generated by the new method. Nevertheless, numerical challenges exist in accurately interpolating the highly nonlinear van der Waals term when solving the finite-difference fluid dynamics equations. We were able to bypass the challenge rigorously by merging the van der Waals potential and pressure together when solving the fluid dynamics equations and by considering its contribution in the free-boundary condition analytically. The multi-scale simulation method was first validated by reproducing the solute-solvent interface of a single atom with analytical solution. Next, we performed the relaxation simulation of a restrained symmetrical monomer and observed a symmetrical solvent interface at equilibrium with detailed surface features resembling those found on the solvent excluded surface. Four typical small molecular complexes were then tested, both volume and force balancing analyses showing that these simple complexes can reach equilibrium within the simulation time window. Finally, we studied the quality of the multi-scale solute-solvent interfaces for the four tested dimer complexes and found that they agree well with the boundaries as sampled in the explicit water simulations.

Current research in cavitating fluid films

NASA Technical Reports Server (NTRS)

Brewe, D. E. (Editor); Ball, J. H. (Editor); Khonsari, M. M. (Editor)

1990-01-01

A review of the current research of cavitation in fluid films is presented. Phenomena and experimental observations include gaseous cavitation, vapor cavitation, and gas entrainment. Cavitation in flooded, starved, and dynamically loaded journal bearings, as well as squeeze films are reviewed. Observations of cavitation damage in bearings and the possibility of cavitation between parallel plates with microasperities were discussed. The transcavity fluid transport process, meniscus motion and geometry or form of the film during rupture, and reformation were summarized. Performance effects were related to heat transfer models in the cavitated region and hysteresis influence on rotor dynamics coefficients. A number of cavitation algorithms was presented together with solution procedures using the finite difference and finite element methods. Although Newtonian fluids were assumed in most of the discussions, the effect of non-Newtonian fluids on cavitation was also discussed.

Implementing a Loosely Coupled Fluid Structure Interaction Finite Element Model in PHASTA

NASA Astrophysics Data System (ADS)

Pope, David

Fluid Structure Interaction problems are an important multi-physics phenomenon in the design of aerospace vehicles and other engineering applications. A variety of computational fluid dynamics solvers capable of resolving the fluid dynamics exist. PHASTA is one such computational fluid dynamics solver. Enhancing the capability of PHASTA to resolve Fluid-Structure Interaction first requires implementing a structural dynamics solver. The implementation also requires a correction of the mesh used to solve the fluid equations to account for the deformation of the structure. This results in mesh motion and causes the need for an Arbitrary Lagrangian-Eulerian modification to the fluid dynamics equations currently implemented in PHASTA. With the implementation of both structural dynamics physics, mesh correction, and the Arbitrary Lagrangian-Eulerian modification of the fluid dynamics equations, PHASTA is made capable of solving Fluid-Structure Interaction problems.

Computational Fluid Dynamics. [numerical methods and algorithm development

NASA Technical Reports Server (NTRS)

1992-01-01

This collection of papers was presented at the Computational Fluid Dynamics (CFD) Conference held at Ames Research Center in California on March 12 through 14, 1991. It is an overview of CFD activities at NASA Lewis Research Center. The main thrust of computational work at Lewis is aimed at propulsion systems. Specific issues related to propulsion CFD and associated modeling will also be presented. Examples of results obtained with the most recent algorithm development will also be presented.

Isentropic fluid dynamics in a curved pipe

NASA Astrophysics Data System (ADS)

Colombo, Rinaldo M.; Holden, Helge

2016-10-01

In this paper we study isentropic flow in a curved pipe. We focus on the consequences of the geometry of the pipe on the dynamics of the flow. More precisely, we present the solution of the general Cauchy problem for isentropic fluid flow in an arbitrarily curved, piecewise smooth pipe. We consider initial data in the subsonic regime, with small total variation about a stationary solution. The proof relies on the front-tracking method and is based on [1].

Nonlinear Hyperbolic Equations - Theory, Computation Methods, and Applications. Volume 24. Note on Numerical Fluid Mechanics

DTIC Science & Technology

1989-01-01

Calculations and Experiments (B.van den Berg/ D.A. Humphreysl E. Krause /J.P. F. Lindhout) Volume 20 Proceedings of the Seventh GAMM-Conference on...GRID METHODS FOR HYPERBOLIC PROBLEMS Wolfgang Hackbusch Sigrid Hagemann Institut fUr Informatik und Praktische Mathematik Christian-Albrechts...Euler Equations. Proceedings of the 8th Inter- national Conference on Numerical Methods in Fluid Dynamics (E. Krause , ed.), Aachen, 1988. Springer

Effect of wettability on two-phase quasi-static displacement: Validation of two pore scale modeling approaches

NASA Astrophysics Data System (ADS)

Verma, Rahul; Icardi, Matteo; Prodanović, Maša

2018-05-01

Understanding of pore-scale physics for multiphase flow in porous media is essential for accurate description of various flow phenomena. In particular, capillarity and wettability strongly influence capillary pressure-saturation and relative permeability relationships. Wettability is quantified by the contact angle of the fluid-fluid interface at the pore walls. In this work we focus on the non-trivial interface equilibria in presence of non-neutral wetting and complex geometries. We quantify the accuracy of a volume-of-fluid (VOF) formulation, implemented in a popular open-source computational fluid dynamics code, compared with a new formulation of a level set (LS) method, specifically developed for quasi-static capillarity-dominated displacement. The methods are tested in rhomboidal packings of spheres for a range of contact angles and for different rhomboidal configurations and the accuracy is evaluated against the semi-analytical solutions obtained by Mason and Morrow (1994). While the VOF method is implemented in a general purpose code that solves the full Navier-Stokes (NS) dynamics in a finite volume formulation, with additional terms to model surface tension, the LS method is optimized for the quasi-static case and, therefore, less computationally expensive. To overcome the shortcomings of the finite volume NS-VOF system for low capillary number flows, and its computational cost, we introduce an overdamped dynamics and a local time stepping to speed up the convergence to the steady state, for every given imposed pressure gradient (and therefore saturation condition). Despite these modifications, the methods fundamentally differ in the way they capture the interface, as well as in the number of equations solved and in the way the mean curvature (or equivalently capillary pressure) is computed. This study is intended to provide a rigorous validation study and gives important indications on the errors committed by these methods in solving more complex geometry and dynamics, where usually many sources of errors are interplaying.

Design and dynamic modeling of electrorheological fluid-based variable-stiffness fin for robotic fish

NASA Astrophysics Data System (ADS)

Bazaz Behbahani, Sanaz; Tan, Xiaobo

2017-08-01

Fish actively control their stiffness in different swimming conditions. Inspired by such an adaptive behavior, in this paper we study the design, prototyping, and dynamic modeling of compact, tunable-stiffness fins for robotic fish, where electrorheological (ER) fluid serves as the enabling element. A multi-layer composite fin with an ER fluid core is prototyped and utilized to investigate the influence of electrical field on its performance. Hamilton's principle is used to derive the dynamic equations of motion of the flexible fin, and Lighthill's large-amplitude elongated-body theory is adopted to estimate the hydrodynamic force when the fin undergoes base-actuated rotation. The dynamic equations are then discretized using the finite element method, to obtain an approximate numerical solution. Experiments are conducted on the prototyped flexible ER fluid-filled beam for parameter identification and validation of the proposed model, and for examining the effectiveness of electrically controlled stiffness tuning. In particular, it is found that the natural frequency is increased by almost 40% when the applied electric field changes from 0 to 1.5× {10}6 {{V}} {{{m}}}-1.

Data use investigation for the magnetic field satellite (MAGSAT) mission: Geomagnetic field forecasting and fluid dynamics of the core

NASA Technical Reports Server (NTRS)

Benton, E. R. (Principal Investigator)

1982-01-01

MAGSAT data were used to construct a variety of spherical harmonic models of the main geomagnetic field emanating from Earth's liquid core at poch 1980. These models were used to: (1) accurately determine the radius of Earth's core by a magnetic method, (2) calculate estimates, of the long-term ange of variation of geomagnetic Gauss coefficients; (3) establish a preferred truncation level for current spherical harmonic models of the main geomagnetic field from the core; (4) evaluate a method for taking account of electrical conduction in the mantle when the magnetic field is downward continued to the core-mantle boundary; and (5) establish that upwelling and downwelling of fluid motion at the top of the core is probably detectable, observationally. A fluid dynamics forecast model was not produced because of insufficient data.

Cellular Biotechnology Operations Support System Fluid Dynamics Investigation

NASA Technical Reports Server (NTRS)

2003-01-01

Aboard the International Space Station (ISS), the Tissue Culture Medium (TCM) is the bioreactor vessel in which cell cultures are grown. With its two syringe ports, it is much like a bag used to administer intravenous fluid, except it allows gas exchange needed for life. The TCM contains cell culture medium, and when frozen cells are flown to the ISS, they are thawed and introduced to the TCM through the syringe ports. In the Cellular Biotechnology Operations Support System-Fluid Dynamics Investigation (CBOSS-FDI) experiment, several mixing procedures are being assessed to determine which method achieves the most uniform mixing of growing cells and culture medium.

Astrophysical fluid dynamics

NASA Astrophysics Data System (ADS)

Ogilvie, Gordon I.

2016-06-01

> These lecture notes and example problems are based on a course given at the University of Cambridge in Part III of the Mathematical Tripos. Fluid dynamics is involved in a very wide range of astrophysical phenomena, such as the formation and internal dynamics of stars and giant planets, the workings of jets and accretion discs around stars and black holes and the dynamics of the expanding Universe. Effects that can be important in astrophysical fluids include compressibility, self-gravitation and the dynamical influence of the magnetic field that is `frozen in' to a highly conducting plasma. The basic models introduced and applied in this course are Newtonian gas dynamics and magnetohydrodynamics (MHD) for an ideal compressible fluid. The mathematical structure of the governing equations and the associated conservation laws are explored in some detail because of their importance for both analytical and numerical methods of solution, as well as for physical interpretation. Linear and nonlinear waves, including shocks and other discontinuities, are discussed. The spherical blast wave resulting from a supernova, and involving a strong shock, is a classic problem that can be solved analytically. Steady solutions with spherical or axial symmetry reveal the physics of winds and jets from stars and discs. The linearized equations determine the oscillation modes of astrophysical bodies, as well as their stability and their response to tidal forcing.

Aeroelastic analysis of bridge girder section using computer modeling

DOT National Transportation Integrated Search

2001-05-01

This report describes the numerical simulation of wind flow around bridges using the Finite Element Method (FEM) and the principles of Computational Fluid Dynamics (CFD) and Computational Structural Dynamics (CSD). Since, the suspension bridges are p...

Actively controlling coolant-cooled cold plate configuration

DOE Office of Scientific and Technical Information (OSTI.GOV)

Chainer, Timothy J.; Parida, Pritish R.

A method is provided to facilitate active control of thermal and fluid dynamic performance of a coolant-cooled cold plate. The method includes: monitoring a variable associated with at least one of the coolant-cooled cold plate or one or more electronic components being cooled by the cold plate; and dynamically varying, based on the monitored variable, a physical configuration of the cold plate. By dynamically varying the physical configuration, the thermal and fluid dynamic performance of the cold plate are adjusted to, for example, optimally cool the one or more electronic components, and at the same time, reduce cooling power consumptionmore » used in cooling the electronic component(s). The physical configuration can be adjusted by providing one or more adjustable plates within the coolant-cooled cold plate, the positioning of which may be adjusted based on the monitored variable.« less

Scaling in two-fluid pinch-off

NASA Astrophysics Data System (ADS)

Pommer, Chris; Suryo, Ronald; Subramani, Hariprasad; Harris, Michael; Basaran, Osman

2009-11-01

Two-fluid pinch-off is encountered when drops or bubbles of one fluid are ejected from a nozzle into another fluid or when a compound jet breaks. While the breakup of a drop in a passive environment and that of a passive bubble in a liquid are well understood, the physics of pinch-off when both the inner and outer fluids are dynamically active is inadequately understood. In this talk, the breakup of a compound jet whose core and shell are both incompressible Newtonian fluids is analyzed computationally by a method of lines ALE algorithm which uses finite elements with elliptic mesh generation for spatial discretization and adaptive finite differences for time integration. Pinch-off dynamics are investigated well beyond the limit of experiments set by the wavelength of visible light and that of various algorithms used in the literature. Simulations show that the minimum neck radius r initially scales with time τ before breakup as &αcirc; where α varies over a certain range. However, depending on the values of the governing dimensionless groups, this initial scaling regime may be transitory and, closer to pinch-off, the dynamics may transition to a final asymptotic regime for which r ˜&βcirc;, where β!=α.

Fluid-fluid interfacial mobility from random walks

NASA Astrophysics Data System (ADS)

Barclay, Paul L.; Lukes, Jennifer R.

2017-12-01

Dual control volume grand canonical molecular dynamics is used to perform the first calculation of fluid-fluid interfacial mobilities. The mobility is calculated from one-dimensional random walks of the interface by relating the diffusion coefficient to the interfacial mobility. Three different calculation methods are employed: one using the interfacial position variance as a function of time, one using the mean-squared interfacial displacement, and one using the time-autocorrelation of the interfacial velocity. The mobility is calculated for two liquid-liquid interfaces and one liquid-vapor interface to examine the robustness of the methods. Excellent agreement between the three calculation methods is shown for all the three interfaces, indicating that any of them could be used to calculate the interfacial mobility.

Computational fluid dynamics (CFD) insights into agitation stress methods in biopharmaceutical development.

PubMed

Bai, Ge; Bee, Jared S; Biddlecombe, James G; Chen, Quanmin; Leach, W Thomas

2012-02-28

Agitation of small amounts of liquid is performed routinely in biopharmaceutical process, formulation, and packaging development. Protein degradation commonly results from agitation, but the specific stress responsible or degradation mechanism is usually not well understood. Characterization of the agitation stress methods is critical to identifying protein degradation mechanisms or specific sensitivities. In this study, computational fluid dynamics (CFD) was used to model agitation of 1 mL of fluid by four types of common laboratory agitation instruments, including a rotator, orbital shaker, magnetic stirrer and vortex mixer. Fluid stresses in the bulk liquid and near interfaces were identified, quantified and compared. The vortex mixer provides the most intense stresses overall, while the stir bar system presented locally intense shear proximal to the hydrophobic stir bar surface. The rotator provides gentler fluid stresses, but the air-water interfacial area and surface stresses are relatively high given its low rotational frequency. The orbital shaker provides intermediate-level stresses but with the advantage of a large stable platform for consistent vial-to-vial homogeneity. Selection of experimental agitation methods with targeted types and intensities of stresses can facilitate better understanding of protein degradation mechanisms and predictability for "real world" applications. Copyright © 2011 Elsevier B.V. All rights reserved.

Advances in numerical and applied mathematics

NASA Technical Reports Server (NTRS)

South, J. C., Jr. (Editor); Hussaini, M. Y. (Editor)

1986-01-01

This collection of papers covers some recent developments in numerical analysis and computational fluid dynamics. Some of these studies are of a fundamental nature. They address basic issues such as intermediate boundary conditions for approximate factorization schemes, existence and uniqueness of steady states for time dependent problems, and pitfalls of implicit time stepping. The other studies deal with modern numerical methods such as total variation diminishing schemes, higher order variants of vortex and particle methods, spectral multidomain techniques, and front tracking techniques. There is also a paper on adaptive grids. The fluid dynamics papers treat the classical problems of imcompressible flows in helically coiled pipes, vortex breakdown, and transonic flows.

Computational Methods for Stability and Control (COMSAC): The Time Has Come

NASA Technical Reports Server (NTRS)

Hall, Robert M.; Biedron, Robert T.; Ball, Douglas N.; Bogue, David R.; Chung, James; Green, Bradford E.; Grismer, Matthew J.; Brooks, Gregory P.; Chambers, Joseph R.

2005-01-01

Powerful computational fluid dynamics (CFD) tools have emerged that appear to offer significant benefits as an adjunct to the experimental methods used by the stability and control community to predict aerodynamic parameters. The decreasing costs for and increasing availability of computing hours are making these applications increasingly viable as time goes on and the cost of computing continues to drop. This paper summarizes the efforts of four organizations to utilize high-end computational fluid dynamics (CFD) tools to address the challenges of the stability and control arena. General motivation and the backdrop for these efforts will be summarized as well as examples of current applications.

Particles size distribution in diluted magnetic fluids

NASA Astrophysics Data System (ADS)

Yerin, Constantine V.

2017-06-01

Changes in particles and aggregates size distribution in diluted kerosene based magnetic fluids is studied by dynamic light scattering method. It has been found that immediately after dilution in magnetic fluids the system of aggregates with sizes ranging from 100 to 250-1000 nm is formed. In 50-100 h after dilution large aggregates are peptized and in the sample stationary particles and aggregates size distribution is fixed.

Lubricated immersed boundary method in two dimensions

NASA Astrophysics Data System (ADS)

Fai, Thomas G.; Rycroft, Chris H.

2018-03-01

Many biological examples of fluid-structure interaction, including the transit of red blood cells through the narrow slits in the spleen and the intracellular trafficking of vesicles into dendritic spines, involve the near-contact of elastic structures separated by thin layers of fluid. Motivated by such problems, we introduce an immersed boundary method that uses elements of lubrication theory to resolve thin fluid layers between immersed boundaries. We demonstrate 2nd-order accurate convergence for simple two-dimensional flows with known exact solutions to showcase the increased accuracy of this method compared to the standard immersed boundary method. Motivated by the phenomenon of wall-induced migration, we apply the lubricated immersed boundary method to simulate an elastic vesicle near a wall in shear flow. We also simulate the dynamics of a vesicle traveling through a narrow channel and observe the ability of the lubricated method to capture the vesicle motion on relatively coarse fluid grids.

Preliminary Numerical Simulations of Nozzle Formation in the Host Rock of Supersonic Volcanic Jets

NASA Astrophysics Data System (ADS)

Wohletz, K. H.; Ogden, D. E.; Glatzmaier, G. A.

2006-12-01

Recognizing the difficulty in quantitatively predicting how a vent changes during an explosive eruption, Kieffer (Kieffer, S.W., Rev. Geophys. 27, 1989) developed the theory of fluid dynamic nozzles for volcanism, utilizing a highly developed predictive scheme used extensively in aerodynamics for design of jet and rocket nozzles. Kieffer's work shows that explosive eruptions involve flow from sub to supersonic conditions through the vent and that these conditions control the erosion of the vent to nozzle shapes and sizes that maximize mass flux. The question remains how to predict the failure and erosion of vent host rocks by a high-speed, multiphase, compressible fluid that represents an eruption column. Clearly, in order to have a quantitative model of vent dynamics one needs a robust computational method for a turbulent, compressible, multiphase fluid. Here we present preliminary simulations of fluid flowing from a high-pressure reservoir through an eroding conduit and into the atmosphere. The eruptive fluid is modeled as an ideal gas, the host rock as a simple incompressible fluid with sandstone properties. Although these simulations do not yet include the multiphase dynamics of the eruptive fluid or the solid mechanics of the host rock, the evolution of the host rock into a supersonic nozzle is clearly seen. Our simulations show shock fronts both above the conduit, where the gas has expanded into the atmosphere, and within the conduit itself, thereby influencing the dynamics of the jet decompression.

An Immersed Boundary Method for Solving the Compressible Navier-Stokes Equations with Fluid Structure Interaction

NASA Technical Reports Server (NTRS)

Brehm, Christoph; Barad, Michael F.; Kiris, Cetin C.

2016-01-01

An immersed boundary method for the compressible Navier-Stokes equation and the additional infrastructure that is needed to solve moving boundary problems and fully coupled fluid-structure interaction is described. All the methods described in this paper were implemented in NASA's LAVA solver framework. The underlying immersed boundary method is based on the locally stabilized immersed boundary method that was previously introduced by the authors. In the present paper this method is extended to account for all aspects that are involved for fluid structure interaction simulations, such as fast geometry queries and stencil computations, the treatment of freshly cleared cells, and the coupling of the computational fluid dynamics solver with a linear structural finite element method. The current approach is validated for moving boundary problems with prescribed body motion and fully coupled fluid structure interaction problems in 2D and 3D. As part of the validation procedure, results from the second AIAA aeroelastic prediction workshop are also presented. The current paper is regarded as a proof of concept study, while more advanced methods for fluid structure interaction are currently being investigated, such as geometric and material nonlinearities, and advanced coupling approaches.

Bioassay of safinamide in biological fluids of humans and various animal species.

PubMed

Dal Bo, Lorenzo; Mazzucchelli, Paolo; Fibbioli, Monia; Marzo, Antonio

2006-01-01

This paper describes three methods to bioassay safinamide (CAS 133865-89-1) in biological fluids of humans and laboratory animals for pharmacokinetic, toxicokinetic and bioavailability studies. Two methods profited from liquid chromatography tandem mass spectrometry (LC-MS-MS) system, one (micro bioassay) working in the low dynamic range 0.5-20 ng/ml, the other in the range 20-6000 ng/ml. A third method used high-performance liquid chromatrography with fluorimetric detection (HPLC-FD), working in the dynamic range 20-1000 ng/ml. All the three methods were validated in compliance with the accreditated views on analytical bioassays. The shorter run time (5.5 min vs 16 min) has led the authors to prefer the two LC-MS-MS methods to the HPLC-FD bioassay, even if all the performances of the three methods were excellent. The methods described in this paper were and are still now extensively used to assay safinamide in more than 10,000 specimens of biological fluids from humans and laboratory animals in the development program of the drug. Main pharmacokinetic results obtained in various Phase I trials on healthy volunteers are briefly reported.

Molecular dynamics simulation of a needle-sphere binary mixture

NASA Astrophysics Data System (ADS)

Raghavan, Karthik

This paper investigates the dynamic behaviour of a hard needle-sphere binary system using a novel numerical technique called the Newton homotopy continuation (NHC) method. This mixture is representative of a polymer melt where both long chain molecules and monomers coexist. Since the intermolecular forces are generated from hard body interactions, the consequence of missed collisions or incorrect collision sequences have a significant bearing on the dynamic properties of the fluid. To overcome this problem, in earlier work NHC was chosen over traditional Newton-Raphson methods to solve the hard body dynamics of a needle fluid in random media composed of overlapping spheres. Furthermore, the simplicity of interactions and dynamics allows us to focus our research directly on the effects of particle shape and density on the transport behaviour of the mixture. These studies are also compared with earlier works that examined molecular chains in porous media primarily to understand the differences in molecular transport in the bulk versus porous systems.

B-spline Method in Fluid Dynamics

NASA Technical Reports Server (NTRS)

Botella, Olivier; Shariff, Karim; Mansour, Nagi N. (Technical Monitor)

2001-01-01

B-spline functions are bases for piecewise polynomials that possess attractive properties for complex flow simulations : they have compact support, provide a straightforward handling of boundary conditions and grid nonuniformities, and yield numerical schemes with high resolving power, where the order of accuracy is a mere input parameter. This paper reviews the progress made on the development and application of B-spline numerical methods to computational fluid dynamics problems. Basic B-spline approximation properties is investigated, and their relationship with conventional numerical methods is reviewed. Some fundamental developments towards efficient complex geometry spline methods are covered, such as local interpolation methods, fast solution algorithms on cartesian grid, non-conformal block-structured discretization, formulation of spline bases of higher continuity over triangulation, and treatment of pressure oscillations in Navier-Stokes equations. Application of some of these techniques to the computation of viscous incompressible flows is presented.

Computational Fluid Dynamics Simulation Study of Active Power Control in Wind Plants

DOE Office of Scientific and Technical Information (OSTI.GOV)

Fleming, Paul; Aho, Jake; Gebraad, Pieter

2016-08-01

This paper presents an analysis performed on a wind plant's ability to provide active power control services using a high-fidelity computational fluid dynamics-based wind plant simulator. This approach allows examination of the impact on wind turbine wake interactions within a wind plant on performance of the wind plant controller. The paper investigates several control methods for improving performance in waked conditions. One method uses wind plant wake controls, an active field of research in which wind turbine control systems are coordinated to account for their wakes, to improve the overall performance. Results demonstrate the challenge of providing active power controlmore » in waked conditions but also the potential methods for improving this performance.« less

Cerebrospinal fluid dynamics study in communicating hydrocephalus

PubMed Central

Ramesh, Vengalathur Ganesan; Narasimhan, Vidhya; Balasubramanian, Chandramouli

2017-01-01

Context: Communicating hydrocephalus often poses a challenge in diagnosis and management decisions. Aims: The objective of this study is to measure the opening pressure (Po), pressure volume index (PVI), and cerebrospinal fluid outflow resistance (Rout), in patients with communicating hydrocephalus using bolus lumbar injection method and to evaluate its diagnostic and prognostic value. Materials and Methods: The study was conducted in 50 patients with communicating hydrocephalus, including normal pressure hydrocephalus (NPH) (19), post-meningitic hydrocephalus (23) and post-traumatic hydrocephalus (8). An improvised bolus lumbar injection method [the Madras Institute of Neurology (MIN) method] was used. Results: In the NPH Group, the CSF dynamics studies correlated well with the clinico-radiological classification. The prediction of shunt responsiveness by CSF dynamics studies correlated with good outcome in 87.5%. In the post-meningitic hydrocephalus group, the value of CSF dynamics studies in predicting patients needing shunt was 89.5%. The CSF dynamics studies detected patients who needed shunt earlier than clinical or radiological indications. In the post-traumatic hydrocephalus group, 62.5% of patients improved with the treatment based on CSF dynamics studies. Conclusions: The improvised bolus lumbar injection method (MIN method) is a very simple test with fairly reliable and reproducible results. Study of CSF dynamics is a valuable tool in communicating hydrocephalus for confirmation of diagnosis and predicting shunt responsiveness. This is the first time that the value of CSF dynamics has been studied in patients with post-meningitic hydrocephalus. It was also useful for early selection of cases for shunting and for identifying patients with atrophic ventriculomegaly, thereby avoiding unnecessary shunt. PMID:28484522

Modeling of viscoelastic properties of nonpermeable porous rocks saturated with highly viscous fluid at seismic frequencies at the core scale

NASA Astrophysics Data System (ADS)

Wang, Zizhen; Schmitt, Douglas R.; Wang, Ruihe

2017-08-01

A core scale modeling method for viscoelastic properties of rocks saturated with viscous fluid at low frequencies is developed based on the stress-strain method. The elastic moduli dispersion of viscous fluid is described by the Maxwell's spring-dash pot model. Based on this modeling method, we numerically test the effects of frequency, fluid viscosity, porosity, pore size, and pore aspect ratio on the storage moduli and the stress-strain phase lag of saturated rocks. And we also compared the modeling results to the Hashin-Shtrikman bounds and the coherent potential approximation (CPA). The dynamic moduli calculated from the modeling are lower than the predictions of CPA, and both of these fall between the Hashin-Shtrikman bounds. The modeling results indicate that the frequency and the fluid viscosity have similar effects on the dynamic moduli dispersion of fully saturated rocks. We observed the Debye peak in the phase lag variation with the change of frequency and viscosity. The pore structure parameters, such as porosity, pore size, and aspect ratio affect the rock frame stiffness and result in different viscoelastic behaviors of the saturated rocks. The stress-strain phase lags are larger with smaller stiffness contrasts between the rock frame and the pore fluid. The viscoelastic properties of saturated rocks are more sensitive to aspect ratio compared to other pore structure parameters. The results suggest that significant seismic dispersion (at about 50-200 Hz) might be expected for both compressional and shear waves passing through rocks saturated with highly viscous fluids.