Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Volume of fluid”

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 469 records · Page 26

Realistic fuel spray modeling for gasoline direct injection engine applications

Fuel spray modeling plays a critical role during modern gasoline direct injection (GDI) engine development due to fuel injection’s dominant impact on engine performance and emissions as well as the complex physical processes involved. In engineering three-dimensional (3D) computational fluid dynamics (CFD) simulations, the liquid-phase fuel atomization, evaporation, and mixing are usually modeled with the discrete droplet model (DDM) adopting a Lagrangian approach for multiphase CFD simulations. To this end, general practices heavily depend on the reduced order characterization of the injector nozzle flow. However, such simplified injector modeling may lead to insufficient representations of the complex spray dynamics. To tackle this problem, this study proposes a novel workflow to numerically evaluate GDI sub-cooled and flash-boiling sprays under engine-relevant conditions using a side-mounted GDI injector together with real gasoline fuel properties. The workflow introduces a one-way coupling (OWC) method leveraging high-fidelity nozzle flow simulations to provide realistic boundary conditions to the Lagrangian injector model. The proposed workflow was first verified in a constant volume chamber (CVC) environment and then implemented in a practical GDI engine setup to study spray morphology, fuel-air mixing, and wall-wetting propensity. In addition, detailed comparison was performed between the OWC method and the conventional rate of injection (ROI) routine. In conclusion, quantitative analysis of spray characteristics was conducted to highlight possible source of discrepancies of the conventional ROI method.

33 ADVANCED PROPULSION SYSTEMS↗

Intravenous Fluid Generation System

The ability to stabilize and treat patients on exploration missions will depend on access to needed consumables. Intravenous (IV) fluids have been identified as required consumables. A review of the Space Medicine Exploration Medical Condition List (SMEMCL) lists over 400 medical conditions that could present and require treatment during ISS missions. The Intravenous Fluid Generation System (IVGEN) technology provides the scalable capability to generate IV fluids from indigenous water supplies. It meets USP (U.S. Pharmacopeia) standards. This capability was performed using potable water from the ISS; water from more extreme environments would need preconditioning. The key advantage is the ability to filter mass and volume, providing the equivalent amount of IV fluid: this is critical for remote operations or resource- poor environments. The IVGEN technology purifies drinking water, mixes it with salt, and transfers it to a suitable bag to deliver a sterile normal saline solution. Operational constraints such as mass limitations and lack of refrigeration may limit the type and volume of such fluids that can be carried onboard the spacecraft. In addition, most medical fluids have a shelf life that is shorter than some mission durations. Consequently, the objective of the IVGEN experiment was to develop, design, and validate the necessary methodology to purify spacecraft potable water into a normal saline solution, thus reducing the amount of IV fluids that are included in the launch manifest. As currently conceived, an IVGEN system for a space exploration mission would consist of an accumulator, a purifier, a mixing assembly, a salt bag, and a sterile bag. The accumulator is used to transfer a measured amount of drinking water from the spacecraft to the purifier. The purifier uses filters to separate any air bubbles that may have gotten trapped during the drinking water transfer from flowing through a high-quality deionizing cartridge that removes the impurities in the water before entering the salt bag and mixing with the salt to create a normal saline solution.

McQuillen, John↗

Pump having pistons and valves made of electroactive actuators

The present invention provides a pump for inducing a displacement of a fluid from a first medium to a second medium, including a conduit coupled to the first and second media, a transducing material piston defining a pump chamber in the conduit and being transversely displaceable for increasing a volume of the chamber to extract the fluid from the first medium to the chamber and for decreasing the chamber volume to force the fluid from the chamber to the second medium, a first transducing material valve mounted in the conduit between the piston and the first medium and being transversely displaceable from a closed position to an open position to admit the fluid to the chamber, and control means for changing a first field applied to the piston to displace the piston for changing the chamber volume and for changing a second field applied to the first valve to change the position of the first valve.

Bar-Cohen, Yoseph↗

Inkjet printing of epitaxially connected nanocrystal superlattices

