Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Energy minimization”

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 955 records · Page 53

An approximate solution for interlaminar stresses in composite laminates

An efficient approximate solution for interlaminar stresses in finite width, symmetric and unsymmetric laminated composites subjected to axial and/or bending loads is presented. The solution is based upon statically admissible stress fields which take into consideration local property mismatch effects and global equilibrium requirements. Unknown constants in the assumed stress states are determined through minimization of the laminate complementary energy. Typical results are presented for through-thickness and interlaminar stress distributions for angle-ply and cross-ply laminates subjected to axial loading. It is shown that the present formulation represents an improved, efficient approximate solution for interlaminar stresses.

Rose, Cheryl A.↗

Short-Term Aging of NeFeB Magnets for Stirling Linear Alternator Applications

NeFeB type magnets have been proposed for use in free piston Stirling engine driven, linear alternators to generate electric power during long duration space missions. These type of materials provide the highest energy product commercial magnets, thus minimizing alternator size or mass, but do not provide the high temperature stability of magnetic properties found in the SmCo type magnets. Therefore, to apply the NeFeB type magnets at elevated temperatures to multiyear space missions, their long-term aging characteristics must be determined. This report presents 200 hr aging data for six types of NeFeB magnets selected from three manufacturers. Aging was performed under vacuum at 150 C, with a steady demagnetizing field of 5 kOe applied. From the data produced by this short-term aging run, candidate magnet types were selected for a planned 12,000 hr long-term run. Depending on the manufacturer's magnet type, remanence losses observed ranged from 0 to 7%, when measured at 120 C on an established recoil line. Also, intrinsic coercivity losses up to about 4% were observed for the M-H curve at 120 C. In some cases, these coercivity losses were not recoverable by recharge of the magnet, indicating a structural change of the material.

Niedra, Janis M.↗

Structures Self-Assembled Through Directional Solidification

Nanotechnology has created a demand for new fabrication methods with an emphasis on simple, low-cost techniques. Directional solidification of eutectics (DSE) is an unconventional approach in comparison to low-temperature biomimetic approaches. A technical challenge for DSE is producing microstructural architectures on the nanometer scale. In both processes, the driving force is the minimization of Gibb's free energy. Selfassembly by biomimetic approaches depends on weak interaction forces between organic molecules to define the architectural structure. The architectural structure for solidification depends on strong chemical bonding between atoms. Constituents partition into atomic-level arrangements at the liquid-solid interface to form polyphase structures, and this atomic-level arrangement at the liquid-solid interface is controlled by atomic diffusion and total undercooling due to composition (diffusion), kinetics, and curvature of the boundary phases. Judicious selection of the materials system and control of the total undercooling are the keys to producing structures on the nanometer scale. The silicon-titanium silicide (Si-TiSi2) eutectic forms a rod structure under isothermal cooling conditions. At the NASA Glenn Research Center, directional solidification was employed along with a thermal gradient to promote uniform rods oriented with the thermal gradient. The preceding photomicrograph shows the typical transverse microstructure of a solidified Si-TiSi2 eutectic composition. The dark and light gray regions are Si and TiSi2, respectively. Preferred rod orientation along the thermal gradient was poor. The ordered TiSi2 rods have a narrow distribution in diameter of 2 to 3 m, as shown. The rod diameter showed a weak dependence on process conditions. Anisotropic etch behavior between different phases provides the opportunity to fabricate structures with high aspect ratios. The photomicrographs show the resulting microstructure after a wet chemical etch and a dry plasma etch. The wet chemical etches the silicon away, exposing the TiSi2 rods, whereas plasma etching preferentially etches the Si-TiSi2 interface to form a crater. The porous architectures are applicable to fabricating microdevices or creating templates for part fabrication. The porous rod structure can serve as a platform for fabricating microplasma devices for propulsion or microheat exchangers and for fabricating microfilters for miniatured chemical reactors. Although more work is required, self-assembly from DSE can have a role in microdevice fabrication.

Dynys, Frederick W.↗

Small Body GN&C Research Report: A G-Guidance Enhancement to Increase Mission Functionality

