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 991 records · Page 55

Multi-physics simulations of particle tracking in arterial geometries with a scalable moving window algorithm

In arterial systems, cancer cell trajectories determine metastatic cancer locations; similarly, particle trajectories determine drug delivery distribution. Predicting trajectories is challenging, as the dynamics are affected by local interactions with red blood cells, complex hemodynamic flow structure, and downstream factors such as stenoses or blockages. Direct simulation is not possible, as a single simulation of a large arterial domain with explicit red blood cells is currently intractable on even the largest supercomputers. To overcome this limitation, we present a multi-physics adaptive window algorithm, in which individual red blood cells are explicitly modeled in a small region of interest moving through a coupled arterial fluid domain. We describe the coupling between the window and fluid domains, including automatic insertion and deletion of explicit cells and dynamic tracking of cells of interest by the window. We show that this algorithm scales efficiently on heterogeneous architectures and enables us to perform large, highly-resolved particle-tracking simulations that would otherwise be intractable.

Herschlag, Gregory↗

Isovalent Alloying and Heterovalent Substititution as Routes to Accelerate the Development and Optimization of Super-Efficient Halide Perovskite PV Solar Cells

Halide perovskites represent one of the fastest growing photovoltaic (PV) device technologies in the last decade, with research efficiencies among the highest confirmed. However, the research needed to advance tandem perovskite solar cells requires a specialized understanding base, but there is no such base for this young group of materials with the same depth as in semiconductor alloy research. Whereas IV-IV, II-V and II-VI semiconductors have a sophisticated science foundation refined over decades due of their use in electronics, lasers, LEDs and night vision applications, this is not yet the case for halide perovskites. At the same time, these materials behave rather differently than conventional semiconductors. Within the rapidly unfolding field of PV, the emerging halide-perovskite alloy research has largely been based on empiricism. Clearly, a concerted and comprehensive application to halide perovskite alloys of the modern quantum theory of solids brought us over the years a broad understanding of the critical factors enabling conventional inorganic PV could leapfrog the Edisonian, trial-and-error search for solutions to critical problems that are poorly understood. These theory approaches are required to accelerate the progress toward commercially ready hybrid-perovskites by providing a theory-guided experiment component to the EERE investments in this high-priority R&D. During this project we have studied three different rungs on the ladder of knowledge about hybrid perovskites. This kind of single halide perovskites represent the first level of development (rung 1) of solar cells, and includes materials such as APbI 3 , ASnI 3 , APbBr 3 , ASnBr 3 with different A cations, such as methylammonium (MA) with formamidinium (FM). Rung 2 is obtained by mixing B cations (Sn and Pb), and X anions (Cl and Br or I) producing a disordered alloy (A, A’)(M,M’)(X,X’) 3 -space. The most advanced rung corresponds to heterovalent-substitution of halide perovskites, creating an ordered double-perovskite compound. Although the focus of the project was on Rungs 2 and 3, we have made important discoveries regarding the nature of pure perovskites. We’ve observed that cubic perovskites should be modeled by a polymorphous network, instead of the basic single unit cell approach. This leads to the conclusion that perovskites have a fluid-like nature. An article regarding this work has been submitted to Nature Materials. Using the approach described above, we have also simulated alloys of halide perovskites, representing rung 2. We have focused on mixtures of [Cs,FA][Sn,Pb]I 3 with different concentrations, being able to disentangle the structural and electronic factors that determines the variation (bowing) of properties such as band gaps, mixing enthalpies and decomposition energies. For rung 3 compounds, we’ve closely interacted with experimental groups clarifying the chemical routes to obtain a certain compound, and what should be the most stable structure for them. As the final product, and consequent dissemination of the research, we’ve submitted two papers for publication, that are currently under evaluation. The interaction with experimental groups was also an important focus. We have participated in several local, national and international meetings, and reached out more closely to groups from CU Boulder and Arkansas.

14 SOLAR ENERGY↗

Variational nonlinear WKB in the Eulerian frame

