Search NASA⌕ Search

SEARCH · Search NASA

Results for “base flow”

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 217 records · Page 12

Model-Based Air Flow Path Optimization of Heat Pump Clothes Dryer

A heat pump clothes dryer (HPCD) uses a vapor compression system to dry clothes. The condenser heats air, which passes through the drum to evaporate moisture out of the clothes, and the evaporator condenses water out of the air stream. As a result, the HPCD can achieve 50% energy savings compared to a conventional electric resistance dryer. In this work we developed a physics-based, quasi-steady-state HPCD system model with detailed heat exchanger and compressor models. The system model can simulate the inherently transient HPCD drying process, to size components, and to reveal trends in key variables (e.g. compressor discharge temperature, power consumption, required drying time, etc.) The system model was calibrated using experimental data from a prototype HPCD. Air leakages, in and out, along the closed air circulation path of HPCD cause varied effects on the performance. Understanding the location, magnitude, and direction of air leakage of the heat pump clothes dryer is critical for accurately characterizing the performance and developing a high-performance design. The system model was used to reveal the impacts. In addition, model-based parametric optimizations were conducted to design the HPCD charge inventory, flow rate, air path and leak points for optimum performance.

Shen, Bo↗

Inflow/outflow boundary conditions for particle-based blood flow simulations: Application to arterial bifurcations and trees

When blood flows through a bifurcation, red blood cells (RBCs) travel into side branches at different hematocrit levels, and it is even possible that all RBCs enter into one branch only, leading to a complete separation of plasma and R- BCs. To quantify this phenomenon via particle-based mesoscopic simulations, we developed a general framework for open boundary conditions in multiphase flows that is effective even for high hematocrit levels. The inflow at the inlet is duplicated from a fully developed flow generated in a pilot simulation with periodic boundary conditions. The outflow is controlled by adaptive forces to maintain the flow rate and velocity gradient at fixed values, while the particles leaving the arteriole at the outlet are removed from the system. Upon valida- tion of this approach, we performed systematic 3D simulations to study plasma skimming in arterioles of diameters 20 to 32 microns. For a flow rate ratio 6:1 at the branches, we observed the \all-or-nothing" phenomenon with plasma only entering the low flow rate branch. We then simulated blood-plasma separation in arteriolar bifurcations with different bifurcation angles and same diameter of the daughter branches. Our simulations predict a significant increase in RBC flux through the main daughter branch as the bifurcation angle is increased. Lastly, we demonstrated the new methodology for simulating blood flow in ves- sels with multiple inlets and outlets, constructed using an angiogenesis model.

59 BASIC BIOLOGICAL SCIENCES↗

Neural network–based pore flow field prediction in porous media using super resolution

Direct pore-scale simulations of fluid flow through porous media are computationally expensive to perform for realistic systems. Previous works have demonstrated using the geometry of the microstructure of porous media to predict the velocity fields therein based on neural networks. However, such trained neural networks do not perform well for unseen porous media with a large degree of heterogeneity. In this study we propose that incorporating a coarse velocity field in the input of neural networks is an effective way to improve the prediction performance. The coarse velocity field can be simulated with a low computational cost and provides global information to regularize the ill-posedness of the learning problem, which is usually caused by the use of local geometries due to the computational resource constraints. We show that incorporating the coarse-mesh velocity field significantly improves the prediction accuracy of the fine-mesh velocity field by comparison to the prediction that relies on geometric information alone, especially for the porous medium with a large interior vuggy pore space. We also show the flexibility of training the network in using coarse velocity fields with various resolutions. The results suggest that even using coarse velocity field with a very low resolution, the predictions are still enhanced and close to the ground truths. The feasibility of the method is further demonstrated by testing the trained network on real rocks. This study highlights the merits of incorporating a coarse-mesh velocity field into the input for neural networks, which provides global, physics-based information for the model, thereby improving the model's generalization capability.

42 ENGINEERING↗

Osborne Reynolds pipe flow: Direct simulation from laminar through gradual transition to fully developed turbulence

