SEARCH · Engineering Papers

Results for “three phase modeling”

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 469 records · Page 26

All sources NASA NTRS DOE OSTI

NASA NTRS1997

[Low-Frequency Flow Oscillation

The results of the research conducted under this grant are presented in detail in three Master theses, by Heinrich, Balow, and Broeren. Additional analysis of the experimental data can be found in two AIAA Journal articles and two conference papers. Citations for all of the studies' publications can be found in the bibliography which is attached. The objective of Heinrich's study was to document the low-frequency flow oscillation on the LRN-1007 airfoil, which had been previously observed at low Reynolds number, to determine its origin, and explore the phenomenon at higher Reynolds number. Heinrich performed detailed flow visualization on the airfoil using surface fluorescent oil and laser-sheet off-body visualization. A large leading-edge separation bubble and trailing-edge separation was identified on the airfoil just prior to the onset of the unsteady stall flow oscillation. From the laser-sheet data, the unsteady flow appeared as a massive boundary-layer separation followed by flow reattachment. Hot-wire data were taken in the wake to identify the presence of the flow oscillation and the dominant frequency. The oscillation was found in the flow from a Reynolds number of 0.3 to 1.3 x 10 exp 6. The Strouhal number based on airfoil projected height was nominally 0.02 and increased slightly with increasing Reynolds number and significantly with increasing airfoil angle of attack. Balow focused his research on the leading-edge separation bubble which was hypothesized to be the origin of the low-frequency oscillation. Initially, experimental measurements in the bubble at the onset of the low-frequency oscillation were attempted to study the characteristics of the bubble and explain possible relationships to the shear-layer-flapping phenomena. Unfortunately, the bubble proved to be extremely sensitive to the probe interference and it drastically reduced the size of the bubble. These detailed measurements were then abandoned by Balow. However, this led to a series of tests where the leading-edge bubble and trailing-edge separation were altered and the affect on the flow-oscillation studied. Balow found that by tripping the airfoil boundary-layer with "zigzag" tape ahead of bubble separation, the bubble was effectively eliminated mid the oscillation suppressed. Wake survey drag measurements showed a drastic reduction in airfoil drag when the bubble and oscillation were eliminated. Using the "zigzag" tape, the trailing-edge separation was moved downstream approximately 5 percent chord. This was found to reduce the amplitude of the oscillation, particularly in the onset stage at low angle of attack (around 14 degrees). Through detailed analysis of the wake behind the airfoil during the unsteady flow oscillation, Balow provided a better understanding of the wake flowfield. Broeren studied the oscillating flowfield in detail at Reynolds number equal 3 x 10 exp 5 and an angle of attack of 15 degrees using laser Doppler velocimetry (LDV). Two-dimensional LDV data were acquired at 687 grid points above the model upper surface while hot-wire data were taken simultaneously in the wake. Using the hot-wire signal, the LDV data were phase averaged into 24 bins to represent a single ensemble average of one oscillation cycle. The velocity data showed a flowfield oscillation that could be divided into three flow regimes. In the first regime, the flow over the airfoil was completely separated initially, the flowfield reattached from the leading edge and the reattachment point moved downstream with increasing time or phase. Broeren referred to this as the reattachment regime. The bubble development regime followed, where a leading-edge separation bubble formed at the leading edge and grew with increasing time. During the initial part of this regime the trailing-edge separation continued to move downstream. However, during the last 30 degrees of phase the trailing-edge separation moved rapidly forward and appeared to merge with the leading-edge bubble. During the third regime, the separation regime, the flow was segmented from the airfoil leading edge and did not reattach to the airfoil surface. The reverse flow was seen to grow in vertical extent up from the model surface as the phase increased. Next reattachment began again at the leading edge signaling the start of the reattachment regime, and so the cycle continued. From Broeren's work, the details of the unsteady flowfield over the airfoil were seen for the first time. From this research a great deal has been learned about the low-frequency flow oscillation which naturally occurs on the LRN-1007 airfoil near stall. The oscillation was seen to persist at higher Reynolds number, the dependence of the Strouhal number on angle of attack and Reynolds number were discovered, the critical role played by the laminar bubble was shown and the entire upper surface flowfield during a flow oscillation cycle was measured and analyzed. What still eludes understanding is the scaling of the flow oscillation and why certain airfoils, such as the LRN, have a very strong low-frequency mode and other airfoils exhibit no organized low-frequency oscillation at all.

Bragg, Michael B.↗

NASA NTRS1993

Three-dimensional solidification and melting using magnetic field control

A new two-fluid mathematical model for fully three dimensional steady solidification under the influence of an arbitrary acceleration vector and with or without an arbitrary externally applied steady magnetic field have been formulated and integrated numerically. The model includes Joule heating and allows for separate temperature dependent physical properties within the melt and the solid. Latent heat of phase change during melting/solidification was incorporated using an enthalpy method. Mushy region was automatically captured by varying viscosity orders of magnitude between liquidus and solidus temperature. Computational results were obtained for silicon melt solidification in a parallelepiped container cooled from above and from a side. The results confirm that the magnetic field has a profound influence on the solidifying melt flow field thus changing convective heat transfer through the boundaries and the amount and shape of the solid accrued. This suggests that development of a quick-response algorithm for active control of three dimensional solidification is feasible since it would require low strength magnetic fields.

Dulikravich, George S.↗

NASA NTRS2022

Comparisons Between the Distributions of Dust and Combustion Aerosols in MERRA-2, FLEXPART and CALIPSO and Implications for Deposition Freezing Over Wintertime Siberia

Aerosol distributions have a potentially large influence on climate-relevant cloud properties but can be difficult to observe over the Arctic given pervasive cloudiness, long polar nights, data paucity over remote regions, and periodic diamond dust events that satellites can misclassify as aerosol. We compared Arctic 2008–2015 mineral dust and combustion aerosol distributions from the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) satellite, the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2) reanalysis products, and the FLEXible PARTicle (FLEXPART) dispersion model. Based on coincident, seasonal Atmospheric Infrared Sounder (AIRS) Arctic satellite meteorological data, diamond dust may occur up to 60 % of the time in winter, but it hardly ever occurs in summer. In its absence, MERRA-2 and FLEXPART each predict the vertical and horizontal distribution of large-scale patterns in combustion aerosols with relatively high confidence (Kendall tau rank correlation > 0.6), although a sizable amount of variability is still unaccounted for. They do the same for dust, except in conditions conducive to diamond dust formation where CALIPSO is likely misclassifying diamond dust as mineral dust and near the surface (< ∼ 2 km) where FLEXPART may be overpredicting local dust emissions. Comparisons to ground data suggest that MERRA-2 Arctic dust concentrations can be improved by the addition of local dust sources. All three products predicted that wintertime dust and combustion aerosols occur most frequently over the same Siberian regions where diamond dust is most common in the winter. This suggests that dust aerosol impacts on ice phase processes may be particularly high over Siberia, although further wintertime model validation with non-CALIPSO observations is needed. This assessment paves the way for applying the model-based aerosol simulations to a range of regional-scale Arctic aerosol–cloud interaction studies with greater confidence.