Access to a blossoming library of colloidal nanomaterials provides building blocks for complex assembled materials. The journey to bring these prospects to fruition stands to benefit from the application of advanced processing methods. Epitaxially connected nanocrystal (or quantum dot) superlattices present a captivating model system for mesocrystals with intriguing emergent properties. The conventional processing approach to creating these materials involves assembling and attaching the constituent nanocrystals at the interface between two immiscible fluids. Processing small liquid volumes of the colloidal nanocrystal solution involves several complexities arising from the concurrent spreading, evaporation, assembly, and attachment. The ability of inkjet printers to deliver small (typically picoliter) liquid volumes with precise positioning is attractive to advance fundamental insights into the processing science, and thereby potentially enable new routes to incorporate the epitaxially connected superlattices into technology platforms. In this study, we identified the processing window of opportunity, including nanocrystal ink formulation and printing approach to enable delivery of colloidal nanocrystals from an inkjet nozzle onto the surface of a sessile droplet of the immiscible subphase. We demonstrate how inkjet printing can be scaled-down to enable the fabrication of epitaxially connected superlattices on patterned sub-millimeter droplets. Furthermore, we anticipate that insights from this work will spur on future advances to enable more mechanistic insights into the assembly processes and new avenues to create high-fidelity superlattices.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

SOMAFOAM: An OpenFOAM based solver for continuum simulations of low-temperature plasmas

Here, we report the development of SOMAFOAM, a finite volume framework for performing continuum simulations of low-temperature plasmas. The primary goal of this work is to discuss the features of SOMAFOAM along with representative results provided as examples for a range of operating conditions and geometries. This includes plasma and plasma–dielectric systems operating in direct current, radio frequency, and microwave regimes from pressures as low as 100 mTorr to atmospheric pressure. The code has several useful features including the ability to run massively parallel simulations using arbitrary geometries, structured/unstructured meshes, choice of models such as drift–diffusion/full-momentum at runtime, and species-dependent timesteps to name a few. The verification/validation studies presented include comparison with previously published continuum simulations (low-pressure direct current and radio frequency plasma), with experiments (Gaseous Electronics Conference Reference Cell and microwave microplasma ignited in a split ring resonator), and previously published kinetic simulations (low-pressure radio frequency plasma). Other examples provided include a direct current atmospheric pressure microplasma bounded by dielectric sidewalls and a helium–nitrogen plasma ignited using a needle electrode facing a dielectric. The performance of the code is also discussed with serial and distributed memory parallel runs demonstrated up to 512 cores. The design and implementation of the code in a modular object-oriented framework allows for easy extension and seamless coupling with other codes and can be expected to play an important role in both academia and industry.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Numerical simulation of falling film sensible heat transfer over round horizontal tubes

Numerical simulations are performed to explore the heat transfer characteristics of falling films over horizontal round tubes with uniform heat flux imposed for a range of Reynolds numbers spanning the droplet, jet, and sheet regimes. Simulations results agree well with available experimental measurements. The study analyzes the local as well as average heat transfer behavior under the different flow modes. The numerical results show that the local Nusselt number (Nu) distribution depends on the flow features in each mode and varies substantially in all directions for the respective mode. In the droplet mode, the Nu value varies significantly as the droplet impinges and the remnant liquid-bridge retracts (peak instantaneous Nu near 6), followed by wave propagation over the tube surface with peak Nu around 0.25. For the jet modes, the local maximum in the heat transfer occurs off-center to the impingement location with magnitudes of peak Nu = 3.1 for the inline jet mode and Nu = 2.7 for the staggered jet mode, while on the rest of the tube surface, it has an inverse relation with the liquid film thickness. Substantial variations in the heat transfer value are also recorded in the middle of the two impinging jets with Nu = 0.95 in the inline jet mode where the neighboring jets do not interact, and Nu = 0.60 in the crest region of the staggered jet mode where the neighboring jets interact with each other. In the sheet mode, the Nu was seen to depend on the thickness of the liquid waves traversing over the tube surface. Lower Nu values were recorded beneath the crest location of the liquid waves, which increases (1.4–11.6% depending on circumferential location) abruptly in magnitude at the advancing fronts of the waves. The temperature distribution in the liquid film in each of the modes was examined to evaluate the mechanism of heat transfer process. This study also compares the local heat transfer coefficient distribution with the analytical heat transfer models derived to predict heat transfer performance over horizontal tube surfaces.