G-Guidance is a robust G&C (guidance and control) algorithm developed under the small-body GN&C task. The G-Guidance scheme utilizes a model predictive control approach, along with a convexification of the governing dynamics, control constraints, and trajectory/state constraints. This report details an enhancement to the FY2005 G-Guidance algorithm; the addition of a fire-second approach helps to nullify velocity errors and hit desired final velocities much more precisely than the original Fire-first scheme developed in FY2005. The enhancement preserves the primary benefit of the algorithm, which is to ensure required thruster silent times during trajectory maneuvers. The fire-second scheme increases the versatility of the algorithm for missions employing G-Guidance. For instance, a landing sequence could employ the fire-second scheme to ensure a null final velocity at landing, whereas an ascent sequence could utilize a fire-first scheme to ensure there is no drift of the spacecraft toward the ground at the beginning of the maneuver. Examples are provided within to demonstrate a fire-first versus fire-second guidance scheme. As in the existing G-Guidance algorithm, the examples and schemes incorporate gravity models and thruster firing times into discrete dynamics that are solved as a optimal control problem to minimize fuel consumption or thruster energy expenditure.

velocity errors↗

Small Body GN&C Research Report: A Guidance and Control Technique for Small-Body Proximity Operations with Guaranteed Guidance Resolvability and Required Thruster Silent Time

The guidance and control (G&C) algorithms for enabling small-body proximity operations are developed by using a model predictive control approach along with a convexification of the governing dynamics, control constraints, and trajectory/state constraints. The open-loop guidance is solved ahead of time or in a resolvable, real-time manner through the use of PWG (Pseudo Way-point Generation), a technique developed in this research. The PWG scheme ensures required thruster silent times during trajectory maneuvers. The feedback control is implemented to track the PWG trajectories in a manner that guarantees the resolvability for the open-loop problem, enabling the ability to update the G&C in a model-predictive manner. The schemes incorporate gravity models and thruster ring times into discrete dynamics that are solved as a optimal control problem to minimize fuel consumption or thruster energy expenditure. The optimal control problem is cast as an LMI (Linear Matrix Inequality) and then solved through Semi-Definite Programming techniques in a computationally efficient manner that provides convergence and constraint guarantees.

guidance and control (G&C)↗

CFD Study of the J2-X LOX Two-Bladed Inducer in the Presence of a Circumferential Groove

Results are presented of a computational fluid dynamics (CFD) study done in support of water flow experiments of the J-2X Oxidizer Turbopump (OTP) 2-bladed inducer with a circumferential groove that were conducted at Marshall Space Flight Center (MSFC). Sub-scale water flow testing results indicate that the circumferential groove greatly reduces synchronous cavitation and subsequent bearing loads at a minimal performance cost, but the energy reappears as high order cavitation (HOC) that spans a relatively large frequency range. Thus, HOC may have implications for the full-scale OTP inducer in terms of reduced structural margin at higher mode frequencies. Simulations using the LOCI-Stream CFD program were conducted in order to explore the root physical cause of the HOC. It was found that the axial recirculation pattern in the circumferential groove causes high-swirl fluid to interact with the nearly-axial incoming fluid just above the inducer blades. The high-shear interface between the fluids is Kelvin-Helmholtz unstable, resulting in trains of low pressure regions or 'pearls' forming near the upstream edge of the groove. When the pressure in these regions becomes low enough and they get cut by the blade leading edge, HOC is thought to occur. Although further work is required, the numerical models indicate that the root cause of HOC is hydrodynamic. That is, the pearls are always present, even when cavitation is not. Comparisons to ongoing water flow experiments will be discussed, as well as predictions for the full-scale OTP inducer.

Coker, Robert↗

Sustainably Sourced, Thermally Resistant, Radiation Hard Biopolymer

