Stability analysis and aerodynamic design optimization of euler equations using variational methods

A shape design sensitivity analysis (DSA) and optimization of structural transient dynamics are proposed for the finite deformation elastoplastic materials under impact with a rigid surface. A shape variation of the structure is considered using the material derivative approach in continuum mechanics. Hyperelasticitybased multiplicatively decomposed elastoplasticity is used for the constitutive model. The implicit Newmark time integration scheme is used for the structural dynamics. The design sensitivity equation is solved at each converged time step with the same tangent stiffness matrix as response analysis without iteration. The cost of sensitivity computation is more efficient than the cost of response analysis for the implicit time integration method. The efficiency and the accuracy of the proposed method are shown through the design optimization of a vehicle bumper. *Communicated by D. Tortorelli.

A shape design sensitivity analysis and optimization procedure is proposed using a meshfree method. A CAD tool connection is established to facilitate the seamless integration of structural modeling, simulation, and design optimization. A new CAD-based method for design velocity computation is proposed, applicable to both the meshfree method and the finite element method. A shape design variable can be chosen among the geometric parameters defined in the CAD model in order to support a concurrent engineering environment. A new CAD-based method for updating the discrete model is proposed, such that the meshfree model is automatically updated in a continuous and consistent way during design optimization. The proposed method is applied to 2D and 3D linear elastic problems, as well as to a large deformation hyperelastic problem.

A unified method for continuum shape design sensitivity analysis and optimal design of mechanical components is developed. A domain method of shape design sensitivity analysis that uses the material derivative concept of continuum mechanics is employed. For numerical implementation of shape optimal design, parameterization of the boundary shape of mechanical components is defined and illustrated using a Bezier surface. In shape design problems, nodal points of the finite element model move as the shape changes. A method of automatic regridding to account for shape change has been developed using a design velocity field in the physical domain that obeys the governing equilibrium equations of the elastic solid. For numerical implementation of the continuum shape design sensitivity analysis and automatic regridding, an established finite element analysis code is used. To demonstrate the feasibility of the method developed, shape design optimization of a main engine bearing cap is carried out as an example.

Variational methods (VM) sensitivity analysis, which is the continuous alternative to the discrete sensitivity analysis, is employed to derive the costate (adjoint) equations, the transversality conditions, and the functional sensitivity derivatives. In the derivation of the sensitivity equations, the variational methods use the generalized calculus of variations, in which the variable boundary is considered as the design function. The converged solution of the state equations together with the converged solution of the costate equations are integrated along the domain boundary to uniquely determine the functional sensitivity derivatives with respect to the design function. The determination of the sensitivity derivatives of the performance index or functional entails the coupled solutions of the state and costate equations. As the stable and converged numerical solution of the costate equations with their boundary conditions are a priori unknown, numerical stability analysis is performed on both the state and costate equations. Thereafter, based on the amplification factors obtained by solving the generalized eigenvalue equations, the stability behavior of the costate equations is discussed and compared with the state (Euler) equations. The stability analysis of the costate equations suggests that the converged and stable solution of the costate equations is possible only if the computational domain of the costate equations is transformed to take into account the reverse flow nature of the costate equations. The application of the variational methods to aerodynamic shape optimization problems is demonstrated for internal flow problems at supersonic Mach number range. The study shows, that while maintaining the accuracy of the functional sensitivity derivatives within the reasonable range for engineering prediction purposes, the variational methods show a substantial gain in computational efficiency, i.e., computer time and memory, when compared with the finite difference sensitivity analysis.

