SEARCH · Engineering Papers

Results for “Gradient Limiters”

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 1,117 records · Page 62

All sources NASA NTRS DOE OSTI

DOE OSTI2016

Gradient twinned 304 stainless steels for high strength and high ductility

Gradient materials often have attractive mechanical properties that outperform uniform microstructure counterparts. It remains a difficult task to investigate and compare the performance of various gradient microstructures due to the difficulty of fabrication, the wide range of length scales involved, and their respective volume percentage variations. We have investigated four types of gradient microstructures in 304 stainless steels that utilize submicrotwins, nanotwins, nanocrystalline-, ultrafine- and coarse-grains as building blocks. Tensile tests reveal that the gradient microstructure consisting of submicrotwins and nanotwins has a persistent and stable work hardening rate and yields an impressive combination of high strength and high ductility, leading to a toughness that is nearly 50% higher than that of the coarse-grained counterpart. Ex- and in-situ transmission electron microscopy indicates that nanoscale and submicroscale twins help to suppress and limit martensitic phase transformation via the confinement of martensite within the twin lamellar. Twinning and detwinning remain active during tensile deformation and contribute to the work hardening behavior. We discuss the advantageous properties of using submicrotwins as the main load carrier and nanotwins as the strengthening layers over those coarse and nanocrystalline grains. Furthermore, our work uncovers a new gradient design strategy to help metals and alloys achieve high strength and high ductility.

36 MATERIALS SCIENCE↗

DOE OSTI2017

Mixing dynamics in bubbling fluidized beds

Solids mixing affects thermal and concentration gradients in fluidized bed reactors and is, therefore, critical to their performance. Despite substantial effort over the past decades, understanding of solids mixing continues to be lacking because of technical limitations of diagnostics in large pilot and commercial‐scale reactors. This study is focused on investigating mixing dynamics and their dependence on operating conditions using computational fluid dynamics simulations. Toward this end, fine‐grid 3D simulations are conducted for the bubbling fluidization of three distinct Geldart B particles (1.15 mm LLDPE, 0.50 mm glass, and 0.29 mm alumina) at superficial gas velocities U/Umf = 2–4 in a pilot‐scale 50 cm diameter bed. The Two‐Fluid Model (TFM) is employed to describe the solids motion efficiently while bubbles are detected and tracked using MS3DATA. Detailed statistics of the flow‐field in and around bubbles are computed and used to describe bubble‐induced solids micromixing: solids upflow driven in the nose and wake regions while downflow along the bubble walls. Further, within these regions, the hydrodynamics are dependent only on particle and bubble characteristics, and relatively independent of the global operating conditions. Based on this finding, a predictive mechanistic, analytical model is developed which integrates bubble‐induced micromixing contributions over their size and spatial distributions to describe the gross solids circulation within the fluidized bed. Finally, it is shown that solids mixing is affected adversely in the presence of gas bypass, or throughflow, particularly in the fluidization of heavier particles. This is because of inefficient gas solids contacting as 30–50% of the superficial gas flow escapes with 2–3× shorter residence time through the bed. This is one of the first large‐scale studies where both the gas (bubble) and solids motion, and their interaction, are investigated in detail and the developed framework is useful for predicting solids mixing in large‐scale reactors as well as for analyzing mixing dynamics in complex reactive particulate systems. © 2017 American Institute of Chemical Engineers AIChE J , 63: 4316–4328, 2017

Bakshi, A.↗

DOE OSTI2025

A minor respiratory process with major global implications: is atmospheric methane oxidation in tree stems driven by stem respiration rather than microbial methanotrophy?