Nonlinear WKB is a multiscale technique for studying locally plane-wave solutions of nonlinear partial differential equations (PDEs). Its application comprises two steps: (1) replacement of the original PDE with an extended system separating the large scales from the small and (2) reduction of the extended system to its slow manifold. In the context of variational fluid theories with particle relabeling symmetry, nonlinear WKB in the mean Eulerian frame is known to possess a variational structure. This much has been demonstrated using, for instance, the theoretical apparatus known as the generalized Lagrangian mean. On the other hand, the variational structure of nonlinear WKB in the conventional Eulerian frame remains mysterious. By exhibiting a variational principle for the extended equations from step (1) above, here we demonstrate that nonlinear WKB in the Eulerian frame is in fact variational. Remarkably, the variational principle for the extended system admits loops of relabeling transformations as a symmetry group. Noether’s theorem therefore implies that the extended Eulerian equations possess a family of circulation invariants parameterized by S 1 . As an illustrative example, we use our results to systematically deduce a variational model of high-frequency acoustic waves interacting with a larger-scale compressible isothermal flow.

97 MATHEMATICS AND COMPUTING↗

Hyperdiffusive dynamics in Newtonian nanoparticle fluids [Hyperdiffusive dynamics in equilibrated nanoparticle fluids]

Hyperdiffusive relaxations in soft glassy materials are typically associated with out-of-equilibrium states, and non-equilibrium physics and aging are often invoked in explaining their origins. Here, we report on hyperdiffusive motion in a model, equilibrium soft material comprised of single-component polymer-tethered-nanoparticles. In these materials, polymer mediated interactions lead to strong nanoparticle correlations, hyperdiffusive relaxations, and unusual variations of properties with temperature. Our experimental observations complement the current hypothesis that hyperdiffusive relaxations in soft materials require the material to exist in out–of–equilibrium states capable of driving structural rearrangements. Lastly, we propose alternatively that hyperdiffusive relaxations in our materials can arise naturally from volume fluctuations brought about by equilibrium thermal forces.

36 MATERIALS SCIENCE↗

Interplay of Surface Energy and Bulk Thermodynamic Forces in Ordered Block Copolymer Droplets

The wetting state of a simple liquid on a solid substrate, as summarized by Young's equation, is dictated by the interfacial energies of the different phases that coexist in the system. For simple fluids, rotational symmetry gives rise to symmetric droplets around the axis perpendicular to the substrate. This is not the case for nanostructured fluids, such as block copolymers, where the inherent thermodynamic ordering forces compete with surface tension. This competition is particularly important in nanoscale droplets, where the size of the droplets is a small multiple of the natural periodicity of the block copolymer in the bulk. In the nanoscale regime, droplet shape and internal structure arise from a subtle interplay between interfacial and bulk contributions to the free energy. In this work, we examine the consequences of surface-polymer interaction energies on droplet morphology through a concerted simulation and experimental effort. When the block copolymer is deposited on a neutral substrate we find non-circular arrangements with perpendicular domains. However, when a preferential substrate is used, the resulting morphology depends on droplet size. In large droplets, we observe bottle cap-shaped structures with a ring of perpendicular domains along the perimeter, while small droplets exhibit stripes of perpendicular domains.

36 MATERIALS SCIENCE↗

Insight into the assembly of lipid-hyaluronan complexes in osteoarthritic conditions

Interactions between molecules in the synovial fluid and the cartilage surface may play a vital role in the formation of adsorbed films that contribute to the low friction of cartilage boundary lubrication. Osteoarthritis (OA) is the most common degenerative joint disease. Previous studies have shown that in OA-diseased joints, hyaluronan (HA) not only breaks down resulting in a much lower molecular weight (MW), but also its concentration is reduced ten times. Here, we have investigated the structural changes of lipid-HA complexes as a function of HA concentration and MW to simulate the physiologically relevant conditions that exist in healthy and diseased joints. Small angle neutron scattering and dynamic light scattering were used to determine the structure of HA-lipid vesicles in bulk solution, while a combination of atomic force microscopy and quartz crystal microbalance was applied to study their assembly on a gold surface. We infer a significant influence of both MW and HA concentrations on the structure of HA-lipid complexes in bulk and assembled on a gold surface. Our results suggest that low MW HA cannot form an amorphous layer on the gold surface, which is expected to negatively impact the mechanical integrity and longevity of the boundary layer and could contribute to the increased wear of the cartilage that has been reported in joints diseased with OA.