<p>The primary objective of this study is to develop new computational methods for robust design optimization (RDO) and reliability-based design optimization (RBDO) of high-dimensional, complex engineering systems. Four major research directions, all anchored in polynomial dimensional decomposition (PDD), have been defined to meet the objective. They involve: (1) development of new sensitivity analysis methods for RDO and RBDO; (2) development of novel optimization methods for solving RDO problems; (3) development of novel optimization methods for solving RBDO problems; and (4) development of a novel scheme and formulation to solve stochastic design optimization problems with both distributional and structural design parameters.</p> <p>The major achievements are as follows. Firstly, three new computational methods were developed for calculating design sensitivities of statistical moments and reliability of high-dimensional complex systems subject to random inputs. The first method represents a novel integration of PDD of a multivariate stochastic response function and score functions, leading to analytical expressions of design sensitivities of the first two moments. The second and third methods, relevant to probability distribution or reliability analysis, exploit two distinct combinations built on PDD: the PDD-SPA method, entailing the saddlepoint approximation (SPA) and score functions; and the PDD-MCS method, utilizing the embedded Monte Carlo simulation (MCS) of the PDD approximation and score functions. For all three methods developed, both the statistical moments or failure probabilities and their design sensitivities are both determined concurrently from a single stochastic analysis or simulation. Secondly, four new methods were developed for RDO of complex engineering systems. The methods involve PDD of a high-dimensional stochastic response for statistical moment analysis, a novel integration of PDD and score functions for calculating the second-moment sensitivities with respect to the design variables, and standard gradient-based optimization algorithms. The methods, depending on how statistical moment and sensitivity analyses are dovetailed with an optimization algorithm, encompass direct, single-step, sequential, and multi-point single-step design processes. Thirdly, two new methods were developed for RBDO of complex engineering systems. The methods involve an adaptive-sparse polynomial dimensional decomposition (AS-PDD) of a high-dimensional stochastic response for reliability analysis, a novel integration of AS-PDD and score functions for calculating the sensitivities of the failure probability with respect to design variables, and standard gradient-based optimization algorithms, resulting in a multi-point, single-step design process. The two methods, depending on how the failure probability and its design sensitivities are evaluated, exploit two distinct combinations built on AS-PDD: the AS-PDD-SPA method, entailing SPA and score functions; and the AS-PDD-MCS method, utilizing the embedded MCS of the AS-PDD approximation and score functions. In addition, a new method, named as the augmented PDD method, was developed for RDO and RBDO subject to mixed design variables, comprising both distributional and structural design variables. The method comprises a new augmented PDD of a high-dimensional stochastic response for statistical moment and reliability analyses; an integration of the augmented PDD, score functions, and finite-difference approximation for calculating the sensitivities of the first two moments and the failure probability with respect to distributional and structural design variables; and standard gradient-based optimization algorithms, leading to a multi-point, single-step design process.</p> <p>The innovative formulations of statistical moment and reliability analysis, design sensitivity analysis, and optimization algorithms have achieved not only highly accurate but also computationally efficient design solutions. Therefore, these new methods are capable of performing industrial-scale design optimization with numerous design variables.</p>

In the direct method of shape optimum design, mesh regridding is a necessary step in avoiding the distortion of finite element grids. The effect of internal-nodal movements on the shape design sensitivity associated with the regridding process has been carefully examined in this paper. It shows that to maintain computational accuracy, the effects of internal-nodal movement should be considered in shape design sensitivity analysis. To take the internal-nodal movement into account, the numerical accuracy of the speed method is investigated. A new approach, which constructs the velocity field over the element domain, is presented. Although it takes more computational time, the new approach does improve the accuracy of shape design sensitivity.

Because of their superior mechanical, structural, and electronic properties, graphene sheets (GSs) are supposed to be components of nanoelectromechanical systems (NEMS). In this study, we carry out shape and structural optimization of GSs with topological defects in natural vibration problem to enhance their fundamental frequencies. At first, we model GSs as continuum frame models based on the molecular mechanics method. Then, we carry out shape design optimization of GSs by using a developed free-form optimization method for frame structures. In the shape design optimization, we use the fundamental frequency as objective function and enhance it considering volume constraint and repeated eigenvalues. Next, we derive the optimal atom structures of GSs with topological defects using a combination of the Phase-Field-Crystal method, the Voronoi tessellation method, and molecular dynamics simulation. The numerical results show that the fundamental frequencies of GSs can be significantly increased according to the shape and structural design optimization, and we can get their optimal atom structures with defects, which is helpful for designing GSs as components in NEMS.

Shape design sensitivity analysis and optimization of three dimensional resonant systems

Shape design sensitivity analysis and optimization of elasto-plasticity with frictional contact

This paper presents results for the study of equivalence between the total form continuum sensitivity equation (CSE) method and the discrete analytic method of shape design sensitivity analysis. For the discrete analytic method, the sensitivity equations are obtained by taking analytic derivatives of the discretized equilibrium equations with respect to the shape design parameters. For the CSE method, the equilibrium equations are firstly differentiated to form a set of linear continuous sensitivity equations and then discretized for solving the shape sensitivities. The sensitivity equations can be derived by taking the material derivatives (total form) or the partial derivatives (local form) of the equilibrium equations. The total form CSE method is shown for the first time to be equivalent, after finite element discretization, to the discrete analytic method for nonlinear second-order differential equations of a particular form with design dependent loads when they use the same: (1) finite element discretization, (2) numerical integration of element matrices, (3) design velocity fields that are linear with respect to the design variable and (4) shape functions for domain transformation and for design velocity field calculations. The shape sensitivity equations for three-dimensional geometric nonlinear elastic structures and linear potential flow are derived by using both total form CSE and discrete analytic method to show the equivalence of the two methods for these specific examples. The accuracy of shape sensitivity analysis is verified by potential flow around an airfoil and a joined beam with an airfoil under gust load. The results show that analytic sensitivity results are consistent with the complex step results.

