SEARCH · Engineering Papers

Results for “hydrodynamic efficiency”

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.

SourceResearch connectionFrom yearThrough year

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

At least 523 records · Page 29

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Heat transfer in directly-irradiated high-temperature solid–gas flows laden with polydisperse particles

Heat transfer in directly-irradiated high-temperature solid–gas flows laden with polydisperse particles is investigated using a novel transient three-dimensional computational fluid dynamics model. The model couples particle–gas hydrodynamics of solid–gas flows laden with polydisperse particles, radiative heat transfer in non-grey absorbing, emitting and anisotropically-scattering multi-component participating media, conduction heat transfer in the gas phase, and interfacial convection heat transfer. The multiphase particle-in-cell method is used to predict high-fidelity solid–gas flow characteristics, such as the local discrete particle size distribution, with increased computational efficiency by combining the advantages of both Eulerian and Lagrangian methods. The multi-component radiative transfer model is implemented using an advanced collision-based Monte Carlo ray-tracing method. The number of the prescribed discrete particle components is found to be the key parameter affecting the computational accuracy and efficiency, which primarily depends on the size distribution of the particles. For the model particle–gas flow featuring free-falling Gamma-distributed ceramic particles exposed to concentrated solar irradiation, the particle volume fraction, radiative, fluid flow and thermal characteristics appear to converge with the increasing number of the discrete particle components. Five particle components are sufficient to obtain physically meaningful results. A further increase in the number of the particle components only slightly increases the accuracy of the numerical predictions at the expense of a rapidly increasing computational time. For five particle components, the particle vertical velocity at the receiver exit for particles with the diameter of 43.4 μm is 57% of that for the particles with the diameter of 202.8 μm. The temperatures of these two particle components increase from the initial ambient values by factors of 2 and 1.2, respectively, during the simulation time. In conclusion, the model developed allows for increased fidelity of particle–gas flow simulations with significant radiative effects.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Parthenon—a performance portable block-structured adaptive mesh refinement framework

On the path to exascale the landscape of computer device architectures and corresponding programming models has become much more diverse. While various low-level performance portable programming models are available, support at the application level lacks behind. To address this issue, we present the performance portable block-structured adaptive mesh refinement (AMR) framework Parthenon, derived from the well-tested and widely used Athena++ astrophysical magnetohydrodynamics code, but generalized to serve as the foundation for a variety of downstream multi-physics codes. Parthenon adopts the Kokkos programming model, and provides various levels of abstractions from multidimensional variables, to packages defining and separating components, to launching of parallel compute kernels. Parthenon allocates all data in device memory to reduce data movement, supports the logical packing of variables and mesh blocks to reduce kernel launch overhead, and employs one-sided, asynchronous MPI calls to reduce communication overhead in multi-node simulations. Using a hydrodynamics miniapp, we demonstrate weak and strong scaling on various architectures including AMD and NVIDIA GPUs, Intel and AMD x86 CPUs, IBM Power9 CPUs, as well as Fujitsu A64FX CPUs. At the largest scale on Frontier (the first TOP500 exascale machine), the miniapp reaches a total of 1.7 × 10 13 zone-cycles/s on 9216 nodes (73,728 logical GPUs) at [Formula: see text] weak scaling parallel efficiency (starting from a single node). In combination with being an open, collaborative project, this makes Parthenon an ideal framework to target exascale simulations in which the downstream developers can focus on their specific application rather than on the complexity of handling massively-parallel, device-accelerated AMR.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2016

Use of Ferrihydrite-Coated Pozzolana and Biogenic Green Rust to Purify Waste Water Containing Phosphate and Nitrate

The activated sludge treatments combined to the addition of ferric chloride is commonly used to eliminate nitrate and phosphate from waste water in urban area. These processes that need costly infrastructures are not suitable for rural areas and passive treatments (lagoons, reed bed filters…) are more frequently performed. Reed bed filters are efficient for removing organic matter but are not suitable for treating phosphate and nitrate as well. Passive water treatments using various materials (hydroxyapatite, slag…) were already performed, but those allowing the elimination of both nitrate and phosphate are not actually available. The goal of this work is to identify the most suitable iron based materials for such treatments and to determine their optimal use conditions, in particular in hydrodynamic mode. The reactivity of the iron based minerals was measured either by using free particles in suspension or by depositing these particles on a solid substrate. Pouzzolana that is characterized by a porous sponge-like structure suits for settling a high amount of iron oxides. The experimental conditions enabling to avoid any ammonium formation when green rust encounters nitrate were determined within the framework of a full factorial design. The process is divided into two steps that will be performed inside two separated reactors. Indeed, the presence of phosphate inhibits the reduction of nitrate by green rust and the dephosphatation process must precede the denitrification process. In order to remove phosphate, ferrihydrite coated pouzzolana is the best materials. The kinetics of reaction of green rust with nitrate is relatively slow and often leads to the formation of ammonium. The recommendation of the identified process is to favor the accumulation of nitrite in a first step, these species reacting much more quickly with green rust and do not transform into ammonium.