We report that the precise dynamics of breakdown in pipe transition is a century-old unresolved problem in fluid mechanics. We demonstrate that the abruptness and mysteriousness attributed to the Osborne Reynolds pipe transition can be partially resolved with a spatially developing direct simulation that carries weakly but finitely perturbed laminar inflow through gradual rather than abrupt transition arriving at the fully developed turbulent state. Our results with this approach show during transition the energy norms of such inlet perturbations grow exponentially rather than algebraically with axial distance. When inlet disturbance is located in the core region, helical vortex filaments evolve into large-scale reverse hairpin vortices. The interaction of these reverse hairpins among themselves or with the near-wall flow when they descend to the surface from the core produces small-scale hairpin packets, which leads to breakdown. When inlet disturbance is near the wall, certain quasi-spanwise structure is stretched into a Lambda vortex, and develops into a large-scale hairpin vortex. Small-scale hairpin packets emerge near the tip region of the large-scale hairpin vortex, and subsequently grow into a turbulent spot, which is itself a local concentration of small-scale hairpin vortices. This vortex dynamics is broadly analogous to that in the boundary layer bypass transition and in the secondary instability and breakdown stage of natural transition, suggesting the possibility of a partial unification. Under parabolic base flow the friction factor overshoots Moody’s correlation. Plug base flow requires stronger inlet disturbance for transition. Finally, accuracy of the results is demonstrated by comparing with analytical solutions before breakdown, and with fully developed turbulence measurements after the completion of transition.

42 ENGINEERING↗

Time-Average Molecular Rayleigh Scattering Technique for Measurement of Velocity, Denisty, Temperature, and Turbulence Intensity in High Speed Nozzle Flows

A molecular Rayleigh scattering based flow diagnostic is developed to measure time average velocity, density, temperature, and turbulence intensity in a 25.4-mm diameter nozzle free jet facility. The spectrum of the Rayleigh scattered light is analyzed using a Fabry-Perot interferometer operated in the static imaging mode. The resulting fringe pattern containing spectral information of the scattered light is recorded using a low noise CCD camera. Nonlinear least squares analysis of the fringe pattern using a kinetic theory model of the Rayleigh scattered light provides estimates of density, velocity, temperature, and turbulence intensity of the gas flow. Resulting flow parameter estimates are presented for an axial scan of subsonic flow at Mach 0.95 for comparison with previously acquired pitot tube data, and axial scans of supersonic flow in an underexpanded screeching jet. The issues related to obtaining accurate turbulence intensity measurements using this technique are discussed.

Mielke, Amy F.↗

File-Based Data Flow in the CMS Filter Farm

During the LHC Long Shutdown 1, the CMS Data Acquisition system underwent a partial redesign to replace obsolete network equipment, use more homogeneous switching technologies, and prepare the ground for future upgrades of the detector front-ends. The software and hardware infrastructure to provide input, execute the High Level Trigger (HLT) algorithms and deal with output data transport and storage has also been redesigned to be completely file- based. This approach provides additional decoupling between the HLT algorithms and the input and output data flow. All the metadata needed for bookkeeping of the data flow and the HLT process lifetimes are also generated in the form of small “documents” using the JSON encoding, by either services in the flow of the HLT execution (for rates etc.) or watchdog processes. These “files” can remain memory-resident or be written to disk if they are to be used in another part of the system (e.g. for aggregation of output data). We discuss how this redesign improves the robustness and flexibility of the CMS DAQ and the performance of the system currently being commissioned for the LHC Run 2.

Andre, J. M.↗

Flow synthesis of iron oxide nanoparticles: using multiple precursor additions to improve size control

Controlling the size of iron oxide nanoparticles while maintaining uniformity in flow-based synthesis systems has been a great challenge in nanoparticle synthesis. Using an extended LaMer mechanism, we improve both the size and shape uniformity as compared to a conventional flow synthesis. The key to this approach is the injection of additional Fe precursor during the flow reaction, providing extra precursor during the growth process leading to a size focusing step. This approach has also allowed for the systematic variation of the nanoparticle size produced in the flow reaction. This work represent a new synthetic concept for improving the size uniformity of iron oxide nanoparticles in flow-based synthesis and has potential for application across a wide range of nanoparticle systems.

36 MATERIALS SCIENCE↗

Evaluating LSM-Based Water Budgets Over a West African Basin Assisted with a River Routing Scheme