This material represents a breakthrough in the production, manufacturing, and application of thermal protection system (TPS) materials and radiation shielding, as this represents the first effort to develop a non-metallic, non-ceramic, biomaterial-based, sustainable TPS with the capability to also act as radiation shielding. Until now, the standing philosophy for radiation shielding involved carrying the shielding at liftoff or utilizing onboard water sources. This shielding material could be grown onboard and applied as needed prior to different radiation landscapes (commonly seen during missions involving gravitational assists). The material is a bioplastic material. Bioplastics are any combination of a biopolymer and a plasticizer. In this case, the biopolymer is a starch-based material and a commonly accessible plasticizer. Starch molecules are composed of two major polymers: amylase and amylopectin. The biopolymer phenolic compounds are common to the ablative thermal protection system family of materials. With similar constituents come similar chemical ablation processes, with the potential to have comparable, if not better, ablation characteristics. It can also be used as a flame-resistant barrier for commercial applications in buildings, homes, cars, and heater firewall material. The biopolymer is observed to undergo chemical transformations (oxidative and structural degradation) at radiation doses that are 1,000 times the maximum dose of an unmanned mission (10-25 Mrad), indicating that it would be a viable candidate for robust radiation shielding. As a comparison, the total integrated radiation dose for a three-year manned mission to Mars is 0.1 krad, far below the radiation limit at which starch molecules degrade. For electron radiation, the biopolymer starches show minimal deterioration when exposed to energies greater than 180 keV. This flame-resistant, thermal-insulating material is non-hazardous and may be sustainably sourced. It poses no hazardous waste threats during its lifecycle. The material composition is radiation-tolerant up to megarad doses, indicating its use as a radiation shielding material. It is lightweight, non-metallic, and able to be mechanically densified, permitting a tunable gradient of thermal and radiation protection as needed. The dual-use (thermal and radiation shielding), sustainable nature of this material makes it suitable for both industrial applications as a sustainable/green building material, and for space applications as thermal protection material and radiation shield.

Pugel, Diane↗

Statistical Characterization of 18650-Format Lithium-Ion Cell Thermal Runaway Energy Distributions

Effective thermal management systems, designed to handle the impacts of thermal runaway (TR) and to prevent cell-to-cell propagation, are key to safe operation of lithium-ion (Li-ion) battery assemblies. Critical factors for optimizing these systems include the total energy released during a single cell TR event and the fraction of the total energy that is released through the cell casing vs. through the ejecta material. A unique calorimeter was utilized to examine the TR behavior of a statistically significant number of 18650-format Li-ion cells with varying manufacturers, chemistries, and capacities. The calorimeter was designed to contain the TR energy in a format conducive to discerning the fractions of energy released through the cell casing vs. through the ejecta material. Other benefits of this calorimeter included the ability to rapidly test of large quantities of cells and the intentional minimization of secondary combustion effects. High energy (270 Wh/kg) and moderate energy (200 Wh/kg) 18650 cells were tested. Some of the cells had an imbedded short circuit (ISC) device installed to aid in the examination of TR mechanisms under more realistic conditions. Other variations included cells with bottom vent (BV) features and cells with thin casings (0.22 1/4m). After combining the data gathered with the calorimeter, a statistical approach was used to examine the probability of certain TR behavior, and the associated energy distributions, as a function of capacity, venting features, cell casing thickness and temperature.

Thermal Runaway Energy Distribution↗

Pale Orange Dots: the Impact of Organic Haze on the Habitability and Detectability of Earthlike Exoplanets

Hazes are common in known planetary atmospheres, and geochemical evidence suggests that early Earth occasionally supported an organic haze with significant environmental and spectral consequences. The UV spectrum of the parent star drives organic haze formation through methane photochemistry. We use a 1D photochemical-climate model to examine production of fractal organic haze on Archean Earth-analogs in the habitable zones of several stellar types: the modern and early Sun, AD Leo (M3.5V), GJ 876 (M4V), epsilon Eridani (K2V), and sigma Bootis (F2V). For Archean-like atmospheres, planets orbiting stars with the highest UV (ultraviolet) fluxes do not form haze because of the formation of photochemical oxygen radicals that destroy haze precursors. Organic hazes impact planetary habitability via UV shielding and surface cooling, but this cooling is minimized around M dwarfs, whose energy is emitted at wavelengths where organic hazes are relatively transparent. We generate spectra to test the detectability of haze. For 10 transits of a planet orbiting GJ 876 observed by the James Webb Space Telescope, haze makes gaseous absorption features at wavelengths less than 2.5 microns times 2-10 sigma shallower than a haze-free planet, and methane and carbon dioxide are detectable at greater than 5 sigma. A haze absorption feature can be detected at 5 sigma near 6.3 microns, but a higher signal-to-noise ratio is needed to distinguish haze from adjacent absorbers. For direct imaging of a planet at 10 parsecs using a coronagraphic 10 meter class ultraviolet-visible-near-infrared telescope, a UV-blue haze absorption feature would be strongly detectable at greater than 12 sigma in 200 hours.