42 ENGINEERING↗

An arbitrarily high-order three-dimensional Cartesian-grid method for reconstructing interfaces from volume fraction fields

Here Tthis work describes a newly developed, arbitrarily high-order Cartesian-grid method for reconstructing material interfaces from a volume fraction field. The method begins by identifying all of the grid cells in the volume fraction field that are intersected by the interface and need to be approximated by the reconstruction scheme. Finite-differences are used to calculate the gradient of the volume fraction field and provide an estimate of the surface normal in all of the interfacial grid cells. Groups of connected grid cells are then identified which all have the same dominant component of the normal vector. This grouping by orientation determines the proper dependent variable to use in the surface reconstruction (e.g. for a 2D curve, this step determines if the surface will be approximated by a function of x or y). A cumulative integral over the surface is constructed and fit using b-splines for two-dimensional problems or tensor-product b-splines for three-dimensional problems. This construction allows for the interface to be recovered through application of the second fundamental theorem of calculus. Fitting the cumulative integral with $\mathscr{N}$ th-order b-splines (or tensor-product b-splines) yields an ($\mathscr{N}$-1) th-order convergence rate of the interface shape. Differentiation of the b-spline interface function(s) allows for the high-order approximation of the normal vector and curvature to be obtained directly anywhere along b-spline. Together, the proposed reconstruction technique can achieve arbitrarily high mesh convergence rates. Validation tests are presented with mesh convergence rates ranging from fourth- to tenth-order.

97 MATHEMATICS AND COMPUTING↗

Identifying and tracking bubbles and drops in simulations: A toolbox for obtaining sizes, lineages, and breakup and coalescence statistics

Knowledge of bubble and drop size distributions in two-phase flows is important for characterizing a wide range of phenomena, including combustor ignition, sonar communication, and cloud formation. The physical mechanisms driving the background flow also drive the time evolution of these distributions. Accurate and robust identification and tracking algorithms for the dispersed phase are necessary to reliably measure this evolution and thereby quantify the underlying mechanisms in interface-resolving flow simulations. The identification of individual bubbles and drops traditionally relies on an algorithm used to identify connected regions. This traditional algorithm can be sensitive to the presence of spurious structures. A cost-effective refinement is proposed to maximize volume accuracy while minimizing the identification of spurious bubbles and drops. An accurate identification scheme is crucial for distinguishing bubble and drop pairs with large size ratios. The identified bubbles and drops need to be tracked in time to obtain breakup and coalescence statistics that characterize the evolution of the size distribution, including breakup and coalescence frequencies, and the probability distributions of parent and child bubble and drop sizes. An algorithm based on mass conservation is proposed to construct bubble and drop lineages using simulation snapshots that are not necessarily from consecutive time steps. These lineages are then used to detect breakup and coalescence events, and obtain the desired statistics. Accurate identification of large-size-ratio bubble and drop pairs enables accurate detection of breakup and coalescence events over a large size range. Accurate detection of successive breakup and coalescence events requires that the snapshot interval be an order of magnitude smaller than the characteristic breakup and coalescence times to capture these successive events while minimizing the identification of repeated confounding events. Together, these algorithms serve as a toolbox for detailed analysis of two-phase simulations, and enable insights into the mechanisms behind bubble and drop formation and evolution in flows of practical importance.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Cost-efficient finite-volume high-order schemes for compressible magnetohydrodynamics

We present an efficient dimension-by-dimension finite-volume method which solves the adiabatic magnetohydrodynamics equations at high discretization order, using the constrained-transport approach on Cartesian grids. Results are presented up to tenth order of accuracy. The algorithmic architecture of this method is very close to that of commonly employed second-order schemes: it requires only one reconstructed value per face for each computational cell, independently of the scheme's order. This property is highly beneficial for the numerical efficiency. It results from reusing the required values already available in neighboring grid cells, in contrast to standard algorithms that require a number of reconstructions and evaluations which increases with the scheme's order of accuracy. At a given resolution, these high-order schemes present significantly less numerical dissipation than commonly employed lower-order approaches. Thus, results of comparable accuracy are achievable at a substantially coarser resolution, yielding overall performance gains. We also present a way to include physical dissipative terms: viscosity, magnetic diffusivity and cooling functions, respecting the finite-volume and constrained-transport frameworks. Benefits of this method are shown through applications in turbulent flows.

