Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fluid Structure Interaction”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Discovery of two-dimensional binary nanoparticle superlattices using global Monte Carlo optimization

Binary nanoparticle (NP) superlattices exhibit distinct collective plasmonic, magnetic, optical, and electronic properties. Here, we computationally demonstrate how fluid-fluid interfaces could be used to self-assemble binary systems of NPs into 2D superlattices when the NP species exhibit different miscibility with the fluids forming the interface. We develop a basin-hopping Monte Carlo (BHMC) algorithm tailored for interface-trapped structures to rapidly determine the ground-state configuration of NPs, allowing us to explore the repertoire of binary NP architectures formed at the interface. By varying the NP size ratio, interparticle interaction strength, and difference in NP miscibility with the two fluids, we demonstrate the assembly of an array of exquisite 2D periodic architectures, including AB-, AB 2 -, and AB 3 -type monolayer superlattices as well as AB-, AB 2 -, A 3 B 5 -, and A 4 B 6 -type bilayer superlattices. Our results suggest that the interfacial assembly approach could be a versatile platform for fabricating 2D colloidal superlattices with tunable structure and properties.

36 MATERIALS SCIENCE↗
NASA NTRS2021

Particle Acceleration by Dense Impulsive Structures Moving in Ambient Magnetospheric Plasma. 3-D Hybrid Kinetic Modeling and MMS Observations

High resolution observations of dense plasma impulsive structures moving through an ambient background magnetospheric flows were captured by the Magnetospheric Multiscale mission. The observations show particle heating and acceleration, shock-like wave formation, and whistler wave excitation inside the interface between the dense impulsive plasma structures and the ambient plasma. A multiscale hybrid kinetic simulation provides an explanation of the observed wave-particle interactions with the assumption that the dense plasma structures may be represented by plasma clouds which are formed at the magnetopause layer due to reconnection processes. Plain Language Summary: Dense, impulsive plasma structures moving through a background plasma were captured by the NASA Magnetospheric Multiscale mission. The observations show that the dense structures can generate strong perturbations in the electromagnetic field and shock-like waves. Interactions between these electromagnetic waves and the particles results in particle acceleration. 3-D hybrid kinetic modeling (particle description for ions and fluid description for electrons) was used to investigate the plasma physics of the observed structures. It was assumed that the plasma clouds were formed by magnetic field reconnection inside the magnetopause, which is the interface between the solar wind particles and the cold low-density magnetospheric plasma. The work helps us understand the plasma environment at the interface between the Earth and solar wind, near planetary moons, within astrophysical explosions, and possibly at the interface between the solar wind and local interstellar medium.

A. S. Lipatov↗
DOE OSTI2014

A new and effective method for thermostatting confined fluids

We present a simple thermostatting method suitable for nanoconfined fluid systems. Two conventional strategies involve thermostatting the fluid directly or employing a thermal wall that couples only the wall atoms with the thermostat. When only a thermal wall is implemented, the temperature control of the fluid is true to the actual experiment and the heat is transferred from the fluid to the walls. However, for large or complex systems it can often be computationally prohibitive to employ thermal walls. To overcome this limitation many researchers choose to freeze wall atoms and instead apply a synthetic thermostat to the fluid directly through the equations of motion. This, however, can have serious consequences for the mechanical, thermodynamic, and dynamical properties of the fluid by introducing unphysical behaviour into the system [Bernardi et al., J. Chem. Phys. 132, 244706 (2010)]. In this paper, we propose a simple scheme which enables working with both frozen walls and naturally thermostatted liquids. This is done by superimposing the walls with oscillating particles, which vibrate on the edge of the fluid control volume. These particles exchange energy with the fluid molecules, but do not interact with wall atoms or each other, thus behaving as virtual particles. Their displacements violate the Lindemann criterion for melting, in such a way that the net effect would not amount to an additional confining surface. One advantage over standard techniques is the reduced computational cost, particularly for large walls, since they can be kept rigid. Another advantage over accepted strategies is the opportunity to freeze complex charged walls such as β-cristobalite. The method furthermore overcomes the problem with polar fluids such as water, as thermalized charged surfaces require higher spring constants to preserve structural stability, due to the effects of strong Coulomb interactions, thus inevitably degrading the thermostatting efficiency.

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