Environmental Molecular Sciences Laboratory↗

DOE OSTI2015

Operational constraints and strategies for systems to effect the sustainable, solar-driven reduction of atmospheric CO 2

The operational constraints for a 6-electron/6-proton CO 2 reduction system that operates at the concentration of CO 2 in the current atmosphere (p CO 2 = 400 ppm) have been evaluated on a variety of scale lengths that span from laboratory scale to global scale. Due to the low concentration of CO 2 in the atmosphere, limitations due to mass transport of CO 2 from the tropopause have been evaluated through five different regions, each with different characteristic length scales: the troposphere; the atmospheric boundary layer (ABL); the canopy layer; a membrane layer; and an aqueous electrolyte layer. The resulting CO 2 conductances, and associated physical transport limitations, will set the ultimate limit on the efficiency and areal requirements of a sustainable solar-driven CO 2 reduction system regardless of the activity or selectivity of catalysts for reduction of CO 2 at the molecular level. At the electrolyte/electrode interface, the steady-state limiting current density and the concomitant voltage loss associated with the CO 2 concentration overpotential in a one-dimensional solar-driven CO 2 reduction cell have been assessed quantitatively using a mathematical model that accounts for diffusion, migration and convective transport, as well as for bulk electrochemical reactions in the electrolyte. At p CO 2 = 400 ppm, the low diffusion coefficient combined with the low solubility of CO 2 in aqueous solutions constrains the steady-state limiting current density to <0.1 mA cm -2 in a typical electrochemical cell with natural convection and employing electrolytes with a range of pH values. Hence, in such a system, the CO 2 capture area must be 100- to 1000-fold larger than the solar photon collection area to enable a >10% efficient solar-driven CO 2 reduction system (based on the solar collection area). This flux limitation is consistent with estimates of oceanic CO 2 uptake fluxes that have been developed in conjunction with carbon-cycle analyses for use in coupled atmosphere/ocean general circulation models. Two strategies to improve the feasibility of obtaining efficient and sustainable CO 2 transport to a cathode surface at p CO 2 = 400 ppm are described and modeled quantitatively. The first strategy employs yet unknown catalysts, analogous to carbonic anhydrases, that dramatically accelerate the chemically enhanced CO 2 transport in the aqueous electrolyte layer by enhancing the acid–base reactions in a bicarbonate buffer system. The rapid interconversion from bicarbonate to CO 2 in the presence of such catalysts near the cathode surface would in principle yield significant increases in the steady-state limiting current density and allow for >10% solar-fuel operation at the cell level. The second strategy employs a thin-layer cell architecture to improve the diffusive transport of CO 2 by use of an ultrathin polymeric membrane electrolyte. Rapid equilibration of CO 2 at the gas/electrolyte interface, and significantly enhanced diffusive fluxes of CO 2 in electrolytes, are required to increase the steady-state limiting current density of such a system. This latter approach however only is feasible for gaseous products, because liquid products would coat the electrode and therefore thicken the hydrodynamic boundary layer and accordingly reduce the diffusive CO 2 flux to the electrode surface. Regardless of whether the limitations due to mass transport to the electrode surface are overcome on the laboratory scale, at global scales the ultimate CO 2 flux limitations will be dictated by mass transport considerations related to transport of atmospheric CO 2 to the boundary plane of the solar-driven reactor system. The transport of CO 2 across the troposphere/ABL interface, the ABL/canopy layer interface, and the canopy layer/electrolyte interface have therefore been assessed in this work, to provide upper bounds on the ultimate limits for the solar-to-fuel (STF) conversion efficiency for systems that are intended to effect the reduction of atmospheric CO 2 in a sustainable fashion at global scale.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2015

On the interaction between the external magnetic field and nanofluid inside a vertical square duct

