Engineering topics
Publications and source records attributed to Lakin, William D..
A system and method of modeling a circulatory system including a regulatory mechanism parameter. In one embodiment, a regulatory mechanism parameter in a lumped parameter model is represented as a logistic function. In another embodiment, the circulatory system model includes a compliant vessel, the model having a parameter representing a change in pressure due to contraction of smooth muscles of a wall of the vessel.
A whole-body mathematical model (10) for simulating intracranial pressure dynamics. In one embodiment, model (10) includes 17 interacting compartments, of which nine lie entirely outside of intracranial vault (14). Compartments (F) and (T) are defined to distinguish ventricular from extraventricular CSF. The vasculature of the intracranial system within cranial vault (14) is also subdivided into five compartments (A, C, P, V, and S, respectively) representing the intracranial arteries, capillaries, choroid plexus, veins, and venous sinus. The body's extracranial systemic vasculature is divided into six compartments (I, J, O, Z, D, and X, respectively) representing the arteries, capillaries, and veins of the central body and the lower body. Compartments (G) and (B) include tissue and the associated interstitial fluid in the intracranial and lower regions. Compartment (Y) is a composite involving the tissues, organs, and pulmonary circulation of the central body and compartment (M) represents the external environment.
The dynamic behavior of flexible rotating beams continues to receive considerable research attention as it constitutes a fundamental problem in applied mechanics. Further, beams comprise parts of many rotating structures of engineering significance. A topic of particular interest at the present time involves the development of techniques for obtaining the behavior in both space and time of a rotor acted upon by a simple airload loading. Most current work on problems of this type use solution techniques based on normal modes. It is certainly true that normal modes cannot be disregarded, as knowledge of natural blade frequencies is always important. However, the present work has considered a computational structural mechanics (CSM) approach to rotor blade dynamics problems in which the physical properties of the rotor blade provide input for a direct numerical solution of the relevant boundary-and-initial-value problem. Analysis of the dynamics of a given rotor system may require solution of the governing equations over a long time interval corresponding to many revolutions of the loaded flexible blade. For this reason, most of the common techniques in computational mechanics, which treat the space-time behavior concurrently, cannot be applied to the rotor dynamics problem without a large expenditure of computational resources. By contrast, the integrating matrix technique of computational mechanics has the ability to consistently incorporate boundary conditions and 'remove' dependence on a space variable. For problems involving both space and time, this feature of the integrating matrix approach thus can generate a 'splitting' which forms the basis of an efficient CSM method for numerical solution of rotor dynamics problems.
Work completed under the current grant comprises the start of a theoretical and computational attack on the subharmonic route to secondary instabilities in compressible flows. The total flow field in this problem is made up of the following components: (1) a steady streamwise mean boundary layer flow which depends only on the normal space component y; (2) a two-dimensional time dependent T-S wave which moves with wavespeed c and has no spanwise dependence; and (3) a fully three-dimensional, time dependent T-S wave whose streamwise wavenumber is half of the streamwise wavenumber associated with the two-dimensional T-S wave in b. If a frame of reference is adopted which moves with the wavespeed c of the 2-D T-S wave, the time dependence of this portion of the flow can be eliminated. The effective steady mean flow in this problem is now the sum of the original parallel steady mean flow and the initial 2-D T-S instability. Dependence on the streamwise coordinate x in this mean flow can be extracted by assuming normal mode expansions involving complex exponentials and the streamwise wavenumber a. However, it is important to note that, because this is a wave-wave interaction problem, unlike the usual linear instability case, both the complex exponential, and its complex conjugate, must be retained in describing the 2-D T-S wave. The role of the perturbation to the steady mean flow is now played by the 3-D time dependent T-S wave. In treating this wave, normal modes in the streamwise and spanwise directions and time may be used. Consistent with the subharmonic nature of this transition route, the streamwise wavenumber is a/2, and complex conjugates of the complex exponential must be employed. This is not the case with the modes giving z and t dependence with wavespeed o and spanwise wavenumber B as the effective mean flow quantities are independent of z and their time dependence is accounted for by the moving frame of reference. Consequently, the wave-wave interaction which will produce mean flow modification occurs only through the streamwise exponentials.
Originally analytical and numerical models were to be developed for noise production in supersonic jets, wakes and free shear layers. While the effort was concentrated initially on these aspects, other topics were also pursued, most were of interest to the Jet Noise Group of the Aeroacoustics Branch. An overview is given of subjects reviewed and the investigations that were carried out. A significant effort was devoted to numerically predicting the flow field of a turbulent supersonic wall jet. This information is necessary for computing the pressure in the far field. The wall jet was selected because it represents a generic flow that can be associated with plug nozzle in supersonic engines. It combines the characteristic of a boundary layer with that of a free shear flow. The spatially evolving flow obtained using Dash's code would form the input for the stability analysis program. This analysis would determine the large scale instability wave within the flow. The far field pressure can be computed from the shape of the evolving large scale structure by asymptotic methods. Flow characteristics obtained from a program that analyses the turbulent downstream supersonic flow in a nozzle are described and compared with experimental results.
The integrating matrix technique of computational mechanics is extended to include the effects of concentrated masses. The stability of a flexible rotating beam with discrete masses is analyzed to determine the critical rotational speeds for buckling in the inplane and out-of-plane directions. In this problem, the beam is subjected to compressive centrifugal forces arising from steady rotation about an axis which does not pass through the clamped end of the beam. To determine the eigenvalues from which stability is assessed, the differential equations of motion are solved numerically by combining the extended integrating matrix method with an eigenanalysis. Stability boundaries for a discrete mass representation of a uniform beam are shown to asymptotically approach the stability boundaries for the corresponding continuous mass beam as the number of concentrated masses is increased. An error in the literature is also noted for the discrete mass problem concerning the behavior of the critical rotational speed for inplane buckling as the radius of rotation of the clamped end of the beam is reduced.
The fully nonlinear development of small wavelength Goertler vortices in a growing boundary layer is investigated using a combination of asymptotic and numerical methods. The starting point for the analysis is the weakly nonlinear theory of Hall (1982b) who discussed the initial development of small amplitude vortices in a neighborhood of the location where they first become linearly unstable. That development is unusual in the context of nonlinear stability theory in that it is not described by the Stuart-Watson approach. In fact the development is governed by a pair of coupled nonlinear partial differential evolution equations for the vortex flow and the mean flow correction. Here the further development of this interaction is considered for vortices so large that the mean flow correction driven by them is as large as the basic state. Surprisingly it is found that such a nonlinear interaction can still be described by asymptotic means. It is shown that the vortices spread out across the boundary layer and effectively drive the boundary layer. In fact the system obtained by writing down the equations for the fundamental component of the vortex generate a differential equation for the basic state. Thus the mean flow adjusts so as to make these large amplitude vortices locally neutral. Moreover in the region where the vortices exist the mean flow has a square-root profile and the vortex velocity field can be written down in closed form. The upper and lower boundaries of the region of vortex activity are determined by a free-boundary problem involving the boundary layer equations. In general it is found that this region ultimately includes almost all of the original boundary layer and much of the free-stream. In this situation the mean flow has essentially no relationship to the flow which exists in the absence of the vortices.