Vibration of semi-infinite cylindrical shells.
Semiinfinite Timoshenko-Love shell vibration, taking double Laplace transform of both Timoshenko and Vlasov equations of motion with respect to time and axial space variables
Berger, B. S.↗
SEARCH · Engineering Papers
Results for “LAPLACE EQUATION”
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.
Unsteady Aerodynamics - Subsonic Compressible Inviscid Case
This paper presents a new analytical treatment of Unsteady Aerodynamics - the linear theory covering the subsonic compressible (inviscid) case - drawing on some recent work in Operator Theory and Functional Analysis. The specific new results are: (a) An existence and uniqueness proof for the Laplace transform version of the Possio integral equation as well as a new closed form solution approximation thereof. (b) A new representation for the time-domain solution of the subsonic compressible aerodynamic equations emphasizing in particular the role of the initial conditions.
Balakrishnan, A. V.↗
Continuum Modeling of the Dynamics of Externally Injection-locked Coupled Oscillator Arrays
Mutually injection-locked arrays of electronic oscillators provide a novel means of controlling the aperture phase of a phased-array antenna, thus achieving the advantages of spatial power combining while retaining the ability to steer the radiated beam. In a number of design concepts, one or more of the oscillators are injection locked to a signal from an external master-oscillator. The behavior of such a system has been analyzed by numerical solution of a system of nonlinear differential equations which, due to its complexity, yields limited insight into the relationship between the injection signals and the aperture phase. In this paper, we develop a continuum model, which results in a single partial differential equation for the aperture phase as a function of time. Solution of the equation is effected by means of the Laplace transformation and yields detailed information concerning the dynamics of the array under the influence of the external injection signals.
injection locked↗
A 6th Order Mehrstellen Finite Volume Discretization of Poisson's Equation in Three Dimensions
We discuss the derivation of a new, sixth-order finite volume scheme for Poisson’s equation on 3D Cartesian equispaced grids. The scheme is based on a discretization of the Laplace operator with a compact (Mehrstellen) 27-point stencil. To achieve sixth order convergence the right hand side of the equation is replaced with a discrete operator that involves the discrete Laplace and Biharmonic operators and the sum of discrete fourth-order cross derivatives applied to the charge function. Numerical tests demonstrate the superiority of the proposed method compared to the well known schemes associated with the 7-point and 19-point discretizations of the Laplacian.
97 MATHEMATICS AND COMPUTING↗
Bayesian Symbolic Regression: Addressing Challenges in Estimating Fractional Bayes Factors and Application to Fatigue Crack Growth Modeling
This research pioneers advancements in computational mechanics by integrating Bayesian-based uncertainty quantification into symbolic regression, specifically focusing on the critical task of accurately estimating the fractional Bayes factor for selecting arbitrary equations. In our exploration, we rigorously study two prominent methods—sequential Monte Carlo and the Laplace approximation—employed for computing the fractional Bayes factor. Our findings underscore the limitations of the Laplace approximation, revealing its diminished accuracy in nonlinear and multimodal scenarios. Specifically, the Laplace approximation is shown to underpredict fractional Bayes factor on a wide set of equations associated with a symbolic regression benchmark. This comparative analysis sheds light on the nuanced performance of these techniques, guiding researchers toward more informed choices in uncertainty quantification within symbolic regression. Furthermore, we showcase the practical utility of these enhanced symbolic regression tools through their application to a real-world problem in fatigue crack growth modeling, emphasizing their efficacy in capturing the complexities of mechanical systems.
Geoffrey Bomarito↗
Automated reaction kinetics and network exploration (Arkane): A statistical mechanics, thermodynamics, transition state theory, and master equation software
The open-source statistical mechanics software described here, Arkane–Automated Reaction Kinetics and Network Exploration–facilitates computations of thermodynamic properties of chemical species, high-pressure limit reaction rate coefficients, and pressure-dependent rate coefficient over multi-well molecular potential energy surfaces (PES) including the effects of collisional energy transfer on phenomenological kinetics. Arkane can use estimates to fill in information for molecules or reactions where quantum chemistry information is missing. The software solves the internal energy master equation for complex unimolecular reaction systems. Inputs to the software include converged electronic structure computations performed by the user using a variety of supported software packages (Gaussian, Molpro, Orca, TeraChem, Q-Chem, Psi4). The software outputs high-pressure limit rate coefficients and pressure-dependent phenomenological rate coefficients, as well as computed thermodynamic properties (enthalpy, entropy, and constant pressure heat capacity) with added energy corrections. Some of the key features of Arkane include treatment of 1D, 2D or ND hindered internal rotation modes, treatment of free internal rotation modes, quantum tunneling effect consideration, transition state theory (TST) and Rice-Ramsperger-Kassel-Marcus (RRKM) rate coefficient computations, master equation solution with four implemented methods, inverse-Laplace transform of high-pressure limit rate coefficients into the energy domain, energy corrections based on bond-additivity or isodesmic reactions, automated and efficient PES exploration, and PES sensitivity analysis. The present work describes the design of Arkane, how it should be used, and refers to the theory that it employs. Arkane is distributed via the RMG-Py software suite (https://github.com/ReactionMechanismGenerator/RMG-Py).
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Lossy radial diffusion of relativistic Jovian electrons
The radial diffusion equation with synchrotron losses was solved by the Laplace transform method for near-equatorially mirroring relativistic electrons. The evolution of a power law distribution function was found and the characteristics of synchrotron burn-off are stated in terms of explicit parameters for an arbitrary diffusion coefficient. Emissivity from the radiation belts of Jupiter was studied. Asymptotic forms for the distribution in the strong synchrotron loss regime are provided.
Barbosa, D. D.↗
Modelling a single phase voltage controlled rectifier using Laplace transforms
The development of a 20 kHz, AC power system by NASA for large space projects has spurred a need to develop models for the equipment which will be used on these single phase systems. To date, models for the AC source (i.e., inverters) have been developed. It is the intent of this paper to develop a method to model the single phase voltage controlled rectifiers which will be attached to the AC power grid as an interface for connected loads. A modified version of EPRI's HARMFLO program is used as the shell for these models. The results obtained from the model developed in this paper are quite adequate for the analysis of problems such as voltage resonance. The unique technique presented in this paper uses the Laplace transforms to determine the harmonic content of the load current of the rectifier rather than a curve fitting technique. Laplace transforms yield the coefficient of the differential equations which model the line current to the rectifier directly.
Kraft, L. Alan↗
Kinetic theory of particle-in-cell simulation plasma and the ensemble averaging technique
Abstract We derive the kinetic theory of fluctuations in physically and numerically stable particle-in-cell (PIC) simulations of electrostatic plasmas. The starting point is the single-time correlation at the start of the simulation between the statistical fluctuations of the weighted densities of macroparticle centers in the plasma particle phase-space. The fluctuations are associated with different initial conditions, typically due to the random initial conditions (in velocity space) of the macroparticles/simulation plasma, assigned according to their initial distribution of probability. The single-time correlations at all time steps and in each spatial grid cell are then determined from the Laplace–Fourier transforms of the discretized Klimontovich-like equation for the macroparticles and Maxwell’s equations for the fields, as computed by modern PIC codes. We recover the expressions for the electrostatic field and the plasma particle density fluctuation autocorrelation spectra as well as the kinetic equations describing the average evolution of PIC-simulated plasma particles, first derived by Langdon (1970b Proc. 4th Conf. Numerical Simulation of Plasmas ) using a test macroparticle approach perturbing a discretized Vlasovian plasma and then averaging the obtained physical quantity over the initial macroparticle velocity distribution. We generalize and extend these results to the modern algorithms in PIC codes using arbitrary macroparticle weights. Analytical estimates of statistical fluctuation amplitudes are derived as a function of the plasma simulation parameters, using the central limit theorem in the limit of a large number of macroparticles per cell. The theory is then used to analyze the ensemble averaging technique of PIC simulations where statistical averages are performed over ensembles of PIC simulations, modeling the same plasma physics problem but using different statistical realizations of the initial distribution functions of the macroparticles. This method is illustrated by linear Landau damping uncovering (from noise, which is usually considered numerical) the physical fluctuations driven by a single small amplitude electrostatic wave perturbing a PIC simulation plasma in equilibrium.
fluctuations correlations↗
Unsteady transonic flow analysis for low aspect ratio, pointed wings.
Oswatitsch and Keune's parabolic method for steady transonic flow is applied and extended to thin slender wings oscillating in the sonic flow field. The parabolic constant for the wing was determined from the equivalent body of revolution. Laplace transform methods were used to derive the asymptotic equations for pressure coefficient, and the Adams-Sears iterative procedure was employed to solve the equations. A computer program was developed to find the pressure distributions, generalized force coefficients, and stability derivatives for delta, convex, and concave wing planforms.
Kimble, K. R.↗
Lossy radial diffusion of relativistic Jovian electrons
The radial diffusion equation with synchrotron losses is solved by the Laplace-transform method for near equatorially mirroring relativistic electrons. The evolution of a power-law distribution function is found, and the characteristics of synchrotron burnoff are stated in terms of explicit parameters for an arbitrary diffusion coefficient of a specific form. The peaking of the 10.4-cm volume emissivity from Jupiter at an L shell of about 1.8 provides an estimate of the diffusion coefficient in the radiation belts; one value is suggested as the appropriate modification, for an equatorial field strength of 4.2 G, of the Birmingham et al. (1974) result. Nonsynchrotron losses are included phenomenologically; from the phase-space densities reported by McIlwain and Fillius (1975), the particle lifetime is estimated. Asymptotic forms for the distribution in the strong synchrotron loss regime are provided.
Barbosa, D. D.↗
Steady, Oscillatory and Unsteady, Subsonic and Supersonic Aerodynamics (SOUSSA) for complex aircraft configurations
The Green's function method and the computer program SOUSSA (Steady Oscillatory and Unsteady Subsonic and Supersonic Aerodynamics) are reviewed. The Green's function method is applied to the fully unsteady potential equation yielding an integro-differential-delay equation. This equation is approximated by a set of differential-delay equations in time using the finite element method. The Laplace transform is used to yield a matrix relating the velocity potential to the normal wash. The matrix of the generalized aerodynamic forces is obtained by premultiplying and postmultiplying the matrices relating generalized forces to the potential and the normal wash by the generalized coordinates. The program SOUSSA is compared with existing numerical results. Results indicate that the program is not only general, flexible, and easy to use, but also accurate and fast.
Morino, L.↗
Time-temperature effect in adhesively bonded joints
The viscoelastic analysis of an adhesively bonded lap joint was reconsidered. The adherends are approximated by essentially Reissner plates and the adhesive is linearly viscoelastic. The hereditary integrals are used to model the adhesive. A linear integral differential equations system for the shear and the tensile stress in the adhesive is applied. The equations have constant coefficients and are solved by using Laplace transforms. It is shown that if the temperature variation in time can be approximated by a piecewise constant function, then the method of Laplace transforms can be used to solve the problem. A numerical example is given for a single lap joint under various loading conditions.
Delale, F.↗
An Overview of the Patch Integral Method (PIM), a New Heat Transfer Analysis Tool for Hypersonic Wind Tunnel Facilities at NASA Langley
NASA Langley’s hypersonic wind tunnels are heavily leveraged for planetary missions. The data collection method in these tunnels is thermography, and surface temperature measurements of the model surface are collected and reduced to produce surface heating data, as seen in Fig 1. However, during model injection, no temperature data are collected, and thus conventional, integral heat transfer methods cannot be used to solve for surface heating. A method was developed in the 1990’s to reduce this data despite the data gap, known as the step approximation method. The method assumes that the film coefficient behaves as a step function, the model is semi-infinite, and thermal properties are constant. With these simplifying assumptions, a Laplace transform can be performed to result in an equation that takes the initial temperature of the model and a temperature at some point in time to back out the film coefficient at that time. This is the method that is used in the current thermographic data reduction software, IHEAT. While computationally light-weight, the step approximation has several issues associated with it. The time-history of temperature is not accounted for, which is vital as heat transfer is an integral process. Additionally, the required semi-infinite assumption is unnecessary and might be violated during runtime. Thermal variation of material properties can have a sizeable impact on heating results and are not modeled by the method. This method also takes multiple seconds to “collapse” to a steady state, which is undesirable from both a facility and data reduction standpoint. The method is also very sensitive to the “effective time” approximation, an approximation of when heating instantaneously starts (which is a nonphysical simplification), and a small variation in this value can result in an error in heating results.
J. S. Cheatwood↗
Fourth order difference methods for hyperbolic IBVP's
Fourth order difference approximations of initial-boundary value problems for hyperbolic partial differential equations are considered. We use the method of lines approach with both explicit and compact implicit difference operators in space. The explicit operator satisfies an energy estimate leading to strict stability. For the implicit operator we develop boundary conditions and give a complete proof of strong stability using the Laplace transform technique. We also present numerical experiments for the linear advection equation and Burgers' equation with discontinuities in the solution or in its derivative. The first equation is used for modeling contact discontinuities in fluid dynamics, the second one for modeling shocks and rarefaction waves. The time discretization is done with a third order Runge-Kutta TVD method. For solutions with discontinuities in the solution itself we add a filter based on second order viscosity. In case of the non-linear Burger's equation we use a flux splitting technique that results in an energy estimate for certain different approximations, in which case also an entropy condition is fulfilled. In particular we shall demonstrate that the unsplit conservative form produces a non-physical shock instead of the physically correct rarefaction wave. In the numerical experiments we compare our fourth order methods with a standard second order one and with a third order TVD-method. The results show that the fourth order methods are the only ones that give good results for all the considered test problems.
Gustafsson, Bertil↗
Finite transform solution of the vibrating annular membrane.
Vibrating annular membrane problem, including load per unit area and asymmetry of load and vibration, using finite transform derivation
LOAD DISTRIBUTION↗
One-dimensional Coulomb-damped wave motion in prismatic bars
Study analyzes wave motions in prismatic bars with Coulomb damping, using Laplace transforms as an aid in solving the partial differential equations. The results are detailed.
Tomlin, D. D., Jr.↗
Start-On-The-Part Transient Model for In-Situ Automated Tape Placement of Thermoplastic Composites
Fabrication of a complex part by automated tape placement (ATP) can require starting up a new tape-end in the part interior, termed start-on-the-part. Careful thermal management of the starting transient is needed to achieve uniform crystallinity and inter-laminar weld strength - which is the objective of this modeling effort. The transient is modeled by a Fourier-Laplace transform solution of the time-dependent thermal transport equation in two spatial dimensions. The solution is subject to a quasi-steady approximation for the speed and length of the consolidation head. Sample calculations are done for the Langley ATP robot applying PEEK/carbon fiber composite and for two upgrades in robot performance. The head starts out almost at rest which meets an engineering requirement for accurate placement of the new tape-end. The head then rapidly accelerates until it reaches its steady state speed. This rapid acceleration, however, violates the quasi-steady approximation, so uniform weld strength and crystallinity during the starting transient are not actually achieved. The solution does give the elapsed time and distance from start-up to validity of the quasi-steady approximation - which quantifies the length of the non-uniform region. The elapsed time was always less than 0.1 s and the elapsed distance less than 1 cm. This quantification would allow the non-uniform region to be either trimmed away or compensated for in the design of a part. Such compensation would require experiments to measure the degree of non-uniformity, because the solution does not provide this information. The rapid acceleration suggests that the consolidation roller or belt be actively synchronized to avoid abrading the tape.
Costen, Robert c.↗