Within the framework of the African Monsoon Multidisciplinary Analysis (AMMA) Land Surface Model Intercomparison Project phase 2 (ALMIP-2), this study evaluates the water balance simulated by the Interactions between Soil, Biosphere, and Atmosphere (ISBA) over the upper Oum River basin, in Benin, using a mesoscale river routing scheme (RRS). The RRS is based on the nonlinear Muskingum Cunge method coupled with two linear reservoirs that simulate the time delay of both surface runoff and base flow that are produced by land surface models. On the basis of the evidence of a deep water-table recharge in that region,a reservoir representing the deep-water infiltration (DWI) is introduced. The hydrological processes of the basin are simulated for the 2005-08 AMMA field campaign period during which rainfall and stream flow data were intensively collected over the study area. Optimal RRS parameter sets were determined for three optimization experiments that were performed using daily stream flow at five gauges within the basin. Results demonstrate that the RRS simulates stream flow at all gauges with relative errors varying from -22% to 3% and Nash-Sutcliffe coefficients varying from 0.62 to 0.90. DWI varies from 24% to 67% of the base flow as a function of the sub-basin. The relatively simple reservoir DWI approach is quite robust, and further improvements would likely necessitate more complex solutions (e.g., considering seasonality and soil type in ISBA); thus, such modifications are recommended for future studies. Although the evaluation shows that the simulated stream flows are generally satisfactory, further field investigations are necessary to confirm some of the model assumptions.

Multidisciplinary↗

Migrational Instabilities in Particle Suspensions

This work deals with an instability arising from the shear-induced migration of particles in dense suspensions coupled with a dependence of viscosity on particle concentration. The analysis summarized here treats the inertialess (Re = O) linear stability of homogeneous simple shear flows for a Stokesian suspension model of the type proposed by Leighton and Acrivos (1987). Depending on the importance of shear-induced migration relative to concentration-driven diffusion, this model admits short-wave instability arising from wave-vector stretching by the base flow and evolving into particle-depleted shear bands. Moreover, this instability in the time-dependent problem corresponds to loss of ellipticity in the associated static problem (Re = O, Pe = O). While the isotropic version of the Leighton-Acrivos model is found to be stable with their experimentally determined parameters for simple shear, it is known that the stable model does not give a good quantitative description of particle clustering in the core of pipe flow (Nott and Brady 1994). This leads to the conjecture that an appropriate variant on the above model could explain such clustering as a two-phase bifurcation in the base flow.

Goddard, Joe D.↗

Absolute vs Convective Instabilities and Front Propagation in Lipid Membrane Tubes

Here we analyze the stability of biological membrane tubes, with and without a base flow of lipids. Membrane dynamics are completely specified by two dimensionless numbers: the well-known Föppl-von Kármán number Γ and the recently introduced Scriven-Love number SL, respectively quantifying the base tension and base flow speed. For unstable tubes, the growth rate of a local perturbation depends only on Γ, whereas SL governs the absolute versus convective nature of the instability. Furthermore, nonlinear simulations of unstable tubes reveal an initially localized disturbance result in propagating fronts, which leave a thin atrophied tube in their wake. Depending on the value of Γ, the thin tube is connected to the unperturbed regions via oscillatory or monotonic shape transitions-reminiscent of recent experimental observations on the retraction and atrophy of axons. We elucidate our findings through a weakly nonlinear analysis, which shows membrane dynamics may be approximated by a model of the class of extended Fisher-Kolmogorov equations. Our study sheds light on the pattern selection mechanism in axonal shapes by recognizing the existence of two Lifshitz points, at which the front dynamics undergo steady-to-oscillatory bifurcations.

59 BASIC BIOLOGICAL SCIENCES↗

Chemical Kinetics of the TPS and Base Bleeding During Flight Test

The present research deals with thermal degradation of polyurethane foam (PUF) during flight test. Model of thermal decomposition was developed that accounts for polyurethane kinetics parameters extracted from thermogravimetric analyses and radial heat losses to the surrounding environment. The model predicts mass loss of foam, the temperature and kinetic of release of the exhaust gases and char as function of heat and radiation loads. When PUF is heated, urethane bond break into polyol and isocyanate. In the first stage, isocyanate pyrolyses and oxidizes. As a result, the thermo-char and oil droplets (yellow smoke) are released. In the second decomposition stage, pyrolysis and oxidization of liquid polyol occur. Next, the kinetics of chemical compound release and the information about the reactions occurring in the base area are coupled to the CFD simulations of the base flow in a single first stage motor vertically stacked vehicle configuration. The CFD simulations are performed to estimate the contribution of the hot out-gassing, chemical reactions, and char oxidation to the temperature rise of the base flow. The results of simulations are compared with the flight test data.

Osipov, Viatcheslav↗

