Search NASA⌕ Search

SEARCH · Search NASA

Results for “FLOW FIELD”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Soot Predictions with a Laminar Flamelet Combustion Model in SIERRA/Fuego on a Coflow Scenario

This report describes an assessment of flamelet based soot models in a laminar ethylene coflow flame with a good selection of measurements suitable for model validation. Overall flow field and temperature predictions were in good agreement with available measurements. Soot profiles were in good agreement within the flame except for near the centerline where imperfections with the acetylene-based soot-production model are expected to be greatest. The model was challenged to predict the transition between non-sooting and sooting conditions with non-negligible soot emissions predicted even down to small flow rates or flame sizes. This suggests some possible deficiency in the soot oxidation models that might alter the amount of smoke emissions from flames, though this study cannot quantify the magnitude of the effect for large fires.

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

Turbulent entrainment in finite-length wind farms

In this article, we present an entrainment-based model for predicting the flow and power output of finite-length wind farms. The model is an extension of the three-layer approach of Luzzatto-Fegiz & Caulfield ( Phys. Rev. Fluids , vol. 3, 2018, 093802) for wind farms of infinite length, and assumes dependence of key flow quantities, such as the wind farm bulk velocity, on the streamwise distance from the farm entrance. To assist our analysis and validate the proposed model, we undertake a series of large-eddy simulations with different turbine spacing arrangements and layouts. Comparisons are also made with the top-down model with entrance effects of Meneveau ( J. Turbul. , vol. 13, 2012, N7) and data from the literature. The finite-length entrainment model is shown to be capable of capturing the power drop between contiguous rows of turbines as well as describing the advection and turbulent transport of kinetic energy in both the entrance and fully developed regions. The fully developed regime is approximated only deep in the wind farm, after approximately 15 rows of turbines. Our data suggest that for the cases considered in this study, the empirical coefficients that can be used to describe turbulent entrainment and transfers above the wind farm exhibit little dependence on the farm layout and may be considered constant for modelling purposes. However, the flow field within the wind farm layer can be strongly modulated by the turbine density (spacing) as well as the array layout, and to that extent it can be argued that they are both primary factors determining the wind farm power output.

17 WIND ENERGY↗
DOE OSTI2022

Combining shallow-water and analytical wake models for tidal array micro-siting

For tidal-stream energy to become a competitive renewable energy source, clustering multiple turbines into arrays is paramount. Array optimisation is thus critical for achieving maximum power performance and reducing cost of energy. However, ascertaining an optimal array layout is a complex problem, subject to specific site hydrodynamics and multiple inter-disciplinary constraints. In this work, we present a novel optimisation approach that combines an analytical-based wake model, FLORIS, with an ocean model, Thetis. The approach is demonstrated through applications of increasing complexity. By utilising the method of analytical wake superposition, the addition or alteration of turbine position does not require re-calculation of the entire flow field, thus allowing the use of simple heuristic techniques to perform optimisation at a fraction of the computational cost of more sophisticated methods. Using a custom condition-based placement algorithm, this methodology is applied to the Pentland Firth for arrays with turbines of 3.05 m/s rated speed, demonstrating practical implications whilst considering the temporal variability of the tide. For a 24-turbine array case, micro-siting using this technique delivered an array 15.8% more productive on average than a staggered layout, despite flow speeds regularly exceeding the rated value. Performance was evaluated through assessment of the optimised layout within the ocean model that treats turbines through a discrete turbine representation. Used iteratively, this methodology could deliver improved array configurations in a manner that accounts for local hydrodynamic effects.

16 TIDAL AND WAVE POWER↗
DOE OSTI2020

Characterization and Scaling of Forced Convective Swirl in Sinusoidal Wavy-Plate-Fin Cores of Compact Heat Exchangers