Arney, Giada N.↗

Fitting strategies of accretion column models and application to the broadband spectrum of Cen X-3

Due to the complexity of modeling the radiative transfer inside the accretion columns of neutron star binaries, their X-ray spectra are still commonly described with phenomenological models, e.g., a cut off power law. While the behavior of these models is well understood and they allow for a comparison of different sources and studying source behavior, the extent to which the underlying physics can be derived from the model parameters is very limited. During recent years, several physically motivated spectral models have been developed to overcome these limitations. Their application, however, is generally computationally much more expensive and they require a high number of parameters which are difficult to constrain. In particular, Becker & Wolff(2007) presented an analytical solution to the radiative transfer equation inside the accretion column assuming a velocity profile that is linear in the optical depth. An implementation of this model that is both fast and accurate enough to be fitted to observed spectra is available in XSPEC. The main difficulty of this implementation is that some solutions violate energy conservation and therefore have to be rejected by the user. We propose a novel fitting strategy that ensures energy conservation during theχ2-minimization which simplifies the application of the model considerably. We demonstrate this approach as well as a study of possible parameter degeneracies with a comprehensive Markov-Chain Monte Carlo analysis of the complete parameter space for a combined NuSTAR and Swift/XRT dataset of Cen X-3.The derived accretion-flow structure features a small column radius of∼63 m and a spectrum dominated by bulk-Comptonization of bremsstrahlung seed photons, in agreement with previous studies.

X-rays: binaries↗

Quantifying the Effect of Cosmic Ray Showers on the X-IFU Energy Resolution

The X-ray Integral Field Unit (X-IFU) will operate an array of more than 3000 Transition Edge Sensor pixels at 90 mK with an unprecedented energy resolution of 2.5 eV at 7 keV. In space, primary cosmic rays and secondary particles produced in the instrument structure will continuously deposit energy on the detector wafer and induce fluctuations on the pixels’ thermal bath. We have investigated through simulations of the X-IFU readout chain how these fluctuations eventually influence the energy measurement of X-ray photons. Realistic timelines of thermal bath fluctuations at different positions in the array are generated as a function of a thermal model and the expected distribution of the deposited energy of the charged particles. These are then used to model the TES response to these thermal perturbations and their influence on the onboard energy reconstruction process. Overall, we show that with adequate heatsinking, the main energy resolution degradation effect remains minimal and within the associated resolution allocation of 0.2 eV. We further study how a dedicated triggering algorithm could be put in place to flag the rarer large thermal events.

P. Peille↗

Fitting Strategies of Accretion Column Models and Application to the Broadband Spectrum of Cen X-3

Due to the complexity of modeling the radiative transfer inside the accretion columns of neutron star binaries, their X-ray spectra are still commonly described with phenomenological models, e.g., a cut off power law. While the behavior of these models is well understood and they allow for a comparison of different sources and studying source behavior, the extent to which the underlying physics can be derived from the model parameters is very limited. During recent years, several physically motivated spectral models have been developed to overcome these limitations. Their application, however, is generally computationally much more expensive and they require a high number of parameters which are difficult to constrain. In particular, Becker & Wolff(2007) presented an analytical solution to the radiative transfer equation inside the accretion column assuming a velocity profile that is linear in the optical depth. An implementation of this model that is both fast and accurate enough to be fitted to observed spectra is available in XSPEC. The main difficulty of this implementation is that some solutions violate energy conservation and therefore have to be rejected by the user. We propose a novel fitting strategy that ensures energy conservation during theχ2-minimization which simplifies the application of the model considerably. We demonstrate this approach as well as a study of possible parameter degeneracies with a comprehensive Markov-Chain Monte Carlo analysis of the complete parameter space for a combined NuSTAR and Swift/XRT data-set of Cen X-3.The derived accretion-flow structure features a small column radius of∼63 m and a spectrum dominated by bulk-Comptonization of bremsstrahlung seed photons, in agreement with previous studies

Philipp Thalhammer↗

Progress on Optimizing Wind Farms and Rotor Designs Using Adjoints