Stability Quantification for Consensus-Based Power Flow Control between Transmission and Distribution Power Systems

With increasing integration of distributed energy resources (DERs), distribution systems (DS) with DERs are expected to provide proactive grid services. The result is that power references, or known as dispatch signals, required for DS can become faster changing than in legacy power system. Consensus-based integral controls have been proposed for the purpose of coordinating power generations of DERs in DS to match the power references. These existing controls assume the integral control signal to be sufficiently slow or constant, and ignore the potential dynamics of the integral controller when tracking more varying power references. Therefore, in this paper we present an improvement for such controls by deriving the stability condition utilizing generalized Nyquist criterion. The stability condition is quantified by a set of integral gains that guarantee stability in closed-loop system without assuming a constant integral signal. A rule-of-thumb criterion is also derived to instruct the design of the consensus topology that can provide faster convergence rate for the closed-loop system. Here, the stability and convergence improvements developed in this paper are demonstrated and verified through numerical examples.

42 ENGINEERING↗

Improving Experience Beyond the Ride: Joint Assignment and Repositioning for Rideshare

In this study, a novel joint assignment and repositioning strategy is designed, implemented, and tested in a simulated traffic environment. The strategy addresses limitations of previously constructed rideshare repositioning algorithms by integrating the consideration of rider pooling into the optimization. By considering pooling, the strategy can minimize excess vehicle movement and increase riders' use of pooled services over a zonal-flow-based repositioning strategy and a fleet operating with no repositioning at all. The new strategy can increase pooling by up to 2.3% compared to no repositioning in a large region, and up to 0.1% in a small region, without compromising overall fleet revenue. Additionally, compared to a baseline, flow-based repositioning strategy, the new strategy can reduce empty vehicle travel by up to 30% in the larger region, and 2.1% in the smaller simulated area.

Paul, Joseph↗

Higher Order Chemistry Models in the CFD Simulation of Laser-Ablated Carbon Plumes

Production of single-walled carbon nanotubes (SWNT) has taken place for a number of years and by a variety of methods such-as laser ablation, chemical vapor deposition, and arc-jet ablation. Yet, little is actually understood about the exact chemical kinetics and processes that occur in SWNT formation. In recent time, NASA Johnson Space Center has devoted a considerable effort to the experimental evaluation of the laser ablation production process for SWNT originally developed at Rice University. To fully understand the nature of the laser ablation process it is necessary to understand the development of the carbon plume dynamics within the laser ablation oven. The present work is a continuation of previous studies into the efforts to model plume dynamics using computational fluid dynamics (CFD). The ultimate goal of the work is to improve understanding of the laser ablation process, and through that improved understanding, refine the laser ablation production of SWNT. Fig. 1 shows a basic schematic of the laser-ablation oven at NASA-JSC. Construction of the facility is simple in concept. Two concentric quartz tubes of 1.5 mm thickness form the inner and outer tubes with inside diameters of 2.2 and 5.08 cm respectively. At one end of the inner tube are located two 60 Hz pulsed lasers operating at 1064 nm and 532 nm wavelength with beam diameters of 5 mm aligned coaxially with the longitudinal axis of the inner quartz tube. For standard nanotube production runs, a 10 ns 532 nm pulse is followed 50 ns later by a 10 ns 1064 nm pulse. Each pulse is of 300 mJ energy. A target of carbon graphite with approximately 1% nickel and cobalt catalysts is located at the other end of the inner quartz tube. In the ordinary processing of SWNT, a base flow of 100 sccm of argon is maintained from the laser location and exits past the carbon target at a pressure of 66.7 kPa. These conditions yield a baseline mass flow through the chamber of 2.723x10(exp -6)kg/s of argon. The whole oven facility is heated to a temperature of 1473 K prior to nanotube production runs. Upon laser irradiation, part of the carbon target ablates immediately and forms a carbon vapor plume that penetrates into the argon base flow towards the laser initially at supersonic velocities. In the previous studies either a single carbon species, C3, was used to model the plume development, or a simplified 11 species carbon model reduced down from the Krestinin and Moravsky model for full fullerene chemistry was used. While both of these options yielded significant results, it was felt that the actual chemistry occ&g in the carbon plume might have a greater affect on the plume than assumed. Indeed, in the earlier attempts at modeling the carbon plume, several thermophysical characteristics could never be matched to experimental observations of plume development - mainly the propagation distance of the plume itself. In the present study, two additional chemistry models will be used to duplicate the previous studies simulations of the carbon plume. The first chemistry model used in this study is again a reduced form of the Krestinin and Moravsky rates. However the highest order carbon species allowed has been increased from C6 to C30 - therefore allowing the simulation of up to a half of the standard C60 fullerene. The second chemistry model investigated is a reduced form of a full carbon nanotube model developed at NASA-JSC. The C30 studies have already been accomplished at the present time, and the reduced SWNT model studies are currently underway. To pursue the current study, one sacrifice had to be made in that the simulation grid spacing had to be increased from 0.5 mm spacing to 1 mm spacing for the sake of computational efficiency since computational effort is proportional to the square of the number of grid points multiplied by the number of species considered. propagation that is far more in line with the experimental results observed by Puretzky et al as shown in Fig. 2. e C6 studies had yielded a far greater propagation in previous studies. In addition, chemical species development with the C30 model indicates that many higher order carbon species are produced outside of the plume proper (indicated by plotting contours of the background argon concentrations in Fig. 3) - this result was not observed in previous studies. In fact, some species primarily occurred outside of the plume itself - as shown for C27 in Fig. 4 when it is compared to Fig. 3. It could be asserted that this has occurred because all of the C27 in the plume had already been consumed in the formation of C30, but this does not seem to be indicated over time. Several other factors that arose in the previous studies have also been made more clear by the use of the higher order chemical models - one being that the use of c6 as an indicator species was mistaken. C6 is the only carbon species in the previous studies that was not injected into the flowfield as a boundary condition; it was therefore hoped that this species would provide insight into the formation of higher order carbon species for comparison to full SWNT production. But, when the plot of total mass in the plume is examined on a species by species basis in Fig. 5, it is seen that Cg was a fairly insignificant contributor to the total carbon mass in the plume and would not provide information on higher order carbon formation. the thermophysical characteristics of the carbon plume as well as simulate the carbon plume using the reduced SWNT model to provide an even better simulation of full chemistry effects upon plume propagation.