59 BASIC BIOLOGICAL SCIENCES↗

Niowave Neutron Source Converter: Lead-Bismuth-Eutectic (LBE) Windowless Target Design and Evaluation

Los Alamos National Laboratory (LANL) is working with Niowave on the design and evaluation of their lead-bismuth-eutectic (LBE) windowless target (i.e., neutron source converter). Niowave plans to use 200 kW electron beam at 40 MeV beam energy to produce neutrons by photonuclear reaction with LBE. Then, the neutrons undergo fission at the surrounding uranium target assembly (UTA) to produce Molybdenum 99 (Mo-99) as a fission product, which eventually decays to Technetium-99 (Tc-99m). Tc-99m is one of important radioisotopes that is used for medical diagnostics. LANL conducted 3D multiphysics analysis for the Niowave neutron converter design and provided design assessment in thermohydraulic aspects. LANL conducted radiation transport calculations using Monte-Carlo N-Particle (MCNP) code with unstructured meshing scheme. The 3D volumetric heating profiles in the LBE and Stainless-Steel (SS) housing were imported into multiphase computational fluid dynamics (CFD) to obtain 3D temperature profiles of LBE and SS through conjugate heat transfer (CHT) analysis. The key findings are: LBE film thickness at the center of the beam is approximately 1.6 cm with a maximum LBE velocity of approximately 1.8 m/s, which is below a 2 m/s limit to avoid erosion issues on supporting structures; Heat deposition in the LBE peaks at ~1 cm depth from the LBE free-surface because of the forward interactions of electron, photon, and neutron with LBE; LBE maximum temperature is ~360 °C, which is below LBE evaporation initiative temperature, ~450 °C; LBE-SS interface temperature is ~350 °C, which is below the safety thermal limit to prevent severe corrosion on SS. The results indicate that Niowave’s neutron converter design satisfies both hydraulic and thermal criteria for safe operation. By virtue of such computational analysis, Niowave can move toward establishing an experimental setup to experimentally test their LBE neutron converter. The following sections describe the detailed work done by LANL.

43 PARTICLE ACCELERATORS↗

Strongly turbulent stabilization of electron beam-plasma interactions

The stabilization of electron beam interactions due to strongly turbulent nonlinearities is studied analytically and numerically for a wide range of plasma parameters. A fluid mode coupling code is described in which the effects of electron and ion Landau damping and linear growth due to the energetic electron beam are included in a phenomenological manner. Stabilization of the instability is found to occur when the amplitudes of the unstable modes exceed the threshold of the oscillating two-stream instability. The coordinate space structure of the turbulent spectrum which results clearly shows that soliton-like structures are formed by this process. Phenomenological models of both the initial stabilization and the asymptotic states are developed. Scaling laws between the beam-plasma growth rate and the fluctuations in the fields and plasma density are found in both cases, and shown to be in good agreement with the results of the simulation.

Freund, H. P.↗

NASA Lewis steady-state heat pipe code users manual

The NASA Lewis heat pipe code was developed to predict the performance of heat pipes in the steady state. The code can be used as a design tool on a personal computer or with a suitable calling routine, as a subroutine for a mainframe radiator code. A variety of wick structures, including a user input option, can be used. Heat pipes with multiple evaporators, condensers, and adiabatic sections in series and with wick structures that differ among sections can be modeled. Several working fluids can be chosen, including potassium, sodium, and lithium, for which monomer-dimer equilibrium is considered. The code incorporates a vapor flow algorithm that treats compressibility and axially varying heat input. This code facilitates the determination of heat pipe operating temperatures and heat pipe limits that may be encountered at the specified heat input and environment temperature. Data are input to the computer through a user-interactive input subroutine. Output, such as liquid and vapor pressures and temperatures, is printed at equally spaced axial positions along the pipe as determined by the user.

