The propagating fracture of constant shape.
Propagating tensile fracture of constant ellipse shape treated as steady state dynamic mixed boundary value problem
SEARCH · Search NASA
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.
Propagating tensile fracture of constant ellipse shape treated as steady state dynamic mixed boundary value problem
This paper briefly reviews convergent finite difference schemes for hyperbolic initial boundary value problems and their applications to boundary control systems of hyperbolic type which arise in the modelling of vibrations. These difference schemes are combined with the primal and the dual approaches to compute the optimal control in the unconstrained case, as well as the case when the control is subject to inequality constraints. Some of the preliminary numerical results are also presented.
A "simplified shear solution" method approximates the through-thickness shear stress distribution within a composite laminate based on an extension of laminated beam theory. The method does not consider the solution of a particular boundary value problem; rather, it requires only knowledge of the global shear loading, geometry, and material properties of the laminate or panel. It is thus analogous to lamination theory in that ply-level stresses can be efficiently determined from global load resultants at a given location in a structure and used to evaluate the margin of safety on a ply-by-ply basis. The simplified shear solution stress distribution is zero at free surfaces, continuous at ply boundaries, and integrates to the applied shear load. The method has been incorporated within the HyperSizer commercial structural sizing software to improve its predictive capability for designing composite structures. The HyperSizer structural sizing software is used extensively by NASA to design composite structures. In the case of through-thickness shear loading on panels, HyperSizer previously included a basic, industry-standard, method for approximating the resulting shear stress distribution in sandwich panels. However, no such method was employed for solid laminate panels. The purpose of the innovation is to provide an approximation of the through-thickness shear stresses in a solid laminate given the through-thickness shear loads (Qx and Qy) on the panel. The method was needed for implementation within the HyperSizer structural sizing software so that the approximated ply-level shear stresses could be utilized in a failure theory to assess the adequacy of a panel design. The simplified shear solution method was developed based on extending and generalizing bi-material beam theory to plate-like structures. It is assumed that the through-thickness shear stresses arise due to local bending of the laminate induced by the through-thickness shear load, and by imposing equilibrium both vertically and horizontally, the through-thickness shear stress distribution can be calculated. The resulting shear stresses integrate to the applied shear load, are continuous at the ply interfaces, and are zero at the laminate-free surfaces. If both Qx and Qy shear loads are present, it is assumed that they act independently and that their effects can be superposed. The calculated shear stresses can be rotated within each ply to the principal material coordinates for use in a ply-level failure criterion. The novelty of the simplified shear solution method is its simplicity and the fact that it does not require solution of a particular boundary value problem. The advantages of the innovation are that an approximation of the though-thickness shear stress distribution can be quickly determined for any solid laminate or solid laminate region within a stiffened panel.
A class of solutions to a model of forced oscillations in a rotating stratified atmospheric layer is derived and analyzed. The basic model is found to reduce to a boundary value problem with a second-order linear partial differential equation of the hyperbolic type for this range of forcing frequencies. The forced solutions are shown to exhibit resonances with the normal modes of oscillation of the layer. The characteristics of the resonant modes are analyzed in terms of mean tropospheric values of temperature, temperature lapse, wind speed, horizontal and vertical wind shears, latitude, and the frequency and horizontal wavelength of the forcing mechanism. These solutions are compared with solutions to the model for a different (subinertial) range of forcing frequencies. This comparison leads to an elliptic boundary value problem. The solutions in that case do not exhibit the same type of resonance and generally decay away from the region of forcing.
The stability of viscous, incompressible flow along a streamwise corner, often called the corner boundary layer problem is examined. The semi-infinite boundary value problem satisfied by small amplitude disturbances in the "bending boundary layer' region is obtained. The mean secondary flow induced by the corner exhibits a flow reversal in this region. Uniformly valid "first approximations' to solutions of the governing differential equations are derived. Uniformity at infinity is achieved by a suitable choice of the large parameter and use of an approximate Langer variable. Approximations to solutions of balanced type have a phase shift across the critical layer which is associated with instabilities in the case of two dimensional boundary layer profiles.
The stability of viscous, incompressible flow along a streamwise corner, often called the corner boundary layer problem is examined. The semi-infinite boundary value problem satisfied by small amplitude disturbances in the 'bending boundary layer' region is obtained. The mean secondary flow induced by the corner exhibits a flow reversal in this region. Uniformly valid 'first approximations' to solutions of the governing differential equations are derived. Uniformity at infinity is achieved by a suitable choice of the large parameter and use of an approximate Langer variable. Approximations to solutions of balanced type have a phase shift across the critical layer which is associated with instabilities in the case of two dimensional boundary layer profiles. Previously announced in STAR as N84-17532
Equivalence of perturbation method, Newtons method, and adjoint-perturbation method for solving boundary value problems from variational calculus treatment of optimal control problems
Oscillating wing aerodynamic load boundary value problem reduction to sequence of steady lifting- surface problems
Explore the source record for details and available documents.
Finite-difference Navier-Stokes calculations for unsteady, three-dimensional, incompressible, viscous flows induced by initial vorticity distributions are presented and discussed in this paper. The initial vorticity distributions are assumed to be embedded in a flow field of infinite extent that is quiescent at infinity. These vorticity distributions are typical of vortex rings and other closed vortical tubes or structures. Such structures are important elements in fluid flows such as jets, atmospheric convection and the far-field wakes of aircraft; studies of their interaction may aid in an understanding of complex fluid flows. The calculations employ a method recently proposed by Ting to approximate the infinite-domain boundary value problem with a finite boundary computational domain, and this method is shown to yield accurate three-dimensional results for reasonable expenditures of computer time. Because of the efficiency of the boundary condition technique and the resulting Navier-Stokes code, a 16-bit minicomputer with virtual memory was capable of performing the calculations for the unsteady motion of two obliquely colliding vortex rings. The results of these calculations are presented in the paper.
Underintegrated methods are investigated with respect to their stability and convergence properties. The focus was on identifying regions where they work and regions where techniques such as hourglass viscosity and hourglass control can be used. Results obtained show that underintegrated methods typically lead to finite element stiffness with spurious modes in the solution. However, problems exist (scalar elliptic boundary value problems) where underintegrated with hourglass control yield convergent solutions. Also, stress averaging in underintegrated stiffness calculations does not necessarily lead to stable or convergent stress states.
The paper presents an application of adaptive finite element methods to the modeling of low-cycle continuum damage and life prediction of high-temperature components. The major objective is to provide automated and accurate modeling of damaged zones through adaptive mesh refinement and adaptive time-stepping methods. The damage modeling methodology is implemented in an usual way by embedding damage evolution in the transient nonlinear solution of elasto-viscoplastic deformation problems. This nonlinear boundary-value problem is discretized by adaptive finite element methods. The automated h-adaptive mesh refinements are driven by error indicators, based on selected principal variables in the problem (stresses, non-elastic strains, damage, etc.). In the time domain, adaptive time-stepping is used, combined with a predictor-corrector time marching algorithm. The time selection is controlled by required time accuracy. In order to take into account strong temperature dependency of material parameters, the nonlinear structural solution a coupled with thermal analyses (one-way coupling). Several test examples illustrate the importance and benefits of adaptive mesh refinements in accurate prediction of damage levels and failure time.
Thermodynamics and higher-order fluid theories - boundary value problems for unsteady simple shear flow
A new inverse method for the aerodynamic design of airfoils is presented for subcritical flows. The pressure distribution in this method can be prescribed as a function of the arclength of the still unknown body. It is shown that this inverse problem is mathematically equivalent to solving only one nonlinear boundary value problem subject to known Dirichlet data on the boundary.
A numerical algorithm is presented for solving laminar, steady, supersonic interacting boundary-layer flows for quasi-three-dimensional configurations. The interaction problem is treated as a boundary-value problem and a salient feature of the scheme is the direct implementation of the downstream boundary condition. Solutions are presented for axisymmetric and swept (yawed) compression ramps for both adiabatic and heat transfer conditions over a Mach number range of 2-6. The results are in good agreement with experimental data and existing theories for axisymmetric cases. For the swept (yawed) configurations, lack of experimental data makes a direct comparison impossible, but the present solutions are found to be in qualitative agreement with earlier studies. In addition it is shown that the trends obtained here for the sweep effects are well predicted by a simple extension of the two-dimensional asymptotic theory.
Quasi-optimum control technique based on Hamiltonian to solve two-point boundary value problem for nonlinear system
Rarefied gas dynamics - solution of boundary value problem of boltzmann equation at arbitrary knudsen number by transfer equations
Matrix algebra solution, using coordinate systems transformations, to boundary value problem for structures composed of shell elements