Scott, C. D.↗

Composite separators and redox flow batteries based on porous separators

Composite separators having a porous structure and including acid-stable, hydrophilic, inorganic particles enmeshed in a substantially fully fluorinated polyolefin matrix can be utilized in a number of applications. The inorganic particles can provide hydrophilic characteristics. The pores of the separator result in good selectivity and electrical conductivity. The fluorinated polymeric backbone can result in high chemical stability. Accordingly, one application of the composite separators is in redox flow batteries as low cost membranes. In such applications, the composite separator can also enable additional property-enhancing features compared to ion-exchange membranes. For example, simple capacity control can be achieved through hydraulic pressure by balancing the volumes of electrolyte on each side of the separator. While a porous separator can also allow for volume and pressure regulation, in RFBs that utilize corrosive and/or oxidizing compounds, the composite separators described herein are preferable for their robustness in the presence of such compounds.

25 ENERGY STORAGE↗

Development of a Water Based, Critical Flow, Non-Vapor Compression cooling Cycle

Expansion of a high-pressure liquid refrigerant through the use of a thermostatic expansion valve or other device is commonplace in vapor-compression cycles to regulate the quality and flow rate of the refrigerant entering the evaporator. In vapor-compression systems, as the condensed refrigerant undergoes this expansion, its pressure and temperature drop, and part of the liquid evaporates. We (researchers at Kansas State University) are developing a cooling cycle that instead pumps a high-pressure refrigerant through a supersonic converging-diverging nozzle. As the liquid refrigerant passes through the nozzle, its velocity reaches supersonic (or critical-flow) conditions, substantially decreasing the refrigerant’s pressure. This sharp pressure change vaporizes some of the refrigerant and absorbs heat from the surrounding conditions during this phase change. Due to the design of the nozzle, a shockwave trips the supersonic two-phase refrigerant back to the starting conditions, condensing the remaining vapor. The critical-flow refrigeration cycle would provide space cooling, similar to a chiller, by running a secondary fluid such as water or glycol over one or more nozzles. Rather than utilizing a compressor to raise the pressure of the refrigerant, as in a vapor-cycle system, the critical-flow cycle utilizes a high-pressure pump to drive refrigerant liquid through the cooling cycle. Additionally, the design of the nozzle can be tailored for a given refrigerant, such that environmentally benign substances can act as the working fluid. This refrigeration cycle is still in early-stage development with prototype development several years away. The complex multi-phase flow at supersonic conditions presents numerous challenges to fully understanding and modeling the cycle. With the support of DOE and venture-capital investors, initial research was conducted at PAX Streamline, and later, at Caitin. We (researchers at Kansas State University) have continued development of the cycle and have gained an in-depth understanding of the governing fundamental knowledge, based on the laws of physics and thermodynamics and verified with our testing results. Through this research, we are identifying optimal working fluid and operating conditions to eventually demonstrate the core technology for space cooling or other applications.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Validating and improving CT ventilation imaging by correlating with ventilation 4D-PET/CT using {sup 68}Ga-labeled nanoparticles