In this paper, we numerically study how the external magnetic field influences the flow and thermal characteristics of nanofluid inside a vertical square duct. The flow is considered to be laminar and hydrodynamically as well as thermally developed, whereas the thermal boundary condition of constant heat flux per unit axial length with constant peripheral temperature at any cross section, is assumed. The governing equations are solved using the spectral method and the finite difference method. Excellent comparison is noted in the numerical results given by the two methods but the spectral method is found to be superior in terms of both efficiency and accuracy. We have noted that the flow reversal due to high Raleigh number may be controlled by applying an external magnetic field of suitable strength. Moreover, the Nusselt number is found to be almost a linear function of the nanoparticle volume fraction parameter, for different values of the Raleigh number and the magnetic parameter.

36 MATERIALS SCIENCE↗

DOE OSTI2019

High-temperature hohlraum designs with multiple laser-entrance holes

Multiple laser-entrance-hole (LEH) designs are proposed which increase the number of LEHs, n, from two in standard designs. This is done to minimize the laser travel distance in the hohlraum and to obtain algebraically vanishing low order radiation moments, thereby allowing smaller case-to-capsule ratios. This leads to higher coupling efficiencies and temperatures. Symmetry is analyzed using group theory for $n$ ≤ 12 with the LEHs placed at points given by solutions to the Thomson problem. Symmetry is improved beyond the standard n = 2 designs for n corresponding to Platonic solids: the tetrahedron ($n$ = 4), octahedron ($n$ = 6), and icosahedron ($n$ = 12). The first, non-vanishing asymmetry present in the radiation drive is given. Two-dimensional radiation-hydrodynamic simulations are then performed using Lasnex to assess energetics for $n$ = 4 and $n$ = 12. The simulation domain is a conical section of a sphere with the total solid angle equal to 4$π/n$. The total LEH area is kept fixed as n increases, reducing the size of an individual LEH and the laser spot size. A foam block is recessed inside the LEH in order to capture all of the incident laser energy and prevent the rays from propagating large distances across the hohlraum. The radiation temperature near the capsule reaches ~335 eV with a 400 TW peak laser power, and the electron density inside the LEH remains below quarter-critical density.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Transfer learning driven design optimization for inertial confinement fusion

Transfer learning is a promising approach to create predictive models that incorporate simulation and experimental data into a common framework. In this technique, a neural network is first trained on a large database of simulations and then partially retrained on sparse sets of experimental data to adjust predictions to be more consistent with reality. Previously, this technique has been used to create predictive models of Omega [Humbird et al., IEEE Trans. Plasma Sci. 48, 61–70 (2019)] and NIF [Humbird et al., Phys. Plasmas 28, 042709 (2021); Kustowski et al., Mach. Learn. 3, 015035 (2022)] inertial confinement fusion (ICF) experiments that are more accurate than simulations alone. Here, in this work, we conduct a transfer learning driven hypothetical ICF campaign in which the goal is to maximize experimental neutron yield via Bayesian optimization. The transfer learning model achieves yields within 5% of the maximum achievable yield in a modest-sized design space in fewer than 20 experiments. Furthermore, we demonstrate that this method is more efficient at optimizing designs than traditional model calibration techniques commonly employed in ICF design. Such an approach to ICF design could enable robust optimization of experimental performance under uncertainty.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2022

Linear Disturbance Amplification Over Blunted Flat Plates in High-Speed Flows

Modal and nonmodal instability characteristics of cylindrically blunted flat plates with varying leading edge radii are described for Mach 4 and Mach 6 freestream conditions. The selection of leading edge radii and freestream parameters is informed by experimental conditions. The investigation of this 2D problem provides a slow entropy layer swallowing which allows for an isolated development of perturbations seeded upstream within different wall-normal regions of the flow. At both Mach numbers, a decrease in modal instability amplification was seen as the leading edge radius was increased. Nonmodal analysis reveals amplifying perturbations in the boundary layer as well as in the entropy layer. The medium bluntness regime exhibits the strongest amplification of nonmodal disturbances that is nonmonotonic in character. Small amplitude boundary forcing at the plate surface or volumetric forcing at various wall-normal heights was used to account for receptivity effects. While wall forcing effectively induced modal instabilities, only an actuation above the boundary layer captured disturbances that amplify within the entropy layer. The optimal nonmodal theory’s entropy-layer disturbance evolution exhibited outstanding agreement with controlled forcing, including receptivity effects. The evolution of entropy-layer disturbances from the optimal nonmodal theory showed excellent agreement with the results of receptivity to controlled forcing. Therefore, the nonmodal optimal growth analysis may provide a useful as well as efficient technique to identify the complete disturbance spectrum in blunt hypersonic configurations, where both modal and nonmodal disturbances can amplify in the boundary-layer and entropy-layer regions.