Tree stem surfaces are widely recognized as sites of carbon dioxide (CO₂) efflux and oxygen (O₂) influx, reflecting the dynamics of aerobic respiration of photosynthate substrates, such as sugars, delivered via the phloem. Stems are also largely considered passive conduits for methane (CH₄) produced in anoxic soils via microbial methanogenesis, where CH₄ is thought to be transported upward through the transpiration stream and/or diffusion and emitted through stem surfaces and the canopy. However, recent observations from dynamic stem chambers suggest that stems may also act as active sinks for atmospheric CH₄. Despite these findings, the extent and drivers of stem CH₄ consumption remain poorly characterized across biomes, species, and environmental gradients, and its quantitative relationship to stem respiration has not been established. Moreover, previous studies captured only snapshot fluxes, leaving diurnal patterns of CH₄ exchange uncharacterized. Here, we address these limitations by combining real-time measurements of stem CH₄ and O₂ uptake under ambient conditions in a California cherry tree, using a dynamic stem gas exchange system with three chambers receiving a continuous flow of ambient air and automated chamber and reference air sampling every 10 min. Our results confirm that stems of upland trees can actively consume both atmospheric CH₄ and O₂, but with decreasing temperature sensitivity as daily temperatures increase. Early mornings were marked by rapid influxes of both gases, followed by declining uptake as temperatures rose further. Methane uptake was tightly coupled with O₂ influx and represented a minor (0.012% ± 0.002%) fraction of stem respiratory activity, as determined by concurrent O₂ uptake. These findings suggest that while atmospheric CH₄ oxidation is a minor respiratory process in stems, it is strongly linked with stem physiological activity. This challenges the current assumption that terrestrial CH₄ uptake is driven solely by microbial methanotrophy and raises the possibility that living stem tissues may contribute to CH₄ oxidation through an as-yet-unidentified plant-based mechanism.

Atmospheric greenhouse gases↗

DOE OSTI2017

Spin caloritronic nano-oscillator

Energy loss due to ohmic heating is a major bottleneck limiting down-scaling and speed of nano-electronic devices, and harvesting ohmic heat for signal processing is a major challenge in modern electronics. Here, we demonstrate that thermal gradients arising from ohmic heating can be utilized for excitation of coherent auto-oscillations of magnetization and for generation of tunable microwave signals. The heat-driven dynamics is observed in Y 3 Fe 5 O 12 /Pt bilayer nanowires where ohmic heating of the Pt layer results in injection of pure spin current into the Y 3 Fe 5 O 12 layer. This leads to excitation of auto-oscillations of the Y 3 Fe 5 O 12 magnetization and generation of coherent microwave radiation. Thus, our work paves the way towards spin caloritronic devices for microwave and magnonic applications.

36 MATERIALS SCIENCE↗

DOE OSTI2013

Five-field simulations of peeling-ballooning modes using BOUT++ code

The simulations of edge localized modes (ELMs) with a 5-field peeling-ballooning (P-B) model using BOUT++ code are reported in this paper. In order to study the particle and energy transport in the pedestal region, the pressure equation is separated into ion density and ion and electron temperature equations. Through the simulations, the length scale L{sub n} of the gradient of equilibrium density n{sub i0} is found to destabilize the P-B modes in ideal MHD model. With ion diamagnetic effects, the growth rate is inversely proportional to n{sub i0} at medium toroidal mode number n. For the nonlinear simulations, the gradient of n{sub i0} in the pedestal region can more than double the ELM size. This increasing effect can be suppressed by thermal diffusivities χ{sub ∥}, employing the flux limited expression. Thermal diffusivities are sufficient to suppress the perturbations at the top of pedestal region. These suppressing effects lead to smaller ELM size of P-B modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2013

Observation of multipactor suppression in a dielectric-loaded accelerating structure using an applied axial magnetic field

Efforts by a number of institutions to develop a Dielectric-Loaded Accelerating (DLA) structure capable of supporting high gradient acceleration when driven by an external radio frequency source have been ongoing over the past decade. Single surface resonant multipactor has been previously identified as one of the major limitations on the practical application of DLA structures in electron accelerators. In this paper, we report the results of an experiment that demonstrated suppression of multipactor growth in an X-band DLA structure through the use of an applied axial magnetic field. This represents an advance toward the practical use of DLA structures in many accelerator applications.

43 PARTICLE ACCELERATORS↗

DOE OSTI2021

High-performance ECRH at W7-X: experience and perspectives

