Integrated beam orientation and scanning-spot optimization in intensity-modulated proton therapy for brain and unilateral head and neck tumors
Not provided.
SEARCH · Engineering Papers
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.
Not provided.
Laser scanning systems have been used for material processing tasks such as welding, cutting, marking, and drilling. However, applications have been limited by the small range of motion and slow speed of the focusing unit, which carries the focusing optics. To overcome these limitations, a dynamic focusing system with a long travel range and high speed is needed. In this study, a dynamic focusing unit for a laser scanning system with a voice coil motor (VCM) mechanism is proposed to enable fast speed and a wide focusing range. The VCM has finer precision and higher speed than conventional step motors and a longer travel range than earlier lead zirconium titanate actuators. The system has a hollow configuration to provide a laser beam path. This also makes it compact and transmission-free and gives it low inertia. The VCM's magnetics are modeled using a permeance model. Its design parameters are determined by optimization using the Broyden–Fletcher–Goldfarb–Shanno method and a sequential quadratic programming algorithm. After the VCM is designed, the dynamic focusing unit is fabricated and assembled. The permeance model is verified by a magnetic finite element method simulation tool, Maxwell 2D and 3D, and by measurement data from a gauss meter. The performance is verified experimentally. The results show a resolution of 0.2 μm and travel range of 16 mm. These are better than those of conventional focusing systems; therefore, this focusing unit can be applied to laser scanning systems for good machining capability.
The linear optimal control problem of an nth-order time-invariant dynamic system with a quadratic performance functional is usually solved by the Hamilton-Jacobi approach. This leads to the solution of the differential matrix Riccati equation with a terminal condition. The bulk of the computation for the optimal control problem is related to the solution of this equation. There are various algorithms in the literature for solving the matrix Riccati equation. However, computational complexities and storage requirements as a function of numbers of state variables, control variables, and sensors are not available for all these algorithms. In this work, the computational complexities and storage requirements for some of these algorithms are given. These expressions show the immensity of the computational requirements of the algorithms in solving the Riccati equation for large-order systems such as the control of highly flexible space structures. The expressions are also needed to compute the speedup and efficiency of any implementation of these algorithms on concurrent machines.
A new approach is discussed for solving dual structural control optimization problems for high-order flexible space structures, where reduced-order structural models are employed and minimum mass designs are sought. For a given initial structural design, a quadratic control cost is minimized subject to a constant-mass constraint. The sensitivity of the optimal control cost with respect to the structural design variables is then determined and used to obtain successive structural redesigns, using a constrained gradient optimization algorithm. This process is repeated until the constrained control cost sensitivity becomes negligible. The minimum mass design is obtained by solving a sequence of neighboring optimal constant mass designs, where the sequence of optimal performance indices has a minimum at the optimal minimum mass design. A numerical example is presented which demonstrates that this new approach effectively addresses the problem of dual optimization for potentially very high-order structures.
A class of systems governed by second order partial differential equations and driven by controllers located at points r sub i, i = 1, ..., k, is considered. A cost functional quadratic in the time derivative of the state and the control is associated with the system. The optimal controller locations are defined as the ones which minimize the maximum of the cost over all possible initial states. An analytical solution to the associated Riccati equation is presented providing a convenient expression for determining the optimal locations.
A formulation is presented for the coupled optimal design of a structural system and its control by defining a composite objective function as a linear combination of two components: a structural objective and a control objective. For the case when the structural objective is a function of the structural design variables only, and when the control objective is represented by the quadratic functional of the response and control energy, one can analytically express the optimal control in terms of any set of admissible structural design variables. The expression for the optimal control is used recursively in an iterative Newton-Raphson search scheme, the goal of which is to determine a corresponding optimal set of structural design variables that minimize the composite objective function. A numerical example is given to illustrate the computational procedure.
Road vehicle lane changes often initiate traffic disturbances and can therefore impact road networks’ energy and time efficiency. Furthermore, unexpected changes in traffic conditions may also render lane changes counterproductive for the lane-changing vehicle. Vehicle-to-vehicle connectivity combined with anticipative control could address these challenges via improved lane change decisions by automated vehicles. In a move toward this objective, receding horizon control cast as a mixed-integer quadratic program is used to plan lane changing and acceleration in a coupled optimization. A long-term pacing module, based on Pontryagin’s minimum principle from optimal control theory, sets terminal and input references for receding horizon control to target a user’s expected travel time. To remove nonlinear vehicle dynamics from the receding horizon controller, lane change commands are passed to a pure pursuit steering module whose response is approximated by a second-order linear model. Here, comparison against a rule-based reactive algorithm in arterial and highway scenarios shows an 8.9%–13.7% reduction in energy consumption and a 5.2%–10.3% reduction in the travel time, along with navigational improvements.
Consideration is given to the selection of velocity feedback gains for individual dampers for the members of a structurally controlled large flexible space structure. The problem is formulated as an optimal output feedback regulator problem, and necessary conditions are derived for minimizing a quadratic performance function. The diagonal nature of the gain matrix is taken into account, along with knowledge of noise covariances. It is pointed out that the method presented offers a systematic approach to the design of a class of controllers for enhancing structural damping, which have significant potential if used in conjunction with a reduced-order optimal controller for rigid-body modes and selected structural modes.
Here, a novel framework is proposed for integrating biomimetic-based advanced control and multi-agent optimization approaches for nonlinear chemical process applications. In particular, a Biologically-Inspired Optimal Control Strategy, denoted as BIO-CS, is combined with multi-agent optimization (MAO) techniques to provide optimal solutions for dynamic systems. In this combined framework, the BIO-CS algorithm employs gradient-based optimal control solvers for the intermediate problems associated with the leader-follower agents’ local interactions. Also, the MAO uses the capabilities of heuristic-based optimization techniques by sharing process information to obtain optimal operating setpoints for the controller considering an overall process objective. The applicability of the proposed method is demonstrated using a nonlinear, multivariable, process model of a fermentation system. Specifically, the optimal operating points are computed by the MAO implementation for setpoint tracking, trajectory tracking and plant-model mismatch scenarios for BIO-CS application. Results of the developed framework are compared to a gradient-based Sequential Quadratic Programming (SQP) technique and a classical proportional-integral (PI) controller in terms of optimization and control studies, respectively. As an additional contribution, BIO-CS is also cast as a model predictive controller (MPC) for the first time and compared to the agent-based BIO-CS approach in terms of computational time and tracking error. Closed-loop control results show up to 46% improvement in tracking performance during transient for the multi-agent BIO-CS when compared to BIO-CS as MPC for additional computational expense. The obtained results illustrate the capabilities of this novel integrated framework including BIO-CS as MPC to achieve desired nonlinear system performance for various scenarios.
Superconducting (SC) magnets for accelerator concepts are often synthesized by numerically optimizing magnetic field waveforms, a process that requires a subsequent solution of a constrained inverse problem to identify suitable SC magnet windings. When the desired field distribution is intuitive, the inverse process is facilitated by seeding preconceived coil distributions into design optimization methods for refinement. With more complex magnetic field distributions, an initial design may be unknown, and topology optimization tools are required to synthesize current distributions without a priori guidance from a subject matter expert. In this work, we develop a constrained inverse Biot-Savart topology optimization methodology that synthesizes optimal distributions of current density in racetrack-like SC coils. The problem structure is exploited through a computationally efficient quadratic programming formulation, and the method is applied to recently published magnetic field waveforms for a recirculating proton phase shifter, a proton therapy gantry, and dipole magnets with sharp field transitions. The method and results herein identify novel winding configurations that can help magnet designers bring accelerator concepts to fruition.
An algorithm has been developed for the analysis and design of controls for nonlinear systems. The technical approach is to use statistical linearization to model the nonlinear dynamics of a system. A covariance analysis is performed to determine the behavior of the dynamical system and a quadratic cost function. Expressions for the cost function and its derivatives are determined so that numerical optimization techniques can be applied to determine optimal feedback laws. The primary application for this report is centered about the design of controls for nominally linear systems but where the controls are saturated or limited by fixed constraints. The analysis is general however and numerical computation requires only that the specific nonlinearity be considered in the analysis.
The design of polymeric membranes for proton or ionic exchange highly depends on the fundamental understanding of the physical and molecular mechanisms that control the formation of the conduction channels. There is an inherent relation between the dynamical structure of the polymeric membrane and the electrostatic forces that drive membrane segregation and proton transport. Here, we used a multi-scale computational approach to analyze the morphology of sulfonated poly(ether ether ketone) membranes at the mesoscale. A self-consistent description of the electrostatic phenomenon was adopted, where discrete polymer chains and a continuum proton field were embedded in a continuum fluid. Brownian dynamics was used for the evolution of the suspended polymer molecules, while a convection-diffusion transport equation, including the Nernst-Planck diffusion mechanism, accounted for the dynamics of the proton concentration field. We varied the polymer concentration, the degree of sulfonation and the level of confinement to find relationships between membrane structure and proton conduction. Our results indicate that the reduced mobility of polymer chains, at concentrations above overlap, and a moderate degree of sulfonation -i.e., 30% - are essential elements for membrane segregation and proton domain connectivity. These conditions also ensure that the membrane structure is not affected by size or by potential gradients. Importantly, our analysis shows that membrane conductivity and current are linearly dependent on polymer concentration and quadratically dependent on the degree of sulfonation. We found that the optimal polymeric membrane design requires a polymer concentration above overlap and a degree of sulfonation around 50%. These conditions promote a dynamical membrane morphology with a constant density of proton channels. Our results and measurements agree with previous experimental works, thereby validating our model and observations.
Optimal linear controller synthesis with minimized trajectory sensitivity achieved by using quadratic cost function
A technique is presented by which one can apply the Minimum Principle of Pontryagin to stochastic optimal control problems formulated around linear systems with Gaussian noises and general cost criteria. Using this technique, the stochastic nature of the problem is suppressed but for two expectation operations, the optimization being essentially deterministic. The technique is applied to systems with quadratic and non-quadratic costs to illustrate its use.
An implicit finite element numerical solution algorithm is derived for the compressible Navier-Stokes equations expressed in generalized coordinates. The theoretical basis utilizes a Galerkin-Weighted Residuals formulation, and extremization of approximation error within the context of a multipole expansion. A von Neumann analysis for a simplified form indicates the algorithm fourth- to sixth-order phase accurate, with third-order dissipation for the elementary linear element construction. Performance is improved for the algorithm constructed using quadratic interpolation. Numerical experiments for shocked duct flows are employed to optimize the several algorithm parameters. Additional numerical solutions validate algorithm accuracy and utility for aerodynamics applications.
An analysis is made of the dynamics of a thin, homogeneous, flexible square plate exposed to solar radiation disturbance. It is shown that only the first three flexible modes of the plate generate a first order net moment about the center of mass and that the solar radiation pressure does not influence the flexible modes of the plate for small-amplitude vibrations. In the absence of control, for a symmetrical homogeneous square platform, the solar radiation induces a steady angular drift about one of the body principal axes. For extremely flexible platforms, nominally oriented in the local horizontal plane, appreciable rigid modal amplitudes can be induced due to solar radiation, even in the presence of both active and passive control. Optimal control law and feedback gain values are obtained using linear quadratic Gaussian methods.
The minimum weight optimization of structural systems, subject to strength and displacement constraints as well as size side constraints, was investigated by the Simultaneous ANalysis and Design (SAND) approach. As an optimizer, the code NPSOL was used which is based on a sequential quadratic programming (SQP) algorithm. The structures were modeled by the finite element method. The finite element related input to NPSOL was automatically generated from the input decks of such standard FEM/optimization codes as NASTRAN or ASTROS, with the stiffness matrices, at present, extracted from the FEM code ANALYZE. In order to avoid ill-conditioned matrices that can be encountered when the global stiffness equations are used as additional nonlinear equality constraints in the SAND approach (with the displacements as additional variables), the matrix displacement method was applied. In this approach, the element stiffness equations are used as constraints instead of the global stiffness equations, in conjunction with the nodal force equilibrium equations. This approach adds the element forces as variables to the system. Since, for complex structures and the associated large and very sparce matrices, the execution times of the optimization code became excessive due to the large number of required constraint gradient evaluations, the Kreisselmeier-Steinhauser function approach was used to decrease the computational effort by reducing the nonlinear equality constraint system to essentially a single combined constraint equation. As the linear equality and inequality constraints require much less computational effort to evaluate, they were kept in their previous form to limit the complexity of the KS function evaluation. To date, the standard three-bar, ten-bar, and 72-bar trusses have been tested. For the standard SAND approach, correct results were obtained for all three trusses although convergence became slower for the 72-bar truss. When the matrix displacement method was used, correct results were still obtained, but the execution times became excessive due to the large number of constraint gradient evaluations required. Using the KS function, the computational effort dropped, but the optimization seemed to become less robust. The investigation of this phenomenon is continuing. As an alternate approach, the code MINOS for the optimization of sparse matrices can be applied to the problem in lieu of the Kreisselmeier-Steinhauser function. This investigation is underway.
The robustness of optimal regulators for nonlinear, deterministic and stochastic, multi-input dynamical systems is studied under the assumption that all state variables can be measured. It is shown that, under mild assumptions, such nonlinear regulators have a guaranteed infinite gain margin; moreover, they have a guaranteed 50 percent gain reduction margin and a 60 degree phase margin, in each feedback channel, provided that the system is linear in the control and the penalty to the control is quadratic, thus extending the well-known properties of LQ regulators to nonlinear optimal designs. These results are also valid for infinite horizon, average cost, stochastic optimal control problems.