Small Angle Neutron Scattering Study of Complex Coacervate Micelles and Hydrogels Formed from Ionic Diblock and Triblock Copolymers

A complex coacervate is a fluid phase that results from the electrostatic interactions between two oppositely charged macromolecules. The nature of the coacervate core structure of hydrogels and micelles formed from complexation between pairs of diblock or triblock copolymers containing oppositely charged end-blocks as a function of polymer and salt concentration was investigated. Both ABA triblock copolymers of poly[(allyl glycidyl ether)-b-(ethylene oxide)-b-(allyl glycidyl ether)] and analogous poly[(allyl glycidyl ether)-b-(ethylene oxide)] diblock copolymers, which were synthesized to be nearly one-half of the symmetrical triblock copolymers, were studied. The poly(allyl glycidyl ether) blocks were functionalized with either guanidinium or sulfonate groups via postpolymerization modification. Mixing of oppositely charged block copolymers resulted in the formation of nanometer-scale coacervate domains. Small angle neutron scattering (SANS) experiments were used to investigate the size and spacing of the coacervate domains. The SANS patterns were fit using a previously vetted, detailed model consisting of polydisperse core–shell micelles with a randomly distributed sphere or body-centered cubic (BCC) structure factor. For increasing polymer concentration, the size of the coacervate domains remained constant while the spatial extent of the poly(ethylene oxide) (PEO) corona decreased. However, increasing salt concentration resulted in a decrease in both the coacervate domain size and the corona size due to a combination of the electrostatic interactions being screened and the shrinkage of the neutral PEO blocks. Additionally, for the triblock copolymers that formed BCC ordered domains, the water content in the coacervate domains was calculated to increase from approximately 16.8% to 27.5% as the polymer concentration decreased from 20 to 15 wt %.

36 MATERIALS SCIENCE↗
DOE OSTI2014

An investigation of the liquid to glass transition using integral equations for the pair structure of coupled replicae

Extensive numerical solutions of the hypernetted-chain (HNC) and Rogers-Young (RY) integral equations are presented for the pair structure of a system of two coupled replicae (1 and 2) of a “soft-sphere” fluid of atoms interacting via an inverse-12 pair potential. In the limit of vanishing inter-replica coupling ε{sub 12}, both integral equations predict the existence of three branches of solutions: (1) A high temperature liquid branch (L), which extends to a supercooled regime upon cooling when the two replicae are kept at ε{sub 12} = 0 throughout; upon separating the configurational and vibrational contributions to the free energy and entropy of the L branch, the Kauzmann temperature is located where the configurational entropy vanishes. (2) Starting with an initial finite coupling ε{sub 12}, two “glass” branches G{sub 1} and G{sub 2} are found below some critical temperature, which are characterized by a strong remnant spatial inter-replica correlation upon taking the limit ε{sub 12} → 0. Branch G{sub 2} is characterized by an increasing overlap order parameter upon cooling, and may hence be identified with the hypothetical “ideal glass” phase. Branch G{sub 1} exhibits the opposite trend of increasing order parameter upon heating; its free energy lies consistently below that of the L branch and above that of the G{sub 2} branch. The free energies of the L and G{sub 2} branches are found to intersect at an alleged “random first-order transition” (RFOT) characterized by weak discontinuities of the volume and entropy. The Kauzmann and RFOT temperatures predicted by RY differ significantly from their HNC counterparts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1998

A Geometry Based Infra-Structure for Computational Analysis and Design