Tower, Leonard K.↗

OVERFLOW Validation for Predicting Plume Impingement of Underexpanded Axisymmetric Jets onto Angled Flat Plates

This report documents how OVERFLOW, a computational fluid dynamics code, predicts plume impingement of underexpanded axisymmetric jets onto both perpendicular and inclined flat plates. The effects of the plume impinging on a range of plate inclinations varying from 90deg to 30deg are investigated and compared to the experimental results in Reference 1 and 2. The flow fields are extremely complex due to the interaction between the shock waves from the free jet and those deflected by the plate. Additionally, complex mixing effects create very intricate structures in the flow. The experimental data is very limited, so these validation studies will focus only on cold plume impingement on flat and inclined plates. This validation study will help quantify the error in the OVERFLOW simulation when applied to stage separation scenarios.

Lee, Henry C.↗

Code for Multiblock CFD and Heat-Transfer Computations

The NASA Glenn Research Center General Multi-Block Navier-Stokes Convective Heat Transfer Code, Glenn-HT, has been used extensively to predict heat transfer and fluid flow for a variety of steady gas turbine engine problems. Recently, the Glenn-HT code has been completely rewritten in Fortran 90/95, a more object-oriented language that allows programmers to create code that is more modular and makes more efficient use of data structures. The new implementation takes full advantage of the capabilities of the Fortran 90/95 programming language. As a result, the Glenn-HT code now provides dynamic memory allocation, modular design, and unsteady flow capability. This allows for the heat-transfer analysis of a full turbine stage. The code has been demonstrated for an unsteady inflow condition, and gridding efforts have been initiated for a full turbine stage unsteady calculation. This analysis will be the first to simultaneously include the effects of rotation, blade interaction, film cooling, and tip clearance with recessed tip on turbine heat transfer and cooling performance. Future plans call for the application of the new Glenn-HT code to a range of gas turbine engine problems of current interest to the heat-transfer community. The new unsteady flow capability will allow researchers to predict the effect of unsteady flow phenomena upon the convective heat transfer of turbine blades and vanes. Work will also continue on the development of conjugate heat-transfer capability in the code, where simultaneous solution of convective and conductive heat-transfer domains is accomplished. Finally, advanced turbulence and fluid flow models and automatic gridding techniques are being developed that will be applied to the Glenn-HT code and solution process.

Fabian, John C.↗

Use of Solvatochromism to Assay Preferential Solvation of a Prototypic Catalytic Site

The composition of the reaction medium near photoactive catalytic sites can be inferred from the solvatochromism of the absorption and emission spectra of the wetted sites, which depend on the polarizability of the fluid. In brief, solvatochromism measures the interaction of the dipole moments of the ground and excited states with the electric field imposed by the solvent shell: a field, which does not relax on the time scale of the absorption or emission events. As a first step in establishing the utility of the technique for inorganic catalysts that operate in complex reaction media, such as encountered in the upgrading of biogenic fuels, we have measured the solvatochromism of a common, structural feature of metal oxide catalysts. In toluene, cyclohexene, chloroform and tetrahydrofuran, POSS-ligated oxometalates exhibit strong ligand-to-metal charge-transfer bands in their UV-visible absorption and emission spectra. The magnitudes of the solvatochromic effects agreed well with the calculated changes in dipole moments that accompany the electronic transitions. From the solvatochromism of the chromophores dissolved in toluene-chloroform mixtures we inferred an unexpectedly strong, preferential solvation of the chromophore even when all three components (oxometalate and the two solvents) were highly miscible. This research was supported in part by the Laboratory Directed Research & Development program at Pacific Northwest National Laboratory. PNNL is operated by Battelle for the US Department of Energy under contract DE-AC05-76RL01830. A portion of the research was performed at EMSL, a national scientific user facility sponsored by the Department of Energy's Office of Biological and Environmental Research and located at PNNL.