The characterization and scaling of the thermal-hydraulic performance in wavy plate-fin compact heat exchanger cores, based on the understanding of the physical phenomena and heat transfer enhancement mechanism is delineated. Experimental data are presented for forced convection in air (Pr = 0.71) with flow rates in the range 50 ≤ Re ≤ 4000. A variety of wavy-fin cores that span viable applications, with geometrical attributes described by the cross section aspect ratio α (ratio of fin spacing to height), fin corrugation aspect ratio γ (ratio of 2× corrugation amplitude to wave pitch), and fin spacing ratio ζ (ratio of fin spacing to wave pitch), are considered. To characterize and correlate the vortex-flow mixing in interfin spaces, a Swirl number Sw is introduced from the balance of viscous, inertial and centrifugal forces. It is shown that the laminar, transitional and turbulent flow regimes can be explicitly identified by this Swirl number. Based on the experimental results and extended analysis, new correlations for Fanning friction factor f and Colburn factor j are developed with Sw, α, γ, and ζ as scaling parameters. Requisite expressions are devised by a superposition of both enhancement components due to the corrugated surface area enlargement and induced swirl flow field, and they are combined to cover the laminar, transitional and turbulent regimes by the method of asymptotic matching. The resulting generalized correlations are further shown to predict all available experimental data for f and j factors to within ±20% and ±15%, respectively.

42 ENGINEERING↗
DOE OSTI2026

Background-Oriented Schlieren Velocimetry of Helium Coolant Flow in Additively Manufactured Channels

High-pressure helium gas cooling is an attractive solution for thermal management of the fusion blanket first wall, as this coolant is chemically and neutronically inert and separable from hydrogenic species. However, due to the low thermal mass of helium, geometric optimization of these channels is required to provide sufficient cooling at manageable flow rates and pumping burdens. Increasingly, analysis and optimization of these coolant channels rely on computational fluid dynamics (CFD) simulations, and these require relevant experimental data for turbulence model validation. Toward this end, a high-pressure helium gas flow visualization system has been employed to image the flow of helium in flow channels with one-sided heating, mimicking the blanket first wall environment. Flow of helium at 4 MPa pressure and flow rates up to 68 g/s (Reynolds number 57 000) is supplied to rectangular channel test sections, with uniform heating applied to the bottom wall of the channel at heat fluxes varied between roughly 50 and 130 kW/m2. A high-speed camera is used to image index of refraction (IOR) gradients in the fluid via background-oriented schlieren (BOS), and temperature and pressure instrumentation are used to characterize thermal-hydraulic performance of each channel. Cross correlation of time-resolved BOS images is then used to calculate time-averaged 2-D helium velocity fields. Flow in additively manufactured (AM) channels is examined in this manner, including both featureless channels and those containing baffling as a heat transfer enhancement. The flow distribution seen in the featureless case differs significantly from that seen in prior simulations, whereas the flow in the baffled case shows the predicted behavior of flow forced along the heated wall. This augmented flow distribution is seen to increase the heat transfer coefficient in the baffled test section. Here, strategies are discussed for ongoing and future validation of these simulations, with the aim of model deployment for blanket cooling design and optimization.

Additive manufacturing↗
DOE OSTI2022

Neural-network learning of SPOD latent dynamics

Here, we aim to reconstruct the latent space dynamics of high dimensional, quasi-stationary systems using model order reduction via the spectral proper orthogonal decomposition (SPOD). The proposed method is based on three fundamental steps: in the first, once that the mean flow field has been subtracted from the realizations (also referred to as snapshots), we compress the data from a high-dimensional representation to a lower dimensional one by constructing the SPOD latent space; in the second, we build the time-dependent coefficients by projecting the snapshots containing the fluctuations onto the SPOD basis and we learn their evolution in time with the aid of recurrent neural networks; in the third, we reconstruct the high-dimensional data from the learnt lower -dimensional representation. The proposed method is demonstrated on two different test cases, namely, a compressible jet flow, and a geophysical problem known as the Madden-Julian Oscillation. An extensive comparison between SPOD and the equivalent POD-based counterpart is provided and differences between the two approaches are highlighted. The numerical results suggest that the proposed model is able to provide low rank predictions of complex statistically stationary data and to provide insights into the evolution of phenomena characterized by specific range of frequencies. The comparison between POD and SPOD surrogate strategies highlights the need for further work on the characterization of the interplay of error between data reduction techniques and neural network forecasts.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Oscillatory and Collective Dynamics of Gold‐Nanoparticle‐Laden Droplets Driven by Photothermal‐Induced Thermocapillarity