MERRA-2↗

NASA NTRS1979

Negative streamer development in FEP teflon

A computational model is developed which describes the evolution and propagation of an ionizing front (negative streamer) in solid materials. The ionization front consists of drifting avalanching electrons moving self-consistently under the influence of their own space-charge field together with an applied external field. The required input information for the model consists of the functional dependence of the macroscopic transport coefficients on the local electric field, the initial conditions for beginning the calculation, and the strength of the applied field. A computational approach for specifying the transport coefficients and initional conditions is also described. The approach has been implemented by constructing three computer codes which sequentially interface, beginning with single electron scattering, and ending with streamer development. Computational results are presented for model calculations in Teflon. The overall model is perceived to provide a picture of the initiation phase of a propagating discharge in electron-irradiated dielectrics.

Beers, B. L.↗

NASA NTRS2004

Microscopic Theory and Simulation of Quantum-Well Intersubband Absorption

We study the linear intersubband absorption spectra of a 15 nm InAs quantum well using the intersubband semiconductor Bloch equations with a three-subband model and a constant dephasing rate. We demonstrate the evolution of intersubband absorption spectral line shape as a function of temperature and electron density. Through a detailed examination of various contributions, such as the phase space filling effects, the Coulomb many-body effects and the non-parabolicity effect, we illuminate the underlying physics that shapes the spectra. Keywords: Intersubband transition, linear absorption, semiconductor heterostructure, InAs quantum well