Purpose: CT ventilation imaging is a novel functional lung imaging modality based on deformable image registration. The authors present the first validation study of CT ventilation using positron emission tomography with{sup 68}Ga-labeled nanoparticles (PET-Galligas). The authors quantify this agreement for different CT ventilation metrics and PET reconstruction parameters. Methods: PET-Galligas ventilation scans were acquired for 12 lung cancer patients using a four-dimensional (4D) PET/CT scanner. CT ventilation images were then produced by applying B-spline deformable image registration between the respiratory correlated phases of the 4D-CT. The authors test four ventilation metrics, two existing and two modified. The two existing metrics model mechanical ventilation (alveolar air-flow) based on Hounsfield unit (HU) change (V{sub HU}) or Jacobian determinant of deformation (V{sub Jac}). The two modified metrics incorporate a voxel-wise tissue-density scaling (ρV{sub HU} and ρV{sub Jac}) and were hypothesized to better model the physiological ventilation. In order to assess the impact of PET image quality, comparisons were performed using both standard and respiratory-gated PET images with the former exhibiting better signal. Different median filtering kernels (σ{sub m} = 0 or 3 mm) were also applied to all images. As in previous studies, similarity metrics included the Spearman correlation coefficient r within the segmented lung volumes, and Dice coefficient d{sub 20} for the (0 − 20)th functional percentile volumes. Results: The best agreement between CT and PET ventilation was obtained comparing standard PET images to the density-scaled HU metric (ρV{sub HU}) with σ{sub m} = 3 mm. This leads to correlation values in the ranges 0.22 ⩽ r ⩽ 0.76 and 0.38 ⩽ d{sub 20} ⩽ 0.68, with r{sup ¯}=0.42±0.16 and d{sup ¯}{sub 20}=0.52±0.09 averaged over the 12 patients. Compared to Jacobian-based metrics, HU-based metrics lead to statistically significant improvements in r{sup ¯} and d{sup ¯}{sub 20} (p < 0.05), with density scaled metrics also showing higher r{sup ¯} than for unscaled versions (p < 0.02). r{sup ¯} and d{sup ¯}{sub 20} were also sensitive to image quality, with statistically significant improvements using standard (as opposed to gated) PET images and with application of median filtering. Conclusions: The use of modified CT ventilation metrics, in conjunction with PET-Galligas and careful application of image filtering has resulted in improved correlation compared to earlier studies using nuclear medicine ventilation. However, CT ventilation and PET-Galligas do not always provide the same functional information. The authors have demonstrated that the agreement can improve for CT ventilation metrics incorporating a tissue density scaling, and also with increasing PET image quality. CT ventilation imaging has clear potential for imaging regional air volume change in the lung, and further development is warranted.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

A Geometric Volume of Fluid-Based Multiphase Flow Solver Extension to the Reacting Flow Solver, PeleLM

A new algorithm is presented to simulate multiphase flows with surface tension in a pathway for spray combustion simulation. The algorithm combines capabilities from two open-source packages, including the interface reconstruction library (IRL), a library of computational geometry routines to enable the volume of fluid (VOF) method, and PeleLM, a solver for the reacting Navier-Stokes equations. Additionally, surface tension is implemented using the continuum surface force (CSF) model with an improved height function technique in the volume fraction field. Spurious errors in volume fraction arising from our combined strategy are corrected through a topology-based method that improves both numerical stability and accuracy. Multiple validation simulations are conducted, including (i) translations and rotations of Zalesak's disk, (ii) a stationary circular droplet with surface tension, (iii) an oscillating elliptical droplet, and (iv) three-dimensional deformation of a spherical droplet. Results indicate that the combined scheme retains the favorable properties of each of the component algorithms.

42 ENGINEERING↗