Modern wind plants are increasingly tasked with multiple performance objectives. In addition to designing plants that maximize power output and minimize the levelized cost of energy (LCOE), the design and operation of wind plants is increasingly influenced by challenges regarding grid integration of variable generation renewables. This places a growing emphasis on making wind plants more controllable and predictable. WindSE is a Reynolds-averaged Navier-Stokes (RANS) model designed around analytical gradient and adjoint methods, with the ability to capture terrain-induced effects, as shown in Figure 1. The recent addition of an unsteady solver with an actuator line method (ALM) and ongoing work to enable massively parallel optimizations gives it a unique niche to explore coupled plant-level controls and design problems. This code is an open source python package built on the FEniCS framework that utilizes fast, parallel PETSc solvers to model fluid flow throughout wind-farm scale domains. Two recent studies performed using WindSE demonstrate the capability to optimize under a wide variety of flow conditions and objective functions. In the first, we present an optimization focused on modifying the layout of a wind farm with a fixed number of turbines for maximum total power output [1]. This study highlights the ability to quickly perform simulations using the steady Navier-Stokes solver combined with rotors represented as actuator disks while also stressing the importance of capturing terrain-induced effects. Gradient-based optimization using the RANS equations is viable due to the inclusion of efficiently computed adjoint derivatives. We interpret the physical results of the optimal layout and also discuss the computational cost of scaling to larger problems. In the second study, we present the capabilities of the unsteady Navier-Stokes solver, where rotor-blade profiles represented by actuator lines are optimized to enhance wake steering effects and overall power production [2]. We quantify the wind plant performance gains obtained from this type of simultaneous control co-design optimization as compared to optimizing the blade design and yaw independently. Figure 2 shows the differences between a baseline two-turbine system and an optimized system where we fine-tune the blade chord profile. Results and challenges from each study are quickly summarized and used to motivate the current development efforts within WindSE. Current and future work is focused on enabling higher-resolution studies with more degrees of freedom through parallelization of both the simulation and optimization algorithms. We present benchmarking results to show that WindSE performs well in both weak- and strong-scaling tests and further demonstrate that the optimizer obtains the same convergence rates in both shared- and distributed-memory environments. Using larger wind farms, we can study deep-array effects within an optimization context, allowing the use of objective functions that have been previously unstudied. As an example, we present ongoing work on a blockage metric which characterizes the loss of available kinetic energy due to wake effects from multiple upstream turbines.

adjoint optimization↗

Demand-adaptive Transit Design for Urban Transportation Hubs

In this study, we proposed a novel three-stage framework for planning the optimal demand-adaptive transit (DAT) at urban transportation hubs. Given the potential trip demand and road traffic condition, the proposed framework sequentially generates the optimal set of candidate routes, combines the outgoing routes and incoming routes at the hub, and derives the optimal fleet size and corresponding route frequency under the fixed budget. In particular, we build the route generation algorithm which maximizes passenger demand coverage with travel time deviation constraint. And a heuristic algorithm is further developed which yields near-optimal operation routes for real-time demand. The fleet optimization problem is formulated to minimize the weighted cost of energy savings, operation cost and trip revenue. We conduct comprehensive numerical experiments for planning DAT with electric buses at JFK airport in NYC using NYC taxi and for-hire vehicle trip data and GoogleMap speed data. The results show the superior performance of the proposed route generation algorithm which is able to cover citywide passenger demand with only 61 DAT routes. The results also suggest that the proposed DAT planning framework may serve over 47% of existing taxi and FHV demand by operating 18 routes using the fleet of 62 electric buses.

demand adaptive↗

Performance assessment of 3D printed multi-material energy absorber for automotive bumper: pedestrian lower extremity protection