Boundary layer transition↗

DOE OSTI2017

Hydrodynamical Coupling of Mass and Momentum in Multiphase Galactic Winds

Using a set of high-resolution hydrodynamical simulations run with the Cholla code, we investigate how mass and momentum couple to the multiphase components of galactic winds. The simulations model the interaction between a hot wind driven by supernova explosions and a cooler, denser cloud of interstellar or circumgalactic media. By resolving scales of ${\rm{\Delta }}x\lt 0.1$ pc over $\gt 100$ pc distances, our calculations capture how the cloud disruption leads to a distribution of densities and temperatures in the resulting multiphase outflow and quantify the mass and momentum associated with each phase. We find that the multiphase wind contains comparable mass and momenta in phases over a wide range of densities and temperatures extending from the hot wind ($n\approx {10}^{-2.5}$ ${\mathrm{cm}}^{-3}$, $T\approx {10}^{6.5}$ K) to the coldest components ($n\approx {10}^{2}$ ${\mathrm{cm}}^{-3}$, $T\approx {10}^{2}$ K). We further find that the momentum distributes roughly in proportion to the mass in each phase, and the mass loading of the hot phase by the destruction of cold, dense material is an efficient process. These results provide new insight into the physical origin of observed multiphase galactic outflows and inform galaxy formation models that include coarser treatments of galactic winds. Our results confirm that cool gas observed in outflows at large distances from the galaxy ($\gtrsim 1$ kpc) likely does not originate through the entrainment of cold material near the central starburst.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2014

Bright x-ray sources from laser irradiation of foams with high concentration of Ti

Low-density foams irradiated by a 20 kJ laser at the Omega laser facility (Laboratory for Laser Energetics, Rochester, NY, USA) are shown to convert more than 5% of the laser energy into 4.6 to 6.0 keV x rays. This record efficiency with foam targets is due to novel fabrication techniques based on atomic-layer-deposition of Ti atoms on an aerogel scaffold. A Ti concentration of 33 at. % was obtained in a foam with a total density of 5 mg/cm{sup 3}. The dynamics of the ionization front through these foams were investigated at the 1 kJ laser of the Gekko XII facility (Institute for Laser Engineering, Osaka, Japan). Hydrodynamic simulations can reproduce the average electron temperature but fail to predict accurately the heat front velocity in the foam. This discrepancy is shown to be unrelated to the possible water adsorbed in the foam but could be attributed to effects of the foam micro-structure.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2017

Fast heating of fuel assembled in a spherical deuterated polystyrene shell target by counter-irradiating tailored laser pulses delivered by a HAMA 1 Hz ICF driver

Fast heating is a method of heating an assembled high-density plasma into a hot state by irradiating it with short-duration (sub-picosecond), high-intensity ($>$ $10^{18}$ W ${\rm cm}^{-2}$ ) laser pulses before the plasma expands and dissolves hydrodynamically. In this paper, we present detailed experimental results of fast heating fuel assembled in a spherical deuterated polystyrene shell target of 500 ${\mu}$ m diameter and 7 μm thickness with counterbeam illumination by using a HAMA 1 Hz, 5.9 J inertial confinement fusion laser driver with pulse tailoring. These tailored pulses contain three pulses in sequence: a 'foot' pulse of 2.4 J/25 ns, a 'spike' pulse of 0.5 J/300 ps and a 'heater' pulse of 0.4 J/110 fs; these pulses are designed to assemble the fuel and heat it. By varying the energy of the foot pulse, we find that fast heating the fuel is achieved only if the fuel is weakly ablated by the foot pulse and then shock-assembled by the spike pulse into the target centre so that the heater pulse can access the fuel with a focal intensity greater than $10^{18}$ W ${\rm cm}^{-2}$ . Without a foot pulse, the heater pulse contributes to assembling the fuel. For higher foot-pulse energies, the heater pulse drives a hydrodynamic motion with speeds of the order 107 cm ${\rm s}^{-1}$ with intensities of the order $10^{17}$ W ${\rm cm}^{-2}$ , resulting in re-assembling and additional heating of the pre-assembled fuel. Once a shock-assembled core is achieved at the target centre, we succeed qualitatively in fast heating the core for shots in sequence with variations of laser energy within 18%. The coupling efficiency from the heating laser to the core is inferred to be $(10 \pm 2)$ % in total: $(8 \pm 1.6)$ % for the ionized bulk electrons and $(2 \pm 0.4)$ % for the bulk ions. The fusion neutron spectrum detected on the laser axis exhibits peaks at 1.0 MeV, 1.7 MeV and 3.8 MeV. These peaks are attributed to the $C(d, n){\hspace{0pt}}^{13}$ N and $d(d, n){\hspace{0pt}}^3$ He reactions induced by counterpropagating fast deuterons accelerated by the photon pressure of the heating pulses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1987

