An accelerated method is presented for the iterative condensation of eigenproblems, The present study was motivated by the improved reduction system and the succession-level approximate reduction (SAR), The reduction ...
详细信息
An accelerated method is presented for the iterative condensation of eigenproblems, The present study was motivated by the improved reduction system and the succession-level approximate reduction (SAR), The reduction procedures are supplemented with the secund-order approximation in the series expansion of the system transformation. The reduced equation of an equivalent system and the transformation matrix are updated in an iterative manner. In addition, systematic derivation and comparison of the equations involved in various condensations have been sought. The matrix update incorporates not only inverse iteration but also subspace transformation implicitly. The series expansion can be considered as repeated updates of the transformation matrix through inverse iteration. The solution accuracy is sensitive to the selection of the degrees of freedom, for which sequential elimination or energy method may be preferable. When a poor selection causes a sudden failure of the update method, the hybrid dynamic condensation fan be used. The method of SAR is closely related to the hybrid dynamic condensation.
A Mach 3 adiabatic flat plate turbulent boundary layer is studied using large-eddy simulation (LES), The filtered compressible Navier-Stokes equations are solved on a three-dimensional unstructured grid of tetrahedral...
详细信息
A Mach 3 adiabatic flat plate turbulent boundary layer is studied using large-eddy simulation (LES), The filtered compressible Navier-Stokes equations are solved on a three-dimensional unstructured grid of tetrahedral cells. A compressible extension of the rescaling-reintroducing process of Lund et al. (Lund, T., Wu, X., and Squires, K., "Generation of Turbulent Inflow Data for Spatially-Developing Boundary Layer Simulations," Journal of Compartational Physics, Vol. 140, No. 2, 1998, pp. 233-2583 is developed to generate the inflow conditions. The effect of the subgrid-scale motion is incorporated using two approaches, namely, monotone integrated LES (MILES) and the Smagorinsky subgrid-scale model. A detailed grid refinement study is performed. The statistical predictions (friction velocity, adiabatic wall temperature, mean velocity profile, normal Reynolds stress, and Reynolds shear stress) obtained from MILES are in close agreement viith experimental data and direct numerical simulation. The results indicate that the subgrid-scale effects can be adequately modeled using MILES without the need for the Smagorinsky model.
A continuous adjoint approach is developed to obtain the sensitivity derivatives for the Euler equations, The complete derivation of the costate equations and their transversality (boundary) conditions are presented. ...
详细信息
A continuous adjoint approach is developed to obtain the sensitivity derivatives for the Euler equations, The complete derivation of the costate equations and their transversality (boundary) conditions are presented. Both the state and the costate equations are second-order finite volume discretized fur unstructured meshes, and they are coupled with a constrained optimization algorithm. Also integrated into the overall methodology are a geometry parameterization method for the shape optimization, and a dynamic unstructured mesh method fur the shape evolution and the consequent volume mesh adaptations. For the proof of concept, three transonic airfoil optimization problems are presented. This method accepts general cost functionals, which are not necessarily functions of pressure only. It is also shown that a switch to the natural coordinate system in conjunction with the reduction of the governing state equation to the control surface results in sensitivity integrals that are only a function of the tangential derivatives of the state variables. This approach eliminates the need fur normal derivative computations that can be erroneous.
Provides information on a study which examined an approach to the general solution of laminar boundary-layer flows discovered by Falkner and Skan in 1931. numerical procedure; Results and discussion; Conclusions.
Provides information on a study which examined an approach to the general solution of laminar boundary-layer flows discovered by Falkner and Skan in 1931. numerical procedure; Results and discussion; Conclusions.
We study the behavior of crystals that undergo martensitic transformations. On cooling, the high-temperature phase (austenite) transforms to the martensite phase changing its crystalline symmetry. The lower crystallin...
详细信息
We study the behavior of crystals that undergo martensitic transformations. On cooling, the high-temperature phase (austenite) transforms to the martensite phase changing its crystalline symmetry. The lower crystalline symmetry of the martensite gives rise to several variants of martensite. Each variant has an associated transformation strain. These variants accommodate themselves (according to the boundary conditions) forming a microstructure that minimizes the elastic energy. This minimum value of the energy is called the the effective energy. We assume that all the material is in the martensite phase (i.e. the material is at low temperatures). We show that, assuming the geometrically linear approximation, the maximum of the effective energy restricted to applied strains in the convex hull of the transformation strains is attained by an applied strain that is a convex combination of only two transformation strains. We derive a recurrence relation to compute the energy corresponding to laminated microstructures of arbitrary rank, under the assumption that the variants of martensite are linearly elastic and their elastic moduli are isotropic. We use this recurrence relation to develop an algorithm that minimizes the energy over microstructures that belong to the class of rank-v laminates. We apply our methods to the case in which the transformation is cubic to monoclinic (corresponding to TiNi). We conclude with some comments on the possible implications of our calculations on the behavior of this shape-memory alloy. (C) 2001 Elsevier Science Ltd. All rights reserved.
A numerical algorithm for nonlinear elastic relaxation of a multilayer woven fabric is introduced and tested. The equilibrium solutions are compared with real samples. An excellent result is obtained in spite of two s...
详细信息
A numerical algorithm for nonlinear elastic relaxation of a multilayer woven fabric is introduced and tested. The equilibrium solutions are compared with real samples. An excellent result is obtained in spite of two simplifications: Bending stiffness of the fibers and friction between the fibers are both neglected. The numerical simulation is very fast and cost efficient in the search for optimal fabrics. (C) 2001 American Institute of Physics.
This note presents a numerically reliable approach for computing the normalized-coprime factorization of an arbitrary (possibly improper) rational matrix. We consider both the cases in which the factorization is copri...
详细信息
This note presents a numerically reliable approach for computing the normalized-coprime factorization of an arbitrary (possibly improper) rational matrix. We consider both the cases in which the factorization is coprime and normalized with respect to either the left-half plane or the unit disk, corresponding to systems with continuous- or discrete-time evolutions. The proposed algorithm remedies the numerical drawbacks of alternative methods available thus far in the literature.
dA model is presented for plastic deformation with grain subdivision into parallel bands. The experimental reference is the subdivision into "cell blocks" observed in rolled aluminium. The model maintains in...
详细信息
dA model is presented for plastic deformation with grain subdivision into parallel bands. The experimental reference is the subdivision into "cell blocks" observed in rolled aluminium. The model maintains intragranular strain continuity between the bands with relaxed constraints. One version of the model maintains intergranular strain continuity by imposing identical strains in all grains. Another version does not provide formal fulfilment of intergranular strain continuity, but it tries to minimize strain discontinuity by selection of the appropriate physical solutions. Part I deals with the initial stage of grain subdivision at low strain. Part II (Leffers, T. 2001. A model for rolling deformation with grain subdivision. Part II: the subsequent stage. Int. J. Plasticity 17, 491-511.) deals with the subsequent stages at higher strains and the resulting rolling texture. (C) 2001 Elsevier Science Ltd. All rights reserved.
Part I of the present work dealt with the initial stage of plastic deformation with grain subdivision into two band families. It covered the situation with different average strains in the individual grains (the non-T...
详细信息
Part I of the present work dealt with the initial stage of plastic deformation with grain subdivision into two band families. It covered the situation with different average strains in the individual grains (the non-Taylor case) and the situation with identical average strains in the individual grains (the Taylor case). For the latter situation part I included solutions with 5 plus 3 and with 4 plus 4 active slip systems in the two band families, respectively. In part II we deal with the subsequent stage of plastic deformation (finite strains) for the Taylor case with four active slip systems in each band family. In the subsequent stage the cooperation between the two band families leads to an energetic advantage for grain subdivision, and it leads to a simulated texture with decreased sharpness in accordance with experimental observations. Finally, the physics behind grain subdivision are discussed on the basis of the results in Parts I and II and various general considerations. (C) 2001 Elsevier Science Ltd. All rights reserved.
This paper deals with a thermodynamically consistent numerical formulation for coupled thermoplastic problems including phase-change phenomena and frictional contact. The final goal is to get an accurate, efficient an...
详细信息
This paper deals with a thermodynamically consistent numerical formulation for coupled thermoplastic problems including phase-change phenomena and frictional contact. The final goal is to get an accurate, efficient and robust numerical model, able for the numerical simulation of industrial solidification processes. Some of the current issues addressed in the paper are the following. A fractional step method arising from an operator split of the governing differential equations has been used to solve the nonlinear coupled system of equations, leading to a staggered product formula solution algorithm. Nonlinear stability issues are discussed and isentropic and isothermal operator splits are formulated. Within the isentropic split, a strong operator split design constraint is introduced, by requiring that the elastic and plastic entropy, as well as the phase-change induced elastic entropy due to the latent heat, remain fixed in the mechanical problem. The formulation of the model has been consistently derived within a thermodynamic framework. All the material properties have been considered to be temperature dependent. The constitutive behavior has been defined by a thermoviscous/ elastoplastic free energy function, including a thermal multiphase change contribution. Plastic response has been modeled by a J2 temperature dependent model, including plastic hardening and thermal softening. The constitutive model proposed accounts for a continuous transition between the initial liquid state, the intermediate mushy state and the final solid state taking place in a solidification process. In particular, a pure viscous deviatoric model has been used at the initial fluid-like state. A thermomecanical contact model, including a frictional hardening and temperature dependent coupled potential, is derived within a fully consistent thermodinamical theory. The numerical model has been implemented into the computational finite element code COMET developed by the authors. numerical s
暂无评论