The second operation phase of W7-X (OP1.2) showed the potential of exclusively electron cyclotron resonance heating (ECRH)-sustained plasma operations in stellarators. Employing multi-pass ECRH scenario in the second harmonic O-mode (O2-ECRH), stationary densities of up to 1.4 × 10 20 m -3 could be achieved. This scenario also made stationary divertor detachment possible, which is a reactor-relevant scenario for power and particle exhaust. At high densities and with sufficiently high density gradients for an improved ion confinement, the coupling between the electrons and ions was strong enough to bring the ion temperature to values above 3 keV and to the neoclassical limit for some magnetic configurations, thus enabling to test the W7-X neoclassical optimization. The planned enhancement of the ECRH performance will enable to advance towards reactor-relevant beta values and to investigate their stability and confinement of fast particles, which is a priority goal of W7-X.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Tokamak level performance in the optimized stellarator Wendelstein 7-X with stable peaked density profiles

Progress towards achieving net gain from nuclear fusion is typically characterized by the triple product 𝑛⁢𝑇⁢𝜏. Although stellarators have several benefits for a reactor, such as the inherent steady-state capability and the lack of disruptions, the maximum achieved triple product has historically lagged behind that of tokamaks. In the large optimized stellarator Wendelstein 7-X, a triple product of (1.10 ± 0.15) × 10 20 m −3 ⁢ keVs has recently been achieved and held stable for 1.9 s. This marks an important step in stellarator core plasma performance and puts it in line with 𝐻-mode performance in similar sized tokamaks. The progress in Wendelstein 7-X performance is due to the suppression of ion-temperature-gradient turbulence by strong density gradients. In the record case, these are maintained by a fine balance of the core particle transport and central fueling using neutral beam injection. The duration of the scenario was limited in time only by technical restrictions of the neutral beam injectors and control of the divertor heat loads. This paper presents the result in the context of the leading magnetic confinement devices along with details of the scenario.