The computational steps traditionally taken for most engineering analysis suites (computational fluid dynamics (CFD), structural analysis, heat transfer and etc.) are: (1) Surface Generation -- usually by employing a Computer Assisted Design (CAD) system; (2) Grid Generation -- preparing the volume for the simulation; (3) Flow Solver -- producing the results at the specified operational point; (4) Post-processing Visualization -- interactively attempting to understand the results. For structural analysis, integrated systems can be obtained from a number of commercial vendors. These vendors couple directly to a number of CAD systems and are executed from within the CAD Graphical User Interface (GUI). It should be noted that the structural analysis problem is more tractable than CFD; there are fewer mesh topologies used and the grids are not as fine (this problem space does not have the length scaling issues of fluids). For CFD, these steps have worked well in the past for simple steady-state simulations at the expense of much user interaction. The data was transmitted between phases via files. In most cases, the output from a CAD system could go to Initial Graphics Exchange Specification (IGES) or Standard Exchange Program (STEP) files. The output from Grid Generators and Solvers do not really have standards though there are a couple of file formats that can be used for a subset of the gridding (i.e. PLOT3D data formats). The user would have to patch up the data or translate from one format to another to move to the next step. Sometimes this could take days. Specifically the problems with this procedure are:(1) File based -- Information flows from one step to the next via data files with formats specified for that procedure. File standards, when they exist, are wholly inadequate. For example, geometry from CAD systems (transmitted via IGES files) is defined as disjoint surfaces and curves (as well as masses of other information of no interest for the Grid Generator). This is particularly onerous for modern CAD systems based on solid modeling. The part was a proper solid and in the translation to IGES has lost this important characteristic. STEP is another standard for CAD data that exists and supports the concept of a solid. The problem with STEP is that a solid modeling geometry kernel is required to query and manipulate the data within this type of file. (2) 'Good' Geometry. A bottleneck in getting results from a solver is the construction of proper geometry to be fed to the grid generator. With 'good' geometry a grid can be constructed in tens of minutes (even with a complex configuration) using unstructured techniques. Adroit multi-block methods are not far behind. This means that a million node steady-state solution can be computed on the order of hours (using current high performance computers) starting from this 'good' geometry. Unfortunately, the geometry usually transmitted from the CAD system is not 'good' in the grid generator sense. The grid generator needs smooth closed solid geometry. It can take a week (or more) of interaction with the CAD output (sometimes by hand) before the process can begin. One way Communication. (3) One-way Communication -- All information travels on from one phase to the next. This makes procedures like node adaptation difficult when attempting to add or move nodes that sit on bounding surfaces (when the actual surface data has been lost after the grid generation phase). Until this process can be automated, more complex problems such as multi-disciplinary analysis or using the above procedure for design becomes prohibitive. There is also no way to easily deal with this system in a modular manner. One can only replace the grid generator, for example, if the software reads and writes the same files. Instead of the serial approach to analysis as described above, CAPRI takes a geometry centric approach. This makes the actual geometry (not a discretized version) accessible to all phases of the analysis. The connection to the geometry is made through an Application Programming Interface (API) and NOT a file system. This API isolates the top-level applications (grid generators, solvers and visualization components) from the geometry engine. Also this allows the replacement of one geometry kernel with another, without effecting these top-level applications. For example, if UniGraphics is used as the CAD package then Parasolid (UG's own geometry engine) can be used for all geometric queries so that no solid geometry information is lost in a translation. This is much better than STEP because when the data is queried, the same software is executed as used in the CAD system. Therefore, one analyzes the exact part that is in the CAD system. CAPRI uses the same idea as the commercial structural analysis codes but does not specify control. Software components of the CAD system are used, but the analysis suite, not the CAD operator, specifies the control of the software session. This also means that the license issues (may be) minimized and individuals need not have to know how to operate a CAD system in order to run the suite.

Haimes, Robert↗
NASA NTRS1992

Coupled finite-difference/finite-element approach for wing-body aeroelasticity

Computational methods using finite-difference approaches for fluids and finite-element approaches for structures have individually advanced to solve almost full-aircraft configurations. However, coupled approaches to solve fluid/structural interaction problems are still in their early stages of development, particularly for complex geometries using complete equations such as the Euler/Navier-Stokes equations. Earlier work demonstrated the success of coupling finite-difference and finite-element methods for simple wing configurations using the Euler/Navier-Stokes equations. In this paper, the same approach is extended for general wing-body configurations. The structural properties are represented by beam-type finite elements. The flow is modeled using the Euler/Navier-Stokes equations. A general procedure to fully couple structural finite-element boundary conditions with fluid finite-difference boundary conditions is developed for wing-body configurations. Computations are made using moving grids that adapt to wing-body structural deformations. Results are illustrated for a typical wing-body configuration.

Guruswamy, Guru P.↗
DOE OSTI2023

In-situ experiment reveals CO2 enriched fluid migration in faulted caprock

Abstract The sealing characteristics of the geological formation located above a CO 2 storage reservoir, the so-called caprock, are essential to ensure efficient geological carbon storage. If CO 2 were to leak through the caprock, temporal changes in fluid geochemistry can reveal fundamental information on migration mechanisms and induced fluid–rock interactions. Here, we present the results from a unique in-situ injection experiment, where CO 2 -enriched fluid was continuously injected in a faulted caprock analogue. Our results show that the CO 2 migration follows complex pathways within the fault structure. The joint analysis of noble gases, ion concentrations and carbon isotopes allow us to quantify mixing between injected CO 2 -enriched fluid and resident formation water and to describe the temporal evolution of water–rock interaction processes. The results presented here are a crucial complement to the geophysical monitoring at the fracture scale highlighting a unique migration of CO 2 in fault zones.

Science & Technology - Other Topics↗
NASA NTRS1989

The critical role of aerodynamic heating effects in the design of hypersonic vehicles

Hypersonic vehicles operate in a hostile aerothermal environment, which has a significant impact on their aerothermostructural performance. Significant coupling occurs between the aerodynamic flow field, structural heat transfer, and structural response, creating a multidisciplinary interaction. The critical role of aerodynamic heating effects in the design of hypersonic vehicles is identified with an example of high localized heating on an engine-cowl leading edge. Recent advances is integrated fluid-thermal-structural finite-element analyses are presented.

Wieting, Allan R.↗
DOE OSTI2021

Microfracture-pore structure characterization and water-rock interaction in three lithofacies of the Lower Eagle Ford Formation

The microfracture-pore structures and water-rock interactions of eight samples of three typical lithofacies (three limestone, three wackestone, and two mudstone) from the Lower Eagle Ford Formation have been investigated using an integrated methodology. The methods used were X-ray diffraction (XRD), total organic carbon (TOC) content, pyrolysis, thin-section petrography, scanning electron microscopy (SEM), mercury intrusion porosimetry (MIP), spontaneous water imbibition, and water vapor & nitrogen physisorption. Petrophysical properties such as organic richness, maturity, mineral composition, types and patterns of microfractures and pores, porosity, pore-throat size distribution, pore size distribution, pore connectivity, water adsorption and desorption behavior were determined from these methods, and the results compared among the three lithofacies. For example, the three limestone EAS samples and three wackestone DM samples show a low porosity range from 2.02–3.35%, whereas the two mudstone VRU samples exhibit a much large average porosity of 11.7% by MIP results. Nitrogen physisorption results show that pore size distributions have a good similarity in N2 adsorption volume for the samples with similar lithology. Microfractures are morphologically divided into seven types, with their definitions & pictorial examples presented and possible formation mechanisms proposed. In addition, the water vapor physisorption behavior and their controlling factors are discussed. Here, the water adsorption capacity shows a positive relationship (R 2 = 0.90) with clay content and a negative relationship (R 2 = 0.99) with carbonate content. In general, this integrated study of microfracture-pore structure characterization and water-rock interaction provides some insights of water vapor adsorption behavior in different lithofacies and whether microfractures could be an important pathway for fluid flow.

58 GEOSCIENCES↗
DOE OSTI2019

Size-dependent thermodynamic structural selection in colloidal crystallization

Nucleation and growth of crystalline phases play an important role in a variety of physical phenomena, ranging from freezing of liquids to assembly of colloidal particles. Understanding these processes in the context of colloidal crystallization is of great importance for predicting and controlling the structures produced. In many systems, crystallites that nucleate have structures differing from those expected from bulk equilibrium thermodynamic considerations, and this is often attributed to kinetic effects. In this work, we consider the self-assembly of a binary mixture of colloids in two dimensions, which exhibits a structural transformation from a non–close-packed to a close-packed lattice during crystal growth. We show that this transformation is thermodynamically driven, resulting from size dependence of the relative free energy between the two structures. We demonstrate that structural selection can be entirely thermodynamic, in contrast to previously considered effects involving growth kinetics or interaction with the surrounding fluid phase.

42 ENGINEERING↗
DOE OSTI2025

A geometric perspective on kinetic matter–radiation interaction and moment systems

Here, in this paper, we provide a geometric perspective on the kinetic interaction of matter and radiation, based on a pair bracket approach. We discuss the interaction of kinetic theories via dissipative brackets, with our fundamental example being the coupling of matter, described by the Boltzmann equation, and radiation, described by the radiation transport equation. We explore the transition from kinetic systems to their corresponding moment systems, provide a Hamiltonian description of such moment systems, and give a geometric interpretation of the moment closure problem for kinetic theories. As an application, we discuss in detail diffusion radiation hydrodynamics as an example of a pair bracket formulation on a space of moments corresponding to kinetic matter–radiation interaction. Additionally, using the variable moment closure framework of Burby [J. W. Burby, Variable-moment fluid closures with Hamiltonian structure, Sci. Rep. 13 (2023) 18286, doi:10.1038/s41598-023-45416-5], we show how to construct Hamiltonian moment closures for kinetic transport equations with arbitrary Hamiltonian. Using this general construction, we derive novel Hamiltonian moment closures for pure radiation transport.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1990

Simulation of unsteady rotational flow over propfan configuration

During the past decade, aircraft engine manufacturers and scientists at NASA have worked on extending the high propulsive efficiency of a classical propeller to higher cruise Mach numbers. The resulting configurations use highly swept twisted and very thin blades to delay the drag divergence Mach number. Unfortunately, these blades are also susceptible to aeroelastic instabilities. This was observed for some advanced propeller configurations in wind tunnel tests at NASA Lewis Research Center, where the blades fluttered at cruise speeds. To address this problem and to understand the flow phenomena and the solid fluid interaction involved, a research effort was initiated at Georgia Institute of Technology in 1986, under the support of the Structural Dynamics Branch of the NASA Lewis Research Center. The objectives of this study are: (1) the development of solution procedures and computer codes capable of predicting the aeroelastic characteristics of modern single and counter-rotation propellers; and (2) the use of these solution procedures to understand physical phenomena such as stall flutter, transonic flutter, and divergence.

Srivastava, Rakesh↗
NASA NTRS1993

Steady state and dynamical structure of a cosmic-ray-modified termination shock

A hydrodynamic model is developed for the structure of a cosmic-ray-modified termination shock. The model is based on the two-fluid equations of diffuse shock acceleration (Drury and Volk, 1981). Both the steady state structure of the shock and its interaction with outer heliospheric disturbances are considered. Under the assumption that the solar wind is decelerated by diffusing interstellar cosmic rates, it is shown that the natural state of the termination shock is a gradual deceleration and compression, followed by a discontinuous jump to a downstream state which is dominated by the pressure contribution of the cosmic rays. A representative model is calculated for the steady state which incorporates both interstellar cosmic ray mediation and diffusively accelerated anomalous ions through a proposed thermal leakage mechanism. The interaction of large-scale disturbances with the equilibrium termination shock model is shown to result in some unusual downstream structure, including transmitted shocks and cosmic-ray-modified contact discontinuities. The structure observed may be connected to the 2-kHz outer heliospheric radio emission (Cairns et al., 1992a, b). The time-dependent simulations also demonstrate that interaction with solar wind compressible turbulence (e.g., traveling interplanetary shocks, etc.) could induce the termination shock to continually fluctuate between cosmic-ray-dominated and gas-dynamic states. This fluctuation may represent a partial explanation of the galactic cosmic ray modulation effect and illustrates that the Pioneer and Voyager satellites will encounter an evolving shock whose structure and dynamic properties are strongly influence by the mediation of interstellar and anomalous cosmic rays.

Donohue, D. J.↗
NASA NTRS1999

Dynamics of Single Chains of Suspended Ferrofluid Particles

We present an experimental study of the dynamics of isolated chains made of super-paramagnetic particles under the influence of a magnetic field. The motivation of this work is to understand if the chain fluctuations exist and, if it does, how does the fluctuation affect chain aggregation. We find that single chains strongly fluctuate and that the characteristic frequency of their fluctuations is inversely proportional to the magnetic field strength. The higher the field the lower the characteristic frequency of the chain fluctuations. In the high magnetic field limit, chains behave like rigid rods without any internal motions. In this work, we used ferrofluid particles suspended in water. These particles do not have any intrinsic magnetization. Once a magnetic field is applied, a dipole moment is induced in each particle, proportional to the magnetic field. A dipolar magnetic interaction then occurs between particles. If dipole-dipole magnetic energy is higher than the thermal energy, the result is a structure change inside the dipolar fluid. The ratio of these two energies is expressed by a coupling constant lambda as: lambda = (pi(a(exp 3))(chi(exp 2))(mu(sub 0))(H(sub 0))(exp 2))/18kT Where a is the particle radius, mu(sub 0) is the vacuum magnetic permeability, H(sub 0) the applied magnetic field, k the Boltzmann constant and T the absolute temperature. If lambda > 1, magnetic particles form chains along the field direction. The lateral coalescence of several chains may form bigger aggregates especially if the particle volume fraction is high. While many studies and applications deal with the rheological properties and the structural changes of these dipolar fluids, this work focuses on the understanding of the chain dynamics. In order to probe the chain dynamics, we used dynamic light scattering (DLS) in self-beating mode as our experimental technique. The experimental geometry is such that the scattering plane is perpendicular to the magnetic field. Therefore, only motions in this plane are probed. A very dilute sample of a ferrofluid emulsion with a particle volume fraction of 10(exp -5) is used in this experiment. We chose such a low volume fraction to avoid multiple light scattering as well as lateral chain-chain aggregation. DLS measures the dynamic structure factor S(q,t) of the sample (q is the scattering wave vector, t is the time). In the absence of the magnetic field, identical particles of ferrofluid droplets are randomly distributed and S(q,t) reduces to exp(-q(exp 2)2D(sub 0)t). D(sub 0)=(kT/(6(pi)(eta)(a)) is the diffusion coefficient of Brownian particles (where Xi = (6(pi)(eta)(a)) is the Stokes frictional coefficient of a spherical particle in a fluid of viscosity eta). If interactions or polydispersity can not be ignored, an effective diffusion coefficient is introduced. Formally, D(sub eff) is defined as: D(sub eff) = - q(exp -2) partial derivative of (ln(S(q,t)) with respect to time, as t goes to 0. D(sub eff) reduces to D(sub 0) if no interactions and only a few particles size are present. Therefore, we can use DLS to measure particle size. The particle radius was found to be a=0.23 mu m with 7% of polydispersity. In this case, if we vary the scattering angle theta (and so q) we do not have any change in the measured diffusion coefficient: it is q-independent. When a magnetic field is applied, particles aggregate into chains if lambda > 1. We first studied the kinetics of the chain formation when lambda = 406. At a fixed scattering angle, we measured diffusion coefficient D(sub eff) as a function of time. Experimentally, we find that D(sub eff) decreases monotonously with time. Physically, this means that chains are becoming longer and longer. Since we are only sensitive to motions in the scattering plane and since chains have their main axis perpendicular to this plane, the measured diffusion coefficient is the trans-verse diffusion coefficient. We can relate D(sub eff) to the mean number of particles per chain N(t) at a given time and to the diffusion coefficient of an isolated particle D(sub 0) as D(sub eff)=f(N(t))D(sub 0). Since f(N) is known from other recent work, N can be expressed as a function of the time. We found a square root dependency: N(t) proportional to the square root of t. As expected for very low volume fraction, this behavior is characteristic of a diffusion-limited aggregation as suggested by several authors and by our previous work. In this study, we focus on the dependence of the effective diffusion coefficient on the scattering angle and the magnetic field strength. After the magnetic field is applied (lambda = 406) for a long time, typically 6 hours, kinetics of chain formation becomes very slow. Chain size does not vary much over the next hour period. Thus, we can perform different interesting experiments. First, at a fixed magnetic field, we measure the effective diffusion coefficient as a function of the scattering angle (from 5 to 130 deg). Our results show that the measured diffusion coefficient increases linearly with the scattering angle: D(sub eff) proportional to q. If we do the same experiment for different lambda values, D(sub eff) depends on lambda as D(sub eff) proportional to lambda(exp -1/2). We also find for different lambda values that the same asymptotic D(sub eff) value is obtained when q approaches zero. The angle dependency of D(sub eff) suggests that an additional motion exists besides chain drifting. Chain size is constant during experiment, which was verified by measuring the same diffusion coefficient at the beginning and at the end of the angle switching. If chains are rigid, D(sub eff) is independent of q. Therefore, we found that D(sub eff) not only measures the motion of the entire chain but also its internal fluctuations. These internal motions are the fluctuations of the particles in the chain. To understand the q dependency of D(sub eff), let us look at the probing length used. In our study, the characteristic length scale probed is l=2pi/q which is in the range of 0.9<l/a<20. When l is much larger than the particle radius (l/a going to infinity) we are mainly sensitive to the center of mass diffusion of the chain. D(sub eff) is then the diffusion of the entire chain and depends only on N but not on q and lambda. The value D(sub eff) allows us to obtain the number of particles per chain N. In the opposite limit (i.e. l/a < 1), we are sensitive to motions on the size of individual particles. The main contributions to the measured diffusion coefficient come from internal motions of the chains (i.e. particles' fluctuations inside the chain). We investigated also the effect of the chain size on D(sub eff). For the same value of the magnetic interaction lambda and different chain sizes, we found that the slope of D(sub eff) versus q decreases when N increases as a power law D(sub eff) proportional to N(exp -0.7). This is an important behavior because this means that the longer the chain the slower its internal fluctuations besides the whole chain motion.

Cutillas, S.↗
DOE OSTI2013

Structure factor of model bidisperse ferrofluids with relatively weak interparticle interactions

In the present manuscript we develop a theoretical approach to describe the pair correlation function of bidisperse magnetic dipolar hard- and soft-spheres. We choose bidisperse system as the first step to allow for polydispersity when studying thermodynamics of magnetic fluids. Using diagram technique we calculate the virial expansion of the pair correlation function up to the first order in density and fourth order in the dipolar strength. Even though, the radial distribution functions are extremely sensitive to the steric potential, we show that the behaviour of the isotropic centre-centre structure factor is almost indifferent to the type of the short-range repulsion. We extensively compare our theoretical results to the data of molecular dynamics simulations, which helps us to understand the range of validity of the virial expansion both on density and magnetic dipolar strength. We also investigate the influence of the granulometric composition on the height, width, and position of the structure factor first peak in order to clarify whether it is possible to extract structural information from experimentally measured small angle neutron scattering intensities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRSDate not supplied

Multiphase Computational Fluid Dynamics Simulation of Ignition Overpressure and Plume-Water Interaction in the SLS Scale Model Acoustics Test

The solid rocket boosters (SRBs) on the Space Launch System (SLS) generate powerful ignition overpressure (IOP) waves that can potentially damage the vehicle and surrounding ground structures. Based on historical Shuttle experience, the IOP will be mitigated by utilizing an IOP/Sound Suppression (SS) water system which sprays water in and around the SRB and RS-25 engine plumes. However, the complex multiphase physics of the resulting plume-water interaction can lead to undesirable effects, such as high velocity water spray onto the engine nozzles and SRB nozzle side loads. Prior to launch, insight into this environment can be achieved through scale model testing and predictive computational fluid dynamics (CFD) modeling. The Scale Model Acoustic Test (SMAT) was conducted at Marshall Space Flight Center to provide the data necessary to develop the SLS liftoff acoustic environments by using a 5% scale SLS vehicle, fully integrated with a scaled IOP/SS water system, mobile launcher (ML), and tower. NASA engineers have put forth significant effort in validating the Loci/CHEM-Real Fluids multiphase CFD solver on a wet hold down test from the SMAT series. In this paper, the IOP acoustics and plume-water interactions predicted by the CFD model are compared and validated against the experiment, enabling relatively inexpensive, rapid turnaround for multiphase launch environment analyses to support future flights and design evolutions of SLS and other launch vehicles.

John E. Soto↗