97 MATHEMATICS AND COMPUTING↗

Go/No-Go Decision Point 1 Report

To satisfy the G/NG point 1 we have employed GEOS hydraulic fracturing module to simulate one of the Diagnostic Fracture Injection Tests conducted at well 58-32. The hydraulic fracturing solver relies on a finite element discretization of the mechanics fully coupled to a finite volume discretization of the fluid flow. Fractures are explicitly modeled as lower dimensional manifolds and discretized by 2D elements sitting at the boundary of cell elements. We have calibrated our model to reproduce the DFIT corresponding to Zone 2 Cycle 4 of the injection tests that took place at well 58-32.

15 GEOTHERMAL ENERGY↗

A novel dry coal feeding concept for high-pressure gasifiers

A novel dry coal feeding concept was developed for injecting ground coal into high-pressure gasifiers. Significant power savings are projected because the coal is injected directly with a ram and there is no requirement for pumping large volumes of gas or fluid against pressure. A novel feature of the concept is that a new seal zone is formed between the ram and injection tube each cycle. The seal zone comprises a mixture of a small quantity of finely ground coal and a fluid. To demonstrate the feasibility of the concept, coal was injected into a 1000-psi chamber with an experimental device having a 7-1/2-inch-diameter ram and a 28-inch-long stroke.

Trumbull, H. E.↗

Computer program for analysis of split-Stirling-cycle cryogenic coolers

A computer program for predicting the detailed thermodynamic performance of split-Stirling-cycle refrigerators has been developed. The mathematical model includes the refrigerator cold head, free-displacer/regenerator, gas transfer line, and provision for modeling a mechanical or thermal compressor. To allow for dynamic processes (such as aerodynamic friction and heat transfer) temperature, pressure, and mass flow rate are varied by sub-dividing the refrigerator into an appropriate number of fluid and structural control volumes. Of special importance to modeling of cryogenic coolers is the inclusion of real gas properties, and allowance for variation of thermo-physical properties such as thermal conductivities, specific heats and viscosities, with temperature and/or pressure. The resulting model, therefore, comprehensively simulates the split-cycle cooler both spatially and temporally by reflecting the effects of dynamic processes and real material properties.

Brown, M. T.↗

Computed tomography used in weld inspections at NASA

This paper describes an application of computed tomography (CT) for an evaluation of a series of welding-end hydraulic accumulators in use at Kennedy Space Center, which were suspected of possibly containing fatigue cracks in their welds, the failure which could result in a high-velocity or high-volume leakage of hydraulic fluid. The accumulator welds were evaluated in terms of incomplete penetration, porosity, and incomplete fusion. Out of 43 accumulators examined, six were given the ranking of 'good', i.e., having no incomplete penetration, and only slight porosity and incomplete fusion. These were selected for use in most critical locations of the hydraulic system used to load the Magellan spacecraft into the payload bay of the Space Shuttle Discovery, with the resulting successful loading and subsequent launching of the Magellan spacecraft.

Murray, S. H.↗

Partitioning of unstructured problems for parallel processing

Many large-scale computational problems are based on unstructured computational domains. Primary examples are unstructured grid calculations based on finite volume methods in computational fluid dynamics, or structural analysis problems based on finite element approximations. The question of how to distribute such unstructured computational domains over a large number of processors in a MIMD machine with distributed memory is addressed. A graph theoretical framework for these problems is established. Based on this framework three decomposition algorithms are introduced. In particular a new decomposition algorithm is discussed, which is based on the computation of an eigenvector of the Laplacian matrix associated with the graph. Numerical comparisons on large-scale two- and three-dimensional problems demonstrate the superiority of the new spectral bisection algorithm.

Simon, H. D.↗

Gigaflop (billion floating point operations per second) performance for computational electromagnetics