Droplets have long intrigued researchers due to their ability to exhibit complex and fascinating behavior when subjected to external stimuli. Here, a coupled oscillatory behavior of gold-nanoparticle-surfactant-laden aqueous droplets is investigated at an oil-oil interface stimulated by light. This study shows that the interaction between light and the droplets gives rise to a range of oscillatory modes, including bouncing and clustering, where droplets exhibit collective movement. From experiment and numerical simulations, this study elucidates the underlying mechanism: upon laser irradiation, gold nanoparticles convert light into heat, generating asymmetric thermal gradients that drive upward thermocapillary flows and a hydrodynamic force from photothermal convection. These forces compete with gravity and buoyancy to induce droplet bouncing, while the resulting asymmetric flow fields bias neighboring droplets toward the illuminated droplet, leading to clustering. These findings not only expand the library of dynamic droplet behavior but also offer insights into the potential applications of light-driven systems in materials science, soft robotics, and microfluidics.

Marangoni effect↗
DOE OSTI2025

Spatially programmed alignment and actuation in printed liquid crystal elastomers

Liquid crystal elastomers (LCEs) exhibit reversible shape morphing behavior when cycled above their nematic-to-isotropic transition temperature. During extrusion-based 3D printing, LCE inks are subjected to coupled shear and extensional flows that can be harnessed to spatially control the alignment of their nematic director along prescribed print paths. Here, we combine experiment and modeling to elucidate the effects of ink composition, nozzle geometry, and printing parameters on director alignment. From rheological measurements, we quantify the dimensionless Weissenberg number ( Wi ) for the flow field each ink experiences as a function of printing conditions and demonstrate that Wi is a strong predictor of LCE alignment. We find that director alignment in LCE filaments printed through a tapered nozzle varies radially when Wi < 1, while it is uniform when Wi ≫ 1. Based on COMSOL simulations and in operando X-ray measurements, we show that LCE inks printed through nozzles with an internal hyperbolic geometry exhibit a more uniform director alignment for a given Wi compared to those through tapered nozzles. Concomitantly, the stiffness along the print direction and actuation strain of printed LCEs increases substantially under such conditions. By varying Wi during printing through adjusting the flow rate “on the fly”, LCE architectures with uniform composition, yet locally encoded shape morphing transitions can be realized.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Optical Imaging for Understanding of Thermal Barrier Coated Piston Engine Performance

Thermal barrier coatings (TBC) applied to pistons have been a recently renewed research topic in the field of internal combustion engines. Single cylinder testing of a conventional C15™ Tier 4 final pro-duction steel piston and a TBC coated piston showed that the Yttrium stabilized Zirconia (YSZ) TBC did not significantly alter engine performance, with the TBC piston having slower combustion rates and higher criteria emissions. For deeper understanding, research continued by studying 3D-piston sections optically in a constant pressure vessel using photodiodes, four high-speed cameras (natural luminosity, CH*, and OH*), and visualization of the flame from two orthogonal directions. Particle Im-aging Velocimetry (PIV) algorithms were adapted to perform Combustion Imaging Velocimetry (CIV) to aid in quantification of the observed visual combustion flow field differences. The optical work showed that the TBC piston flame development is slower, spatially different, with lower mixing energy for the TBC variant. Geometric profiles of engine pistons and optical vessel 3D-piston sections were measured using a Coordinate Measurement Machine (CMM), and surface roughness was measured with a stylus surface profilometer. CMM results showed that the TBC piston bowl had reasonable macroscopic accuracy, but the rim had a larger radius with less reentrancy. A discussion of the results and analyses in relation to future TBC effort and direction concludes the work.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Continuous caustic aqueous phase electrochemical reforming (CAPER) of ethanol for producing compressed hydrogen

Caustic Aqueous Phase Electrochemical Reforming (CAPER) of ethanol produces pure compressed H 2 at lower voltages, temperatures, and energy requirements than conventional electrolysis, while capturing carbonaceous products. The CAPER flow reactor consists of stainless steel endplates, stainless steel current collectors, and a Teflon flow field. The electrodes are prepared using commercial palladium nanoparticles on carbon support for the anode and platinum nanoparticles on carbon support for the cathode. The catalyst metal loading was 1 mgcm –2 on a 10cm 2 carbon cloth for both electrodes. Our CAPER reactor produces pure H 2 at 80bar and 80°C applying voltages at ≤0.5V. The only gas-phase product detected was H 2 and any carbonaceous products were contained in the liquid phase. The H 2 was produced with 97±4% faradaic efficiency. Further, operating at smaller electrode separation distances (4mm) and 80°C increased the current density. The highest lower heating value efficiency achieved was 28% at a flow rate of 0.05mLmin –1 and an applied voltage of 0.5V. The compressed H 2 was produced at ≤20 kWh kg –1 , which is much lower than electrolysis systems that require upwards of 50 kWh kg –1 .

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