chromium↗

Mono- and bilayer smectic liquid crystal ordering in dense solutions of “gapped” DNA duplexes

Significance Liquid crystalline behavior of DNA has been studied for decades, yet only recently has one of the most fundamental liquid crystal phases—the smectic phase in which rod-like objects assemble into stacks of fluid layers—been demonstrated in this system. Here, we describe distinct bilayer and monolayer smectic phases in solutions of “gapped” DNA (GDNA), consisting of two DNA duplexes connected by a single strand. The smectic order is amphotropic (sensitive to temperature as well as concentration), and the layer structure is stable against significant perturbations of the GDNA construct (e.g., “gap” length variation and asymmetry in duplex lengths). The temperature sensitivity may provide a means to compare end-to-end interaction strengths between DNA duplexes at biologically relevant concentrations.

36 MATERIALS SCIENCE↗

Vidyut3d: A Non-Equilibrium Plasma Modeling Tool [SWR-24-101]

Vidyut3d is a massively-parallel plasma-fluid solver for low-temperature plasmas (LTPs) that supports both local field (LFA) and local mean energy (LMEA) approximations, as well as complex gas and surface-phase chemistry. The solver supports 2D and 3D domains, and uses AMReX's adaptive mesh refinement capabilities to increase the grid resolution around complex structures (e.g. streamer heads and sheaths) while maintaining a tractable problem size. Vidyut specializes in simulating various types of gas-phase discharges, as well as plasma/surface interactions and surface chemistry (e.g. for plasma-mediated catalysis applications). The solver also supports hybrid CPU/GPU parallelization strategies, and has demonstrated excellent scaling on various HPC architectures for problem sizes consisting of O(100 M) control volumes.

Sitaraman, Hariswaran↗

N-fold Darboux transformation and double-Wronskian-typed solitonic structures for a variable-coefficient modified Kortweg-de Vries equation

Under investigation in this paper is a variable-coefficient modified Kortweg-de Vries (vc-mKdV) model describing certain situations from the fluid mechanics, ocean dynamics and plasma physics. N-fold Darboux transformation (DT) of a variable-coefficient Ablowitz-Kaup-Newell-Segur spectral problem is constructed via a gauge transformation. Multi-solitonic solutions in terms of the double Wronskian for the vc-mKdV model are derived by the reduction of the N-fold DT. Three types of the solitonic interactions are discussed through figures: (1) Overtaking collision; (2) Head-on collision; (3) Parallel solitons. Nonlinear, dispersive and dissipative terms have the effects on the velocities of the solitonic waves while the amplitudes of the waves depend on the perturbation term. - Highlights: Black-Right-Pointing-Pointer N-fold DT is firstly applied to a vc-AKNS spectral problem. Black-Right-Pointing-Pointer Seeking a double Wronskian solution is changed into solving two systems. Black-Right-Pointing-Pointer Effects of the variable coefficients on the multi-solitonic waves are discussed in detail. Black-Right-Pointing-Pointer This work solves the problem from Yi Zhang [Ann. Phys. 323 (2008) 3059].

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Anisotropy and probe-medium interactions in the microrheology of nematic fluids.