Bannmann, Sebastian [Max Planck Institute for Plas↗

DOE OSTI2014

SU-E-J-251: Fast MR-Based DRR Generation Using Highly Undersampled 3D Radial Trajectories

Purpose: The construction of a digitally reconstructed radiograph (DRR) from a magnetic resonance image (MRI) is possible if the cortical bone signal can be acquired and separated from air and soft tissue. This may be accomplished by subtracting a long echo-time, in-phase, gradient echo (GRE) image volume from an ultra-short echo time free induction decay (FID) image to produce a bone-enhanced (BE) image that reveals cortical bone. One limitation of this approach is the length of time required for data acquisition, which can limit the quality of the DRRs due to patient and organ motion. This study aimed to significantly reduce the acquisition time without compromising DRR quality. Methods: Brain data were acquired from two volunteers using a 3T MR scanner (Ingenia, Philips Healthcare). The FID and GRE images were acquired in a single acquisition using a 3D radial readout sequence with the following parameters: TE1=0.142ms (ultra-short), TE2=2.197ms (nearly in-phase), 2*2*2mm3 isotropic voxels, 250*250*250mm3 FOV. To reduce the acquisition time, k-space was sampled at 75, 50 and 25% of a full 3D sphere . The TE2 image was subtracted from the TE1 image to generate the BE images. The BE images were used to generate DRRs using the Pinnacle treatment planning system (Philips-version 9.2). The quality of the DRRs was evaluated qualitatively by 5 board certified medical physicists for clinical usefulness. Results: The acquisition time for 75, 50 and 25% sampling schemes were 219s, 146s, and 73s, respectively, the latter of which was a four-fold reduction in scan time compared to a 300s fully-sampled acquisition. All DRRs obtained were of acceptable quality and were shown to have sufficient information for clinical 2D image matching. Conclusion: Undersampling k-space while maintaining the same range of frequency information results in significantly reduced scan time and clinically acceptable DRR image quality. Drs. B Traughber and R Muzic have research support from Philips Healthcare. Drs. M Traughber and L Hu are employees of Philips Healthcare.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2025

Irradiation Vehicles for Evaluating SiC/SiC Cladding Bowing Under Neutron Flux Gradients

Silicon carbide fiber–reinforced silicon carbide matrix (SiC/SiC) composites are among the most promising candidates for long term accident-tolerant nuclear fuel cladding. A key challenge related to their deployment is lateral bowing caused by differential radiation-induced swelling under dose or temperature gradients, which could obstruct coolant flow or interfere with control rod/blade movements. Although previous modeling efforts have predicted bowing behavior in light-water reactor (LWR) environments, experimental validation remains limited, especially at prototypic LWR temperatures. This study addresses that gap by irradiating six reduced-length SiC/SiC cladding tubes (~600 mm) in the High Flux Isotope Reactor (HFIR) at ~300°C, which is representative of LWR conditions. The tubes will be housed in a sealed vessel with an inert gas gap to maintain target temperatures and prevent direct coolant contact. Arranged in three pairs, each set will receive a different radiation dose (low, medium, high), with the central pair receiving ~0.1 displacements per atom (dpa)—the expected dose for peak bowing. The experiment will determine the dose-dependent bowing behavior and validate predictive finite element models. In this work, the tubes are freely suspended from pins to allow for unconstrained bowing; however, we present a concept for introducing localized constraints to represent grid spacer effects. Post-irradiation examination will include swelling measurements and profilometry to assess bowing and compare the results with model predictions. This work aims to confirm the conditions under which maximum bowing occurs so as to improve the reliability of SiC/SiC performance models in nuclear applications.

Russell, Nick [ORNL] (ORCID:0000000296099820)↗

DOE OSTI2013

Global latitudinal-asymmetric vegetation growth trends and their driving mechanisms: 1982-2009

Using a recent Leaf Area Index (LAI) dataset and the Community Land Model version 4 (CLM4), we investigate percent changes and controlling factors of global vegetation growth for the period 1982 to 2009. Over that 28-year period, both the remote-sensing estimate and model simulation show a significant increasing trend in annual vegetation growth. Latitudinal asymmetry appeared in both products, with small increases in the Southern Hemisphere (SH) and larger increases at high latitudes in the Northern Hemisphere (NH). The south-to-north asymmetric land surface warming was assessed to be the principal driver of this latitudinal asymmetry of LAI trend. Heterogeneous precipitation functioned to decrease this latitudinal LAI gradient, and considerably regulated the local LAI change. CO2 fertilization during the last three decades, was simulated to be the dominant cause for the enhanced vegetation growth. Our study, though limited by observational and modeling uncertainties, adds further insight into vegetation growth trends and environmental correlations. These validation exercises also provide new quantitative and objective metrics for evaluation of land ecosystem process models at multiple spatio-temporal scales.

Mao, Jiafu↗

DOE OSTI2019

Spatial gradients in the characteristics of soil-carbon fractions are associated with abiotic features but not microbial communities

Coastal terrestrial-aquatic interfaces (TAIs) are dynamic zones of carbon biogeochemical cycling influenced by landscape-scale salinity gradients. However, there is significant cross-system heterogeneity in salinity influences over TAI soil biogeochemical function, such as gas fluxes and dissolved organic carbon concentrations. One potential explanation for the lack of consistency is unrecognized variation in soil organic carbon chemistry and microbial communities. To provide a foundation for understanding the role of C chemistry and microbial communities on biogeochemical cycling in coastal soils, we evaluated associations among a landscape scale soil salinity gradient, molecular-level soil carbon chemistry, and microbial community assembly processes. Based on observations in related systems we hypothesized systematic shifts along the salinity gradient in organic carbon thermodynamics, biochemical transformation profiles, and heteroatom content. While the thermodynamic favorability of water soluble organic compounds in shallow soils did decrease with increasing salinity (R2=0.78, p<0.00001), this was opposite the expected direction. The number of inferred biochemical transformations and total heteroatom content also decreased with increasing salinity. These patterns indicate lower microbial activity at higher salinity that is potentially constrained by less favorable organic carbon. Furthermore, organic compound classes appeared to be primarily marine/algal derived in the floodplains with greater terrestrial signatures in upland soil. Based on a recent simulation-based study we further hypothesized a relationship between microbial community assembly processes and carbon chemistry. Null modelling revealed particularly strong influences of dispersal limitation over microbial composition, which may be due limited hydrologic connectivity within the clay-rich soils. Signatures of dispersal limitation indicate stochastically assembled microbial communities, which was further reflected in the lack of an association between community assembly processes and carbon chemistry. This suggests a disconnect between microbial composition and carbon biogeochemistry, thereby indicating that the salinity-associated gradient in carbon chemistry was driven by a combination of spatially-structured inputs and salinity-associated metabolic responses of microbial communities that are independent of community composition. Our results collectively indicate that soil carbon chemistry varies systematically with salinity at the landscape scale. This variation has not been accounted for in studies examining impacts of salinity to biogeochemical rates and pools. We propose that impacts of salinity in coastal soil biogeochemistry need to be understood in context of carbon chemistry and microbial physiology, while microbial composition may have less influence.

soil carbon, Coastal Biogeochemistry, microbial co↗

DOE OSTI2022

Towards 2+1 Flavor Lattice QCD Results for the Heavy Quark Diffusion Coefficient

We apply and extend a novel approach to non-perturbatively estimate the heavy-quark momentum diffusion coefficient κ, which is a key input for the theoretical description of heavy quarkonium production in heavy ion collisions, and is important for the understanding of the elliptic flow and nuclear suppression factor of heavy flavor hadrons. In the heavy-quark limit, this coefficient is encoded in the spectral functions of color-electric and color-magnetic correlators that we calculate on the lattice to high pre cision by applying gradient flow. In a recent study we have considered quenched QCD at 1.5 T c , where we performed a detailed study of the lat tice spacing and flow time dependence of the color-electric correlator, and, using theoretically well-established model fits for the spectral reconstruc tion, we estimated the heavy-quark diffusion coefficient. Equipped with the experience obtained in quenched QCD, we estimate $κ$ from 2+1 flavor QCD ensembles at small but finite lattice spacing and flow time without increasing systematic errors significantly.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1993

Why convective heat transport in the solar nebula was inefficient

The radial distributions of the effective temperatures of circumstellar disks associated with pre-main sequence (T Tauri) stars are relatively well-constrained by ground-based and spacecraft infrared photometry and radio continuum observations. If the mechanisms by which energy is transported vertically in the disks are understood, these data can be used to constrain models of the thermal structure and evolution of solar nebula. Several studies of the evolution of the solar nebula have included the calculation of the vertical transport of heat by convection. Such calculations rely on a mixing length theory of transport and some assumption regarding the vertical distribution of internal dissipation. In all cases, the results of these calculations indicate that transport by radiation dominates that by convection, even when the nebula is convectively unstable. A simple argument that demonstrates the generality (and limits) of this result, regardless of the details of mixing length theory or the precise distribution of internal heating is presented. It is based on the idea that the radiative gradient in an optically thick nebula generally does not greatly exceed the adiabatic gradient.

Cassen, P.↗

NASA NTRS1994

Escarpment evolution on high-elevation rifted margins: Insights derived from a surface processes model that combines diffusion, advection and reaction

Experiments with a surface processes model of large-scale (1-1000 km) long-term (1-100 m.y.) erosional denudation are used to establish the controls on the evolution of a model escarpment that is related to the rifting of a continent. The mdoel describes changes in topographic form as a result of sumultaneous short- and long-range mass transport representing hillslope (diffusive) processes and fluvial transport (advection), repsectively. Fluvial entrainment is modeled as a first-order kinetic reaction which reflects the erodibility of the substrate, and therefore the fluvial system is not necessarily carrying at capacity. One dimensional and planform models demonstrate that the principal controls on the evolution of an initially steep model escarpment are (1) antecedent topography/drainage; (2) the timesale (or equivalently a length scale) in the fluvial entrainment reaction; (3) the flexural response of the lithosphere to denudation; and (4) the relative efficiencies of the short- and long-range transport processes. When rainfall and substrate lithology are uniform, a significant amount of discharge draining over the escarpment top causes it to degrade. Only when the top of the model escarpment coincides with a drainage divide can escarpment retreat occur for these conditions. An additional requirement for retreat of a model escarpment without decline is a long reaction time scale for fluvial entrainment. This corresponds to a substrate that is hard to detach by flucial erosion, and therefore to fluvial erosion that is not transport limited. Coninuous backtilting of an escarpment due ot flexural isostatic uplift in response to denudational unloading helps maintain the scarp top as a divide. It is essntial if the escarpment gradient is to be preserved during retreat in a uniform lithology. Low flexural rigidieties propote steep and slowly retreating escarpments. For given rainfall and substrate conditions, the morphology of a retraeating model escarpment is determined by the ratio of hte short-range diffusive and long range advective transport efficiencies. A low ration (which is interpreted to correspond to a more humid, temperate climate) produces a convex upper slope, and concave lower slope morphology and only major escarpments are predicted to preserve a high scarp gradient. Lithological contrasts in hte model produce more complex morphologies and predict the formation of scarps crowned by an erosionaly resistant caprock. However, resistant caprocks are not an essential requirement for model scarps to retreat. We conclude that the inferred controls and model behavior are both consistent with the present-day morphology of rifted continental margins and with modern conceptual models of landscape evolution.

Kooi, Henk↗

NASA NTRS2003

Oxidation of Carbon Fibers in a Cracked Ceramic Matrix Composite Modeled as a Function of Temperature

The oxidation model simulates the oxidation of the reinforcing carbon fibers within a ceramic matrix composite material containing as-fabricated microcracks. The physics-based oxidation model uses theoretically and experimentally determined variables as input for the model. The model simulates the ingress of oxygen through microcracks into a two-dimensional plane within the composite material. Model input includes temperature, oxygen concentration, the reaction rate constant, the diffusion coefficient, and the crack opening width as a function of the mechanical and thermal loads. The model is run in an iterative process for a two-dimensional grid system in which oxygen diffuses through the porous and cracked regions of the material and reacts with carbon in short time steps. The model allows the local oxygen concentrations and carbon volumes from the edge to the interior of the composite to be determined over time. Oxidation damage predicted by the model was compared with that observed from microstructural analysis of experimentally tested composite material to validate the model for two temperatures of interest. When the model is run for low-temperature conditions, the kinetics are reaction controlled. Carbon and oxygen reactions occur relatively slowly. Therefore, oxygen can bypass the carbon near the outer edge and diffuse into the interior so that it saturates the entire composite at relatively high concentrations. The kinetics are limited by the reaction rate between carbon and oxygen. This results in an interior that has high local concentrations of oxygen and a similar amount of consumed carbon throughout the cross section. When the model is run for high-temperature conditions, the kinetics are diffusion controlled. Carbon and oxygen reactions occur very quickly. The carbon consumes oxygen as soon as it is supplied. The kinetics are limited by the relatively slow rate at which oxygen is supplied in comparison to the relatively fast rate at which carbon and oxygen reactions occur. This results in a sharp gradient in oxygen concentration from the edge where it is supplied to the nearest source of carbon, which is where the oxygen is quickly consumed. A moving reaction front is seen in which the outlaying carbon is consumed before the next inner layer of carbon begins to react.

Halbig, Michael C.↗

NASA NTRS2016

Hybrid Differential Dynamic Programming with Stochastic Search

Differential dynamic programming (DDP) has been demonstrated as a viable approach to low-thrust trajectory optimization, namely with the recent success of NASA's Dawn mission. The Dawn trajectory was designed with the DDP-based Static/Dynamic Optimal Control algorithm used in the Mystic software.1 Another recently developed method, Hybrid Differential Dynamic Programming (HDDP),2, 3 is a variant of the standard DDP formulation that leverages both first-order and second-order state transition matrices in addition to nonlinear programming (NLP) techniques. Areas of improvement over standard DDP include constraint handling, convergence properties, continuous dynamics, and multi-phase capability. DDP is a gradient based method and will converge to a solution nearby an initial guess. In this study, monotonic basin hopping (MBH) is employed as a stochastic search method to overcome this limitation, by augmenting the HDDP algorithm for a wider search of the solution space.

Aziz, Jonathan↗

NASA NTRS2016

Hybrid Differential Dynamic Programming with Stochastic Search

Differential dynamic programming (DDP) has been demonstrated as a viable approach to low-thrust trajectory optimization, namely with the recent success of NASAs Dawn mission. The Dawn trajectory was designed with the DDP-based Static Dynamic Optimal Control algorithm used in the Mystic software. Another recently developed method, Hybrid Differential Dynamic Programming (HDDP) is a variant of the standard DDP formulation that leverages both first-order and second-order state transition matrices in addition to nonlinear programming (NLP) techniques. Areas of improvement over standard DDP include constraint handling, convergence properties, continuous dynamics, and multi-phase capability. DDP is a gradient based method and will converge to a solution nearby an initial guess. In this study, monotonic basin hopping (MBH) is employed as a stochastic search method to overcome this limitation, by augmenting the HDDP algorithm for a wider search of the solution space.

Aziz, Jonathan↗