A moving discontinuous Galerkin finite element method with interface condition enforcement for compressible flows

A variation of moving discontinuous Galerkin finite element method with interface condition enforcement (MDG-ICE) is developed for solving the compressible Euler equations. The MDG-ICE method, originating from the work of Corrigan et al. [1], [2], [3], [4], is based on the space-time DG formulation, where both flow field and grid geometry are considered as independent variables and the conservation laws are enforced both on discrete elements and element interfaces. The element conservation laws are solved in the standard discontinuous solution space to determine conservative quantities, while the interface conservation is enforced using a variational formulation in a continuous space to determine discrete grid geometry. The resulting over-determined system of nonlinear equations arising from the MDG-ICE formulation can then be solved in a least-squares sense, leading to an unconstrained nonlinear least-squares problem that is regularized and solved by Levenberg-Marquardt method. A number of numerical experiments for both 1D unsteady and 2D steady state compressible flow problems are conducted to assess the accuracy and robustness of the MDG-ICE method. Numerical results obtained indicate that the MDG-ICE method is able to implicitly detect and track all types of discontinuities via interface conservation enforcement and satisfy the conservation law on both elements and interfaces via grid movement and grid management, demonstrating that an exponential rate of convergence for Sod and Lax-Harden shock tube problems can be achieved and highly accurate solutions without overheating to both double-rarefaction wave and Noh problems can be obtained.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Recommendations on setup in simulating atmospheric gravity waves under conventionally neutral boundary layer conditions

Wind farm-induced atmospheric gravity waves have been the subject of recent research as they can impact wind farm performance. Pressure variations associated with gravity waves can contribute to the global blockage effect and wind farm wake recovery. Therefore, accurate numerical simulation of flow fields, where wind-farm-induced gravity waves may be produced, is important. Three main considerations in such simulations are the overall domain size, the use of Rayleigh damping near domain boundaries to dampen gravity waves, and advection damping at the inlet to prevent spurious oscillations. Often these considerations are treated ad hoc rather than systematically. This work aims to test and extend the systematic modelling of internal gravity waves proposed in a preliminary investigation to modelling of both internal and trapped gravity waves. The preliminary study identifies the length scales to set the domain and damping layer sizes and the time scale to configure the Rayleigh damping coefficient but under linearly stratified conditions. Large eddy simulations of flow through a wind farm canopy are performed under conventionally neutral boundary layer (CNBL) conditions to test the validity of proposed setups for CNBL conditions. Background atmospheric parameters, such as Froude number (Fr), inversion height (H i ), and inversion layer Froude number (Fr i ) control most of the atmospheric gravity wave characteristics. We validated for CBNL conditions that the effective wavelengths of the internal gravity waves are the correct length scale to configure the domain size and damping layer thickness. Likewise, the optimum damping coefficient to dampen the internal gravity waves relates to the free atmosphere's buoyancy frequency or buoyant perturbations' time scale. We infer that the damping coefficient in the inversion layer may relate to the inversion buoyancy frequency to effectively dampen the trapped gravity waves. Moreover, the advection damping length is linked to the horizontal wavelength of the trapped gravity waves in the inversion layer to prevent spurious waves at the inlet by retaining wave energy accumulation.

17 WIND ENERGY↗
DOE OSTI2023

Large Eddy Simulation of Convective Heat Transfer in a Random Pebble Bed Using the Spectral Element Method

The development of fluoride-cooled high-temperature reactors has drastically increased the demand for an in-depth understanding of the heat transfer (HT) in packed beds cooled by liquid salts. The complex flow fields and space-dependent porosity found in a pebble bed require a detailed understanding to ensure the proper cooling of the reactor core during normal and accident conditions. As detailed experimental data are complicated to obtain for these configurations, high-fidelity simulation such as large eddy simulation and direct numerical simulation (DNS) can be employed to create a high-resolution heat transfer numerical database that can assist in addressing industrial-driven issues associated with the heat transfer behavior of fluoride-cooled high-temperature reactors. In this paper, we performed a series of large eddy simulation using computational fluid dynamics (CFD) code NekRS to investigate the heat transfer for a bed of 1741 pebbles. Further, the characteristics of the flow, such as average, rms, and time series of velocity and temperature, have been analyzed. Porous media averages have also been performed. The simulation results show a good agreement between non-conjugate heat transfer and conjugate heat transfer. The generated data will be used to benchmark heat transfer modeling methods and local maxima/minima of heat transfer parameters. It will also be used for supporting convective heat transfer quantification for Kairos Power and benchmarking lower fidelity models.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Numerical Modeling of a Two-Stage Ocean Current Turbine

The Equinox Ocean Turbines (EQOT) current energy converter has a unique design with power generation in two small-diameter turbines attached to the tips of a large-diameter passive rotor. This configuration offers some key advantages for capturing ocean currents. With no centrally placed generator, almost no reaction torque is required at the nacelle of the main large-diameter rotor, and the small-diameter tip turbine generators operate at a higher speed and lower torque. The physics that determine the performance and loads on the turbine are also unique. The interactions of the flow field between the two stages and the general architecture of the system cannot be captured with traditional mid-fidelity modeling tools. For design iterations and large sets of load cases, it is important to have mid-fidelity models that can capture the important phenomenon with enough accuracy to identify global trends. This work uses a limited set of high-fidelity computational fluid dynamics (CFD) simulations to help inform the selection of and construction of a custom mid-fidelity model. Mid-fidelity modeling approaches were verified by comparing key turbine performance quantities to those found with the CFD model. Hydrodynamic interactions of the two-stage rotor were identified through high-fidelity CFD modeling. This highlighted the impact of the main rotor tip vortex and wake on the secondary rotor apparent inflow. This results in a relative flow rotation and sharp deficit, that change the optimal secondary rotor rotation speed and adds unsteadiness to the blade loading respectively. Multiple mid-fidelity approaches were evaluated for their ability to capture these effects. A simple approximation of the combined-stage performance based on single-stage BEM provides a reasonable rough prediction, especially near the peak TSR values, with some larger discrepancy at higher TSRs. Predicting the combined-stage performance based on single-stage CFD data improves this prediction across the TSR range. Although the combined-stage modeling in OLAF was not successful in this stage of the project, it showed promise as a mid-fidelity method, assuming the parameters can be tuned to account for the significant differences in time and length scales between the main and secondary rotors. This may be addressed through code changes in future work. A significant finding from the OLAF work was the agreement between the vortex core radius values found independently via a parameter space search and via CFD. The technique of using single-stage secondary rotor BEM, with a custom inflow taken from single-stage main rotor CFD or OLAF, provides an efficient method to capture one-way coupled flow interactions. This method provided generally good predictions of the impact of the flow rotation on the secondary rotor but struggled to accurately predict the peaks of the unsteady load progression. Future work could include some superposition of a tuned main rotor trailing edge viscous wake into the custom inflow to better predict this interaction.

16 TIDAL AND WAVE POWER↗
DOE OSTI2020

Analysis of combustion acoustic phenomena in compression–ignition engines using large eddy simulation

As computational capabilities continue to grow, exploring the limits of computational fluid dynamics to capture complex and elusive phenomena, which are otherwise difficult to study by experimental techniques, is one of the main targets for the research community. This paper presents a detailed analysis of the physical processes that lead to combustion noise emissions in internal combustion engines. In particular, diesel combustion in a compression-ignition (CI) engine is studied in order to understand the singular behavior of the in-cylinder flow field responsible for the acoustic emissions. The main objective is, therefore, to improve the understanding of the phenomena involved in CI engine noise using large eddy simulations. Several visualization methods are employed to investigate the connection between combustion behavior and its effects on the pressure field. In addition, proper orthogonal decomposition is used to analyze the modal energy distribution among all the acoustic modes. The results show that the acoustic signature is fundamentally conditioned by the intensity of the premixed combustion rather than by the pressure oscillations generated by turbulent fluctuations in the flame surface established during the diffusion stage.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Reproducibility of Radiokrypton in Deep Desert Aquifers: Insights from a Decade of Research

Great technical advances have been achieved since the first atom-trap trace analysis (ATTA) -based radiokrypton application in Egypt, where 1 Myr old groundwater was discovered. Beyond advances in ATTA measurement capabilities, including reduction in sample size, analysis duration, and analytical uncertainty, major progress has been achieved over the past two decades in the sample collection and preparation techniques. These advances paved the expansion of ATTA-based noble gas applications to many other aquifers worldwide, illuminating the nature and flow pattern of deep groundwater systems. While the potential of this new analytical technique for old groundwater dating is well recognized, another important aspect yet to be examined is the reproducibility of radiokrypton in aquifers over time, i.e., how representative is a discrete groundwater sample, collected at a specific time and location, for the natural groundwater system? The likelihood of a negative answer is increased by flow-field disturbance in aquifers following massive groundwater abstraction. Here, in this work, we present repeated 81 Kr sampling and measurements in twenty-one sites over Israel, mostly of deep (up to 1 km) wells tapping confined aquifers in the arid to hyperarid Negev desert. The results demonstrate that radiokrypton measurements are indeed reproducible, even in cases where samples were collected as long as nine years apart and from highly productive (∼1 Mm 3 /yr order) pumping wells. Furthermore, many of the repeated measurements in this study (17 out of the 21 sites) were conducted with different ATTA Instruments in two different laboratories using slightly different sampling, preparation, and analysis techniques, yet with an overall good agreement. The consistency in the ATTA-based 81 Kr-dating results over time highlights the robustness of this state-of-the-art technique as a tool to unravel groundwater flow patterns and encourages further applications to many other yet-to-be-explored deep aquifers.

atom-trap trace analysis↗
DOE OSTI2024

Reversible Electron-Beam Patterning of Colloidal Nanoparticles at Fluid Interfaces

The directed self-assembly of colloidal nanoparticles (NPs) using external fields guides the formation of sophisticated hierarchical materials but becomes less effective with decreasing particle size. As an alternative, electron-beam-driven assembly offers a potential avenue for targeted nanoscale manipulation, yet remains poorly controlled due to the variety and complexity of beam interaction mechanisms. Here, we investigate the beam–particle interaction of silica NPs pinned to the fluid–vacuum interface of ionic liquid droplets. In these experiments, scanning electron microscopy of the droplet surface resolves NP trajectories over space and time while simultaneously driving their reorganization. With this platform, we demonstrate the ability to direct particle transport and create transient, reversible colloidal patterns on the droplet surface. By tuning the beam voltage, we achieve precise control over both the strength and sign of the beam–particle interaction, with low voltages repelling particles and high voltages attracting them. This response stems from the formation of well-defined solvent flow fields generated from trace radiolysis of the ionic liquid, as determined through statistical analysis of single-particle trajectories under varying solvent composition. Altogether, electron-beam-guided assembly introduces a versatile strategy for nanoscale colloidal manipulation, offering new possibilities for the design of dynamic, reconfigurable systems with applications in adaptive photonics and catalysis.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Application of passive vortex generators to enhance vertical mixing in an open raceway pond

A novel use of vortex generators in a raceway pond is described here, with the flow field quantitatively simulated using computational fluid dynamics using the large eddy simulation turbulence model. Persistence lengths of the swirling motion generated by the vortex generators indicate that significant vertical mixing can be achieved by placing vortex generators in the straight section opposite the paddle wheel, downstream of the first hairpin bend. Relatively simple vortex generators are capable of creating stronger swirling motions that persist for a longer distance than those created by the paddle wheel. For optimal performance, vortex generators are positioned side by side but in opposite directions, and their diameters should be equal to or slightly less than the desired liquid depth. The optimal length of a 0.18 m diameter vortex generator in a 0.2 m deep pond was determined to be 0.3 m. As a result, it has been demonstrated that a longer persistence length is achieved by inducing a swirling motion with its rotational axis parallel to the primary flow direction.

09 BIOMASS FUELS↗