A theoretical formalism is presented to analyze and interpret microrheology experiments in anisotropic fluids with nematic order. The predictions of that approach are examined in the context of a simple coarse-grained molecular model which is simulated using nonequilibrium molecular dynamics calculations. The proposed formalism is used to study the effect of confinement, the type of anchoring at the probe-particle surface, and the strength of the nematic field on the rheological response functions obtained from probe-particle active microrheology. As expected, a stronger nematic field leads to increased anisotropy in the rheological response of the material. It is also found that the defect structures that arise around the probe particle, which are determined by the type of anchoring and the particle size, have a significant effect on the rheological response observed in microrheology simulations. Independent estimates of the bulk dynamic modulus of the model nematic fluid considered here are obtained from small-amplitude oscillatory shear simulations with Lees Edwards boundary conditions. The results of simulations indicate that the dynamic modulus extracted from particle-probe microrheology is different from that obtained in the absence of the particle, but that the differences decrease as the size of the defect also decreases. Importantly, the results of the nematic microrheology theory proposed here are in much closer agreement with simulations than those from earlier formalisms conceived for isotropic fluids. As such, it is anticipated that the theoretical framework advanced in this study could provide a useful tool for interpretation of microrheology experiments in systems such as liquid crystals and confined macromolecular solutions or gels.

Cordoba, Andres↗

On the dynamics of a shock-bubble interaction

We present a detailed numerical study of the interaction of a weak shock wave with an isolated cylindrical gas inhomogenity. Such interactions have been studied experimentally in an attempt to elucidate the mechanisms whereby shock waves propagating through random media enhance mixing. Our study concentrates on the early phases of the interaction process which are dominated by repeated refractions of acoustic fronts at the bubble interface. Specifically, we have reproduced two of the experiments performed by Haas and Sturtevant : M(sub s) = 1.22 planar shock wave, moving through air, impinges on a cylindrical bubble which contains either helium or Refrigerant 22. These flows are modelled using the two-dimensional, compressible Euler equations for a two component fluid (air-helium or air-Refrigerant 22). Although simulations of shock wave phenomena are now fairly commonplace, they are mostly restricted to single component flows. Unfortunately, multi-component extensions of successful single component schemes often suffer from spurious oscillations which are generated at material interfaces. Here we avoid such problems by employing a novel, nonconservative shock-capturing scheme. In addition, we have utilized a sophisticated adaptive mesh refinement algorithm which enables extremely high resolution simulations to be performed relatively cheaply. Thus we have been able to reproduce numerically all the intricate mechanisms that were observed experimentally (e.g., transitions from regular to irregular refraction, cusp formation and shock wave focusing, multi-shock and Mach shock structures, jet formation, etc.), and we can now present an updated description for the dynamics of a shock-bubble interaction.

Quirk, James J.↗

Experimental measurement of the effective contact angle for solvent/packing interactions in a structured packed column for CO 2 capture

The contact angle is a critical factor for determining the effective mass transfer area for carbon dioxide (CO 2 ) capture via the chemical absorption process in a packed column, and thus the overall capture efficiency of the packed column. Many widely used commercial packings involve microscale features (perforation, corrugation, etc.) that may also affect the wetting behavior. This study proposes a systematic method of using a modified Wilhelmy plate to measure the effective contact angle to characterize the solvent and featured packing interaction. In lieu of computational efforts relying on assumptions guided by semi-informed correlation, the proposed method directly measures the effective contact angle as a function of the solvent and packing thermophysical and hydrodynamic properties. The characterization of the effective contact angle is then integrated in the computational fluid dynamics modeling to reduce uncertainty in the prediction of effective mas transfer area. Experiments were conducted for stainless steel coupons using water and aqueous sodium hydroxide (NaOH) solvent for verification of the proposed experimental protocol. The surface tension of the aqueous NaOH solvent was altered using surfactant and antifoam. The effective contact angle increases with the increased value of the surface tension for flat stainless steel sheets. On the other hand, effective contact angles do not vary in Mellapak coupons for aqueous monoethanolamine (MEA) and NaOH solvent. In this case, surface textures and sheet design play a dominant role in the surface tension of solvents. Furthermore, CO 2 loading has a significant effect on the contact angle for Mellapak coupons. As expected, the contact angle decreases with the increasing temperature of CO 2 capture solvents (MEA, EEMPA). The effective contact angle measurement using the Wilhelmy plate method can provide more accurate and efficient solvents for characterizing the solvent-packing surface interactions. Subsequently, the present method enhances the accuracy of the prediction of the effective mass transfer area in carbon capture by solvent absorption.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