The role of turbulent convection in the primitive solar nebula. I - Theory. II - Results

A model of convective turbulence which takes radiative dissipation, rotation, and convective motion anisotropy into account, on the basis of a closure for the nonlinear interactions that employs the growth rates of hydrodynamic instabilities, is used to obtain a theoretical framework for modeling the primordial solar nebula. It is assumed that convection is the sole source of turbulence causing the solar nebula to evolve. Vertical structure equations in the thin disk approximation are developed and a detailed comparison with the previous solar nebula convective models of such workers as Lin et al. (1981, 1982) is undertaken. The present values for the turbulent efficiency are much lower and more sensitive to opacity and surface density, resulting in low turbulent speeds, a more massive disk, a lower accretion rate 'best value', and a longer characteristic dispersal time for the disk. It is concluded that convection may not be the dominant source of turbulence needed to evolve young solar/stellar nebulae.

Cabot, W.↗

DOE OSTI2015

Relationship between star formation rate and black hole accretion at z=3: the different contributions in quiescent, normal, and starburst galaxies

We investigate the co-evolution of the black hole accretion rate (BHAR) and the star formation rate (SFR) in 1.5<z<2.5 galaxies displaying a greater diversity of star-forming properties compared to previous studies. We combine X-ray stacking and far-IR photometry of stellar mass-limited samples of normal star-forming, starburst, and quiescent/quenched galaxies in the COSMOS field. We corroborate the existence of a strong correlation between BHAR (i.e., the X-ray luminosity, L{sub X}) and stellar mass (M{sub *}) for normal star-forming galaxies, though we find a steeper relation than previously reported. We find that starbursts show a factor of three enhancement in BHAR compared to normal SF galaxies (against a factor of six excess in SFR), while quiescents show a deficit of a factor times 5.5 at a given mass. One possible interpretation of this is that the starburst phase does not coincide with cosmologically relevant BH growth, or that starburst-inducing mergers are more efficient at boosting SFR than BHAR. Contrary to studies based on smaller samples, we find that the BHAR/SFR ratio of main-sequence (MS) galaxies is not mass invariant, but scales weakly as M{sub ∗}{sup 0.43±0.09}, implying faster BH growth in more massive galaxies at z∼2. Furthermore, BHAR/SFR during the starburst is a factor of two lower than in MS galaxies, at odds with the predictions of hydrodynamical simulations of merger galaxies that foresee a sudden enhancement of L{sub X}/SFR during the merger. Finally, we estimate that the bulk of the accretion density of the universe at z∼2 is associated with normal star-forming systems, with only ∼6(±1)% and ∼11(±1)% associated with starburst and quiescent galaxies, respectively.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

High-Order Mesh hr-adaptivity for Surface Fitting to Implicit Geometries

We present an ℎ𝑟-adaptivity framework for morphing a given mesh to fit a target surface prescribed as the zero isocontour of a discrete function. In this framework, high-order meshing is posed as a variational minimization problem that depends on the mesh quality prescribed via the target matrix optimization paradigm (TMOP) and position of a subset of mesh nodes with respect to the target surface. The proposed formulation ensures that the variational problem is converged and mesh quality degradation near the surface is limited, even when the mesh topology is incompatible with the target surface. Additionally, a mesh subset-based approach and ℎ-refinement is introduced to efficiently increase fitting accuracy while reducing the computational cost of the mesh morphing problem. The ℎ𝑟-adaptivity technique extends to different element types in two- and three-dimensions, and can be used in existing finite element and spectral element frameworks to obtain high-order body-fitted meshes. Various numerical experiments demonstrate the robustness and accuracy of the fitting approach for problems of practical interest such as Lagrangian hydrodynamics and topology optimization.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Simulations of three-layer Richtmyer–Meshkov mixing in a shock tube