Li, Jianzhong↗

NASA NTRS2006

Direct Numerical Simulation of Transitional Multicomponent-Species Gaseous and Multicomponent-Liquid Drop-Laden Mixing

A model of multicomponent-liquid (MC-liquid) drop evaporation in a three-dimensional mixing layer is here exercised at larger Reynolds numbers than in a previous study, and transitional states are obtained. The gas phase is followed in an Eulerian frame and the multitude of drops is described in a Lagrangian frame. Complete coupling between phases is included with source terms in the gas conservation equations accounting for the drop/flow interaction in terms of drop drag, drop heating and species evaporation. The liquid composition, initially specified as a single-Gamma (SG) probability distribution function (PDF) depending on the molar mass is allowed to evolve into a linear combination of two SGPDFs, called the double-Gamma PDF (DGPDF). The compositions of liquid and vapor emanating from the drops are calculated through four moments of the DGPDFs, which are drop-specific and location-specific, respectively. The mixing layer is initially excited to promote the double pairing of its four initial spanwise vortices into an ultimate vortex in which small scales proliferate. Simulations are performed for four liquids of different compositions and the effect of the initial mass loading and initial free-stream gas temperature are explored. For reference, Simulations are also performed for gaseous multicomponent mixing layers for which the effect of Reynolds number is investigated. The results encompass examination of the global layer characteristics, flow visualizations and homogeneous-plane statistics at transition. Comparisons are performed with previous pre-transitional MC-liquid simulations and with transitional single-component (SC) liquid studies. It is found that MCC flows at transition, the classical energy cascade is of similar strength, but that the smallest scales contain orders of magnitude less energy than SC flows, which is confirmed by the larger viscous dissipation in the former case. Contrasting to pre-transitional MC flows, the vorticity and drop organization depend on the initial gas temperature, this being due to the drop/turbulence coupling. The vapor-composition mean molar mass and standard deviation distributions strongly correlate with the initial liquid-composition PDF; such a correlation only exists for the magnitude of the mean but not for that of the standard deviation. Unlike in pre-transitional situations, regions of large composition standard deviation no longer necessarily coincide with regions of large mean molar mass. The kinetic energy, rotational and composition characteristics, and dissipation are liquid specific and the variation among liquids is amplified with increasing free-stream gas temperature. Eulerian and Lagrangian statistics of gas-phase quantities show that the different. Observation framework may affect the perception of the flow characteristics. The gas composition, of which the first four moments are calculated, is shown to be close to, but distinct from a SGPDF. The PDF of the scalar dissipation rate is calculated for drop-laden layers and is shown to depart more significantly from the typically assumed Gaussian in gaseous flows than experimentally measured gaseous scalar dissipation rates, this being attributed to the increased heterogeneity due to drop/flow interactions.

tubulent multicompenent-species↗

NASA NTRS2017

CFD Evaluation of a 3rd Generation LDI Combustor

An effort was undertaken to perform CFD analysis of fluid flow in Lean-Direct Injection (LDI) combustors with axial swirl-venturi elements for next-generation LDI-3 combustor design. The National Combustion Code (NCC) was used to perform non-reacting and two-phase reacting flow computations for a nineteen-element injector array arranged in a three-module, 7-5-7 element configuration. All computations were performed with a consistent approach of mesh-optimization, spray-modeling, ignition and kinetics-modeling with the NCC. Computational predictions of the aerodynamics of the injector were used to arrive at an optimal injector design that meets effective area and fuel-air mixing criteria. LDI-3 emissions (EINOx, EICO and UHC) were compared with the previous generation LDI-2 combustor experimental data at representative engine cycle conditions.

injectors↗

NASA NTRS2017

The Substorm Cycle as Reproduced by Global MHD Models