Accurate and rapid evaluation of radar signature for alternative aircraft/store configurations would be of substantial benefit in the evolution of integrated designs that meet radar cross-section (RCS) requirements across the threat spectrum. Finite-volume time domain methods offer the possibility of modeling the whole aircraft, including penetrable regions and stores, at longer wavelengths on today's gigaflop supercomputers and at typical airborne radar wavelengths on the teraflop computers of tomorrow. A structured-grid finite-volume time domain computational fluid dynamics (CFD)-based RCS code has been developed at the Rockwell Science Center, and this code incorporates modeling techniques for general radar absorbing materials and structures. Using this work as a base, the goal of the CFD-based CEM effort is to define, implement and evaluate various code development issues suitable for rapid prototype signature prediction.

Shankar, V.↗

Point to Point Control of the Hydrogen Mixer

We continue to build on the theoretical modelling of a liquid hydrogen (LH2) and gas hydrogen (GH2) mixer subsystem. The mixer described in this work is responsible for combining high pressure LH2 and GH2 to produce a hydrogen flow that meets certain thermodynamic properties. The desired properties are maintained by precise control of the LH2 and GH2 flows. The mixer is modelled as a general multi-flow lumped volume for single constituent fluids using density and internal energy as states. The set of nonlinear differential equations is modelled in the SIMULINK environment including a table look-up feature of the fluid thermodynamic properties. A small signal (linear) model is developed based on the nonlinear model and simulated as well. Pulse disturbances are introduced to the valve positions and the quality of the linear model is ascertained by comparing its behavior against the nonlinear model simulations. Valve control strategies that simulate an operator-in-the-loop scenario are then explored demonstrating the need for automatic feedback control. Finally, optimal single-output and multi-output Proportional/Integral controllers are designed based on the linear model and applied to the nonlinear model with excellent results to track simultaneous, constant setpoint changes in desired exit flow, exit temperature, and mixer pressure, as well as to reject unmeasureable but bounded additive step perturbations in the valve positions.

Barbieri, Enrique↗

Small Signal Point-to-Point Tracking of a Propellant Mixer

This paper addresses some theoretical modelling and control issues for a mixing chamber used in rocket engine testing at NASA Stennis Space Center. The mixer is responsible for combining high pressure LH2 and GH2 to produce a hydrogen flow that meets certain thermodynamic properties before it is fed into a test article. The desired properties are maintained by precise control of the LH2 and GH2 flows. The mixer is modelled as a general multi-flow lumped volume for single constituent fluids using density and internal energy as states. The set of nonlinear differential equations is modelled in the SIMULINK environment including a table look-up feature of the fluid thermodynamic properties. a small-signal (linear) model is developed based on the nonlinear model and simulated as well. Pulse disturbances are introduced to the valve positions and the quality of the linear model is ascertained by comparing its behavior against the nonlinear model simulations. Valve control strategies that simulate an operator-in-the-loop scenario are then explored demonstrating the need for automatic feedback control. Finally, classical optimal single-output and multi-output Proportional/Integral controllers are designed based on the linear model and applied to the nonlinear model with excellent results to track simultaneous, constant setpoint changes in desired exit flow, exit temperature, and mixer pressure, as well as to reject unmeasurable but bounded additive step perturbations in the valve positions.

Barbieri, Enrique↗

Small Signal Modelling and Control of the Hydrogen Mixer for Facility E1

We have undertaken the theoretical modelling of an existing liquid hydrogen (LH2) and gas hydrogen (GH2) mixer subsystem of the E1 Ground Test Facility at NASA John C. Stennis Space Center. The E1 test facility carries out comprehensive ground-based testing and certification of various liquid rocket engines and their components. The mixer described in this work is responsible for combining high pressure LH2 and GH2 to produce a hydrogen flow that meets certain thermodynamic properties before it is fed into a test article. The desired properties are maintained by precise control of the mixture of LH2 and GH2 flows. The mixer is modelled as a general multi-flow lumped volume for single constituent fluids using density and internal energy as states. The set of nonlinear differential equations is linearized about an equilibrium point and the resulting two-state, 3-input linear model is analyzed as a possible candidate for control design.

Barbieri, Enrique↗