Using the bond-based peridynamics and the parallel computation with binary decomposition, an adjoint shape design sensitivity analysis(DSA) method is developed for the dynamic crack propagation problems. The peridynamics includes the successive branching of cracks and employs the explicit scheme of time integration. The adjoint variable method is generally not suitable for path-dependent problems but employed since the path of response analysis is readily available. The accuracy of analytical design sensitivity is verified by comparing it with the finite difference one. The finite difference method is susceptible to the amount of design perturbations and could result in inaccurate design sensitivity for highly nonlinear peridynamics problems with respect to the design. It turns out that C1-continuous volume fraction is necessary for the accurate evaluation of shape design sensitivity in peridynamic discretization.

A continuum-based sizing design sensitivity analysis (DSA) method is presented for the transient dynamic response of non-linear structural systems with elastic–plastic material and large deformation. The methodology is aimed for applications in non-linear dynamic problems, such as crashworthiness design. The first-order variations of the energy forms, load form, and kinematic and structural responses with respect to sizing design variables are derived. To obtain design sensitivities, the direct differentiation method and updated Lagrangian formulation are used since they are more appropriate for the path-dependent problems than the adjoint variable method and the total Lagrangian formulation, respectively. The central difference method and the finite element method are used to discretize the temporal and spatial domains, respectively. The Hughes–Liu truss/beam element, Jaumann rate of Cauchy stress, rate of deformation tensor, and Jaumann rate-based incrementally objective stress integration scheme are used to handle the finite strain and rotation. An elastic–plastic material model with combined isotropic/kinematic hardening rule is employed. A key development is to use the radial return algorithm along with the secant iteration method to enforce the consistency condition that prevents the discontinuity of stress sensitivities at the yield point. Numerical results of sizing DSA using DYNA3D yield very good agreement with the finite difference results. Design optimization is carried out using the design sensitivity information. Copyright © 2000 John Wiley & Sons, Ltd.

In this paper, a novel optimal shape design method is proposed using the finite-difference time-domain (FDTD) method and the design sensitivity analysis to obtain broad-band characteristics of microwave devices. In shape design problem, the nodes that describe the shape of geometry to be optimized are taken as design variables. The design sensitivity is evaluated using the adjoint variable equation that is obtained from a terminal-value problem. The adjoint equation can be solved by the FDTD technique with the backward time scheme. With this method, a Ka-band unilateral fin line is tested to show validity.

Shape design optimization of thermoelasticity problems using isogeometric boundary element method

Defining and identifying the aeroacoustic sources in a turbulent flow is a great challenge especially for noise control strategy. The purpose of the present study consists in proposing a new methodology to localize regions associated with sound generation. These regions are associated, in the present work, with those of high sensitivity of the acoustic field, using the heuristic argument that modifying the flow in these regions would lead to a very significant change in the radiated noise. The proposed method relies on the efficient coupling between the time-reversal theory applied to the Euler equations and the complex differentiation method to compute the sensitivity variable. To the knowledge of the authors, this is the first time that the time-reversal technique is applied to vectorial hydrodynamic equations, in place of the classical scalar wave equation. Subsequently, regions associated with sound generation are related to spatiotemporal events which exhibit the maximum of sensitivity to acoustical disturbances measured in far field. The proposed methodology is then successively tested on three cases for which the nature of the source is different: injection of mass, vibrating surfaces and flow instabilities arising in a plane mixing layer flow. For each test case, the two-dimensional Euler equations are solved using a numerical solver based on a pseudo-characteristics formulation. During these computations flow, variables are stored only at the computational boundaries. These variables are time reversed and relevant information concerning the acoustical disturbances is tagged using complex differentiation in order to lead the sensitivity analysis. The same numerical solver is used to access the evolution of the time-reversed variables. In each test case, the proposed methodology allows to localize successfully zones associated with noise generation.