Recently, Gordeev et al. (2015) suggested a method to test global MHD models against statistical empirical data. They showed that four community-available global MHD models supported by the Community Coordinated Modeling Center (CCMC) produce a reasonable agreement with reality for those key parameters (the magnetospheric size, magnetic field, and pressure) that are directly related to the large-scale equilibria in the outer magnetosphere. Based on the same set of simulation runs, here we investigate how the models reproduce the global loading-unloading cycle. We found that in terms of global magnetic flux transport, three examined CCMC models display systematically different response to idealized2 h north then 2 h south interplanetary magnetic field (IMF) Bz variation. The LFM model shows a depressed return convection and high loading rate during the growth phase as well as enhanced return convection and high unloading rate during the expansion phase, with the amount of loaded unloaded magnetotail flux and the growth phase duration being the closest to their observed empirical values during isolated substorms. Two other models exhibit drastically different behavior. In the BATS-R-US model the plasma sheet convection shows a smooth transition to the steady convection regime after the IMF southward turning. In the Open GGCM a weak plasma sheet convection has comparable intensities during both the growth phase and the following slow unloading phase. We also demonstrate potential technical problem in the publicly available simulations which is related to post processing interpolation and could affect the accuracy of magnetic field tracing and of other related procedures.

Gordeev, E.↗

NASA NTRS2020

Airborne Science Mission Capabilities of the NASA DC-8 and Possible Alternative Aircraft

Age-related concerns associated with operating the 50-year-old NASA Airborne Science Program DC-8 Flying Laboratory aircraft are raising concerns about the prospects of someday retiring the NASA DC-8. Given that eventuality, this study examines how current and anticipated science requirements can be incorporated with aeronautical performance and cost metrics for assessing, comparing, and down-selecting candidate replacement aircraft. A literature review of NASA DC-8 missions and interviews with current DC-8 science users provides context for the initial phase of this study and offers a pathway for collecting aircraft specifications. Performance models are developed and calibrated against published range-payload data for various mission fuel and payload weight combinations. The models are extended to three airborne science missions, with payload weights ranging from 30,000-52,000 lb and mission ranges up to 5,050 nmi. Results point towards the 767-200ER aircraft as a top contender offering a bit more experimenter floor-space, payload capacity, and range than the NASA DC-8, and to the U.S. Navy P-8 aircraft as a viable option offering perhaps not too much less of those items, with the relative benefit of using less fuel per mission.

Thomas A Ozoroski↗

NASA NTRS2019

NASA Centennial Challenge: 3D Printed Habitat, Phase 3 Final Results

NASA's Centennial Challenges program uses prize competitions with the goal of accelerating innovation in the aerospace industry. Competitions in the Centennial Challenges portfolio have previously focused on advancements in space robotics, regolith excavation, bio-printing, astronaut suit design, small satellites, and solar-powered vehicles. NASA's Three Dimensional (3D) Printed Habitat Centennial Challenge represents a partnership between NASA and the non-profit partner: Bradley University, with co-sponsors Caterpillar, Bechtel, Brick and Mortar Ventures, the American Concrete Institute, and the United States Army Corps of Engineers (USACE) Engineer Research and Development Center (ERDC) to spur development in automated additive construction technologies. The challenge asks teams to design and construct a scaled and simulated Martian habitat using indigenous materials and large scale 3D automated printing systems. Phase 1 of the competition, held in 2015, was an architectural design competition for habitat concepts that could be 3D printed. Phase 2, completed in 2017, asked teams to develop feedstocks from indigenous materials and hydrocarbon polymer recyclables, and demonstrate automated printing systems to manufacture these feedstocks into test specimens to assess mechanical strength. This paper will discuss the Phase 3 competition, focusing on technology outcomes that can potentially be infused into both terrestrial and planetary construction applications. The Phase 3 competition was divided into two sub-competitions: 1) virtual construction, where teams created a high fidelity building information model (BIM) of their 3D-printed habitat design and 2) the construction competition, which required teams to 3D print a structural foundation and subject materials samples to freeze/thaw testing and impact testing (level 1), produce a habitat element and complete a hydrostatic test (level 2), and additively manufacture a 1:3 scale habitat onsite in a head to head competition at Caterpillar, inc.'s Edwards Demonstration & Learning Center near Peoria, Illinois over the course of three days (level 3). While the Phase 2 competition focused primarily on the development of novel feedstocks and robotic printing systems, Phase 3 emphasized the scale-up of these systems and autonomous operation (demonstrating the capability to operate systems on precursor missions prior to the arrival of crew, or terrestrially in field operation settings where human tending of a manufacturing system may be limited). The Phase 3 virtual construction levels yielded a number of novel habitat designs, including both modular habitats and vertically-oriented habitat concepts. The Phase 3 construction competition also challenged teams to autonomously place penetrations and interfacing elements in additively manufactured structures. The paper will emphasize potential applications for the new materials and technologies developed under the umbrella of the competition within NASA's portfolio and in Earth-based applications such as disaster response and infrastructure improvement.