Designing an energy absorber for automotive bumpers involves balancing low-speed and high-speed impacts to ensure safety, reduce repair costs, and meet regulatory standards. Here, this study explores a novel design using multi-material 3D printing and structural optimization to fabricate a lightweight and cost-efficient energy absorber. The design effectively dissipates energy in low-speed collisions and minimizes force transmission in high-speed pedestrain impacts, helping to meet both safety and performance requirements. The energy absorber design combines 20% carbon fiber-reinforced acrylonitrile butadiene styrene (CF-ABS) and thermoplastic polyurethane (TPU) for optimal stiffness and flexibility. It uses 3D-printed lattice structures optimized through finite element simulations to help meet both low-speed and high-speed impact requirements. Full-scale energy absorbers were 3D-printed using optimized CF-ABS/TPU blends and tested under high-speed impact using the Flexible Pedestrian Legform Impactor (Flex-PLI). For fair comparison, a baseline bumper with a traditional triangular lattice structure, also 3D-printed from the same CF-ABS/TPU materials, was similarly tested. Interestingly, both the optimized and baseline 3D-printed energy absorbers showed nearly identical performance, successfully meeting injury limits. Their performances were also benchmarked against an injection-molded energy absorber. While both 3D-printed and injection-molded designs met injury limits, the 3D-printed absorber exhibited a higher tibia bending moment, indicating an opportunity for further optimization. A Techno-Economic Analysis compared the costs of producing energy absorbers using traditional manufacturing and 3D printing. The analysis highlighted that 3D printing offers cost benefits for low to medium production volumes, with the total cost per energy absorber at ∼ $\$$74, compared to traditional methods that become economical beyond 2000 units.

Additive manufacturing↗

Benchmarking intra-nuclear cascade models for neutrino scattering with relativistic optical potentials

The description of final-state interactions (FSI) in the large phase space probed in neutrino experiments poses a great challenge. In neutrino experiments, which operate under semi-inclusive conditions, cascade models are commonly used for this task, while under exclusive conditions FSI can be treated with relativistic optical potentials (ROP). We formulate conditions under which the ROP approach and cascade model can be directly compared. We feed the NEUT cascade with events from a relativistic distorted-wave impulse approximation calculation that uses the real part of an optical potential. Cuts on the missing energy of the resulting events are applied to define a set of events that can be directly compared to RDWIA calculations with the full optical potential. The NEUT cascade and ROP agree for proton kinetic energies $T_p > 150$ MeV for carbon, oxygen and calcium nuclei when a realistic nuclear density is used to introduce events in the cascade. For $T_p < 100$ MeV the ROP and NEUT cross sections differ in shape and differences in magnitude are larger than 50 %. Single transverse variables allow to distinguish different approaches to FSI, but due to a large non-QE contribution the comparison to T2K data does not give an unambiguous view of FSI. We discuss electron scattering and argue that with a cut in missing energy FSI can be studied with minimal confounding factors in e.g. $e4\nu$. The agreement of the ROP and NEUT for T2K conditions lends confidence to these models as a tool in oscillation analyses for sufficiently large nucleon kinetic energies. These results urge for caution when a cascade model is applied for small nucleon energies. The assessment of model assumptions relevant to this region are strongly encouraged. This paper provides novel constraints on cascade models from proton-nucleus scattering that can be easily applied to other neutrino event generators.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A Consistent BGK Model with Velocity-Dependent Collision Frequency for Gas Mixtures

Abstract We derive a multi-species BGK model with velocity-dependent collision frequency for a non-reactive, multi-component gas mixture. The model is derived by minimizing a weighted entropy under the constraint that the number of particles of each species, total momentum, and total energy are conserved. We prove that this minimization problem admits a unique solution for very general collision frequencies. Moreover, we prove that the model satisfies an H-Theorem and characterize the form of equilibrium.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Interleaved Cuk Converter Wave Energy System With Advanced Control and Grid Support Functions for Seamless Integration

This paper presents an innovative wave energy conversion system that integrates an interleaved Cuk converter with advanced nonlinear control for seamless grid integration. The system efficiently extracts power on the DC side using the interleaved Cuk converter, while an inverter manages power transfer to the grid or load on the AC side. A nonlinear control architecture, based on the Lyapunov energy function, ensures stable and optimal operation under varying wave conditions, effectively addressing the challenges of variability and unpredictability inherent in wave energy. The proposed system also incorporates features to enhance power quality and minimize losses, making it a robust solution for renewable energy integration. The system's effectiveness in harnessing wave energy and achieving smooth grid integration is validated through comprehensive computer simulations in MATLAB/Simulink, with case study results demonstrating its capabilities.

HYDRO ENERGY,TIDAL AND WAVE POWER↗