Here, the Richtmyer–Meshkov instability causes perturbations to grow after a shock traverses a fluid density interface. This increases the mixing rate between fluid from either side of the interface. We use the Flash Eulerian hydrodynamic code to investigate alterations when a thin third layer of intermediate density is placed along the interface, effectively creating two adjacent unstable interfaces. This is a common occurrence in engineering applications where a thin barrier initially separates two materials. We find that the width of the mixing layer is similar or slightly reduced; however, the total mass of mixed material can actually increase. The mixing layer becomes more compact and efficient. However, the normalized mixed mass decreases, meaning that finger entrainment becomes more important than in the simple two-layer case. The effect of adding the central layer appears to decrease when the Atwood number is decreased. The Flash results are also benchmarked against two-layer experimental data from a shock tube at the University of Arizona.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2018

Transient Cold Flow Simulation of Fast Fluidized Bed Fuel Reactors for Chemical-Looping Combustion

Circulating fluidized bed in chemical-looping combustion (CLC) is a recent technology that provides great advantage for gas–solid interaction and efficiency. In order to obtain a thorough understanding of this technology and to assess its effectiveness for industrial scale deployment, numerical simulations are conducted. Computational fluid dynamics (CFD) simulations are performed with dense discrete phase model (DDPM) to simulate the gas–solid interactions. CFD commercial software ansysfluent is used for the simulations. Two bed materials of different particle density and diameter, namely the molochite and Fe100, are used in studying the hydrodynamics and particle behavior in a fuel reactor corresponding to the experimental setup of Haider et al. (2016, “A Hydrodynamic Study of a Fast-Bed Dual Circulating Fluidized Bed for Chemical Looping Combustion,” Energy Technol., 4(10), pp. 1254–1262.) at Cranfield University in the UK. Both the simulations show satisfactory agreement with the experimental data for both the static pressure and volume fraction at various heights above the gas inlet in the reactor. It is found that an appropriate drag law should be used in the simulation depending on the particle size and flow conditions in order to obtain accurate results. Furthermore, the simulations demonstrate the ability of CFD/DDPM to accurately capture the physics of circulating fluidized bed-based CLC process at pilot scale which can be extended to industrial scale projects.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1991

A comparison of numerical methods for the Rayleigh equation in unbounded domains

A second-order finite difference and two spectral methods, including a Chebyshev tau and a Chebyshev collocation method were implemented to determine the linear hydrodynamic stability of an unbounded shear flow. The velocity profile of the basic flow in the stability analysis mimicks that of a one-stream free mixing layer. Local and global eigenvalue solution methods are used to determine individual eigenvalues and the eigenvalue spectrum, respectively. The calculated eigenvalue spectrum includes a discrete mode, a continuous spectrum associated with the equation singularity and a continuous spectrum associated with the domain unboundedness. The efficiency and the accuracy of these discretization methods in the prediction of the eigensolutions of the discrete mode were evaluated by comparison with a conventional shooting procedure. Their capabilities in mapping out the continuous eigenvalue spectra are also discussed.

Liou, W. W.↗

DOE OSTI2014

Aeroelastic Modeling of Offshore Turbines and Support Structures in Hurricane-Prone Regions (Poster)

US offshore wind turbines (OWTs) will likely have to contend with hurricanes and the associated loading conditions. Current industry standards do not account for these design load cases (DLCs), thus a new approach is required to guarantee that the OWTs achieve an appropriate level of reliability. In this study, a sequentially coupled aero-hydro-servo-elastic modeling technique was used to address two design approaches: 1.) The ABS (American Bureau of Shipping) approach; and 2.) The Hazard Curve or API (American Petroleum Institute) approach. The former employs IEC partial load factors (PSFs) and 100-yr return-period (RP) metocean events. The latter allows setting PSFs and RP to a prescribed level of system reliability. The 500-yr RP robustness check (appearing in [2] and [3] upcoming editions) is a good indicator of the target reliability for L2 structures. CAE tools such as NREL's FAST and Bentley's' SACS (offshore analysis and design software) can be efficiently coupled to simulate system loads under hurricane DLCs. For this task, we augmented the latest FAST version (v. 8) to include tower aerodynamic drag that cannot be ignored in hurricane DLCs. In this project, a 6 MW turbine was simulated on a typical 4-legged jacket for a mid-Atlantic site. FAST-calculated tower base loads were fed to SACS at the interface level (transition piece); SACS added hydrodynamic and wind loads on the exposed substructure, and calculated mudline overturning moments, and member and joint utilization. Results show that CAE tools can be effectively used to compare design approaches for the design of OWTs in hurricane regions and to achieve a well-balanced design, where reliability levels and costs are optimized.

17 WIND ENERGY↗