Construction↗

NASA NTRS2019

Nasa Centennial Challenge: Three Dimensional (3d) Printed Habitat, Phase 3

Centennial challenge↗

NASA NTRS2012

Modeling Phase-Aligned Gamma-Ray and Radio Millisecond Pulsar Light Curves

Since the discovery of the first eight gamma-ray millisecond pulsars (MSPs) by the Fermi Large Area Telescope, this population has been steadily expanding. Four of the more recent detections, PSR J00340534, PSR J1939+2134 (B1937+21; the first MSP ever discovered), PSR J1959+2048 (B1957+20; the first discovery of a black widow system), and PSR J2214+3000, exhibit a phenomenon not present in the original discoveries: nearly phase-aligned radio and gamma-ray light curves (LCs). To account for the phase alignment, we explore models where both the radio and gamma-ray emission originate either in the outer magnetosphere near the light cylinder or near the polar caps. Using a Markov Chain Monte Carlo technique to search for best-fit model parameters, we obtain reasonable LC fits for the first three of these MSPs in the context of altitude-limited outer gap (alOG) and two-pole caustic (alTPC) geometries (for both gamma-ray and radio emission). These models differ from the standard outer gap (OG)/two-pole caustic (TPC) models in two respects: the radio emission originates in caustics at relatively high altitudes compared to the usual conal radio beams, and we allow both the minimum and maximum altitudes of the gamma-ray and radio emission regions to vary within a limited range (excluding the minimum gamma-ray altitude of the alTPC model, which is kept constant at the stellar radius, and that of the alOG model, which is set to the position-dependent null charge surface altitude). Alternatively, phase-aligned solutions also exist for emission originating near the stellar surface in a slot gap scenario (low-altitude slot gap (laSG) models). We find that the alTPC models provide slightly better LC fits than the alOG models, and both of these give better fits than the laSG models (for the limited range of parameters considered in the case of the laSG models). Thus, our fits imply that the phase-aligned LCs are likely of caustic origin, produced in the outer magnetosphere, and that the radio emission for these pulsars may come from close to the light cylinder. In addition, we were able to constrain the minimum and maximum emission altitudes with typical uncertainties of 30% of the light cylinder radius. Our results therefore describe a third gamma-ray MSP subclass, in addition to the two previously found by Venter et al.: those with LCs fit by standard OG/TPC models and those with LCs fit by pair-starved polar cap models.

Venter, C.↗

NASA NTRS2012

Acoustic Surveys of a Scaled-Model CESTOL Transport Aircraft in Static and Forward Speed Conditions

An 11% scale-model of a Cruise-Efficient Short Take-off and Landing (CESTOL) scalemodel test was recently completed. The test was conducted in the AEDC National Full-Scale Aerodynamic Complex (NFAC) 40- by 80-Foot Wind Tunnel at NASA Ames Research Center. The model included two over-wing pod-mounted turbine propulsion simulators (TPS). The hybrid blended wing-body used a circulation control wing (CCW) with leadingand trailing-edge blowing. The bulk of the test matrix included three forward velocities (40 kts, 60 kts, and 100kts), angle-of-attack variation between -5 and 25 , and CCW mass flow variation. Seven strut-mounted microphones outboard of the left wing provided source directivity. A phased microphone array was mounted outboard of the right wing for source location. The goal of this paper is to provide a preliminary look at the acoustic data acquired during the Advanced Model for Extreme Lift and Improved Aeroacoustics (AMELIA) test for 0 angle-of-attack and 0 sideslip conditions. Data presented provides a good overview of the test conditions and the signal-to-noise quality of the data. TPS height variation showed a difference of 2 dB to 3 dB due to wing shielding. Variation of slot mass flow showed increases of 12 dB to 26 dB above the airframe noise and the TPS increased the overall levels an additional 5 dB to 10 dB.

Burnside, Nathan↗

NASA NTRS1981

Information and display requirements for aircraft terrain following

The display design procedure for manned vehicle systems, is applied and validated, for a particular scenario. The scenario chosen is that of zero visibility high speed terrain following (V = 466 ft/sec, H = 200 ft) with an A-10 aircraft. The longitudal (linearized) dynamics are considered. The variations in (command path over) terrain pi(t) are modeled as a third order random process. The display design methodology is based on the optimal control model of pilot response, and employs this model in various ways in different phases of the design process. The overall methodology indicates that the design process is intended as a precursor to manned simulation. It provides a rank ordering of candidate displays through a three level process.

Kleinman, D. L.↗

NASA NTRS1985

Modeling and analysis of pinhole occulter experiment: Initial study phase

The feasibility of using a generic simulation, TREETOPS, to simulate the Pinhole/Occulter Facility (P/OF) to be tested on the space shuttle was demonstrated. The baseline control system was used to determine the pointing performance of the P/OF. The task included modeling the structure as a three body problem (shuttle-instrument pointing system- P/OP) including the flexibility of the 32 meter P/OF boom. Modeling of sensors, actuators, and control algorithms was also required. Detailed mathematical models for the structure, sensors, and actuators are presented, as well as the control algorithm and corresponding design procedure. Closed loop performance using this controller and computer listings for the simulator are also given.

Vandervoort, R. J.↗

NASA NTRS1984

Observations of electron beams in the low-latitude boundary layer

ISEE 1 electron spectrometer observations are used to study the electron distribution function in the low latitude boundary layer (LLBL) of the magnetosphere in the sunward direction. The study included 43 examples of well-defined layers. Using three-dimensional, model-independent distribution functions of electrons observed in the boundary layer, beams of field-aligned electrons with typically 100 eV energy maximums, streaming in either one, or both, directions have been identified. The phase densities in the LLBL and at lower altitudes are similar, and the beams observed in the LLBL are identified with those observed at lower altitudes. They form an important component of the plasma electrons in the LLBL, which otherwise resembles a mixture of magnetosheath-like and magnetospheric electrons.

Ogilvie, K. W.↗

NASA NTRS1992

Interplanetary program to optimize simulated trajectories (IPOST). Volume 4: Sample cases

The Interplanetary Program to Optimize Simulated Trajectories (IPOST) is intended to support many analysis phases, from early interplanetary feasibility studies through spacecraft development and operations. The IPOST output provides information for sizing and understanding mission impacts related to propulsion, guidance, communications, sensor/actuators, payload, and other dynamic and geometric environments. IPOST models three degree of freedom trajectory events, such as launch/ascent, orbital coast, propulsive maneuvering (impulsive and finite burn), gravity assist, and atmospheric entry. Trajectory propagation is performed using a choice of Cowell, Encke, Multiconic, Onestep, or Conic methods. The user identifies a desired sequence of trajectory events, and selects which parameters are independent (controls) and dependent (targets), as well as other constraints and the cost function. Targeting and optimization are performed using the Standard NPSOL algorithm. The IPOST structure allows sub-problems within a master optimization problem to aid in the general constrained parameter optimization solution. An alternate optimization method uses implicit simulation and collocation techniques.

Hong, P. E.↗

NASA NTRS1992

Interplanetary Program to Optimize Simulated Trajectories (IPOST). Volume 1: User's guide

IPOST is intended to support many analysis phases, from early interplanetary feasibility studies through spacecraft development and operations. The IPOST output provides information for sizing and understanding mission impacts related to propulsion, guidance, communications, sensor/actuators, payload, and other dynamic and geometric environments. IPOST models three degree of freedom trajectory events, such as launch/ascent, orbital coast, propulsive maneuvering (impulsive and finite burn), gravity assist, and atmospheric entry. Trajectory propagation is performed using a choice of Cowell, Encke, Multiconic, Onestep, or Conic methods. The user identifies a desired sequence fo trajectory events, and selects which parameters are independent (controls) and dependent (targets), as well as other constraints and the coat function. Targeting and optimization is performed using the Stanford NPSOL algorithm. IPOST structure allows sub-problems within a master optimization problem to aid in the general constrained parameter optimization solution. An alternate optimization method uses implicit simulation and collocation techniques.

Hong, P. E.↗