Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Moment equations”

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.

At least 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Toward UV models of kinetic mixing and portal matter. V. Indirect probes of the new physics scale

Kinetic mixing of the dark photon, the gauge boson of a hidden U ( 1 ) D , with the Standard Model (SM) gauge fields to induce an interaction between ordinary matter and dark matter (DM) at 1-loop requires the existence of portal matter (PM) fields having both dark and SM charges. As discussed in earlier work, these same PM fields can also lead to other loop-level mechanisms besides kinetic mixing that can generate significant interactions between SM fermions and the dark photon in a manner analogous to those that can be generated between a Dirac neutrino and a SM photon, i.e., dark moments. In either case, there are reasons to believe, e.g., due to the renormalization group equation running of the U ( 1 ) D gauge coupling, that PM fields may have ∼ TeV -scale masses that lie at or above those directly accessible to the HL-LHC. If they lie above the reach of the HL-LHC, then the only way to possibly explore the physics at this high scale in the short term is via indirect measurements made at lower energies, e.g., at lepton colliders operating in the m Z to 1 TeV range. In particular, processes such as e + e − → γ + DM or e + e − → f ¯ f , where f is a SM fermion, may be most useful in this regard. Here we explore these possibilities within the framework of a simple toy PM model, introduced in earlier work, based on a non-Abelian dark gauge group completion operating at the PM scale. In the kinetic mixing setup, we show these efforts fail due to the inherently tiny cross sections in the face of substantial SM backgrounds. However, in the case of interactions via induced dark moments, since they necessarily take the form of higher dimensional operators whose influence grows with energy, we show that access to PM-scale information may become possible for certain ranges of the toy model parameters for both of these e + e − processes at a 1 TeV collider. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

A hybrid Monte Carlo-deterministic second moment method with efficient variance reduction

In this work, we present a hybrid method that combines Monte Carlo with deterministic finite element methods to solve a linear Boltzmann transport equation. Our hybrid method runs orders of magnitude faster than Monte Carlo, without sacrificing accuracy, for a proxy problem from radiative transfer that contains both optically-thick and optically-thin material. We believe that this is the first demonstration of a hybrid Second Moment Method in more than one spatial dimension, the first to consider more than one material, and the first to use variance reduction. Our variance reduction approach arises from an asymptotic analysis in which we show that the magnitude of the scattering source grows without bound. We transform the problem to compute the deviation of the radiation intensity from isotropy. The magnitude of the source in the transformed problem is bounded, and the quality of the hybrid method solution is dramatically improved by a substantial reduction in the variance.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Coupling turbulent flow with blade aeroelastics and control modules in large-eddy simulation of utility-scale wind turbines

We present a large-eddy simulation framework capable of control co-design of large wind turbines, coupling the turbulent flow environment with blade aeroelastics and turbine controllers. The geometry and aerodynamics of the rotor blades and the turbine nacelle are parameterized using an actuator surface model. The baseline collective pitch control and individual pitch control (IPC) algorithms, consisting of a single-input, single-output proportional–integral controller and two integral controllers, respectively, are incorporated into the simulation framework. Furthermore, a second-order model based on the Euler–Bernoulli beam theory is implemented to describe the blade deformation. Simulations are carried out to investigate the impact of collective and individual pitch control strategies on the deflection of turbine blades. Our results show that the IPC reduces the blade tip deflection fluctuations in the out-of-plane direction, while the fluctuations of the blade tip deflection along the in-plane direction are barely affected by the IPC. Furthermore, the blade out-of-plane deformation fluctuation is underestimated by the one-way coupling approach compared to the two-way coupling approach. Furthermore, the findings of this study reveal the importance of advanced control systems in reducing the dynamic loads on wind turbine blades and underscore the potential of control co-design to reduce the levelized cost of wind energy.

42 ENGINEERING↗
DOE OSTI2023

ANN-based ground motion model for Turkey using stochastic simulation of earthquakes

SUMMARY Turkey is characterized by a high level of seismic activity attributed to its complex tectonic structure. The country has a dense network to record earthquake ground motions; however, to study previous earthquakes and to account for potential future ones, ground motion simulations are required. Ground motion simulation techniques offer an alternative means of generating region-specific time-series data for locations with limited seismic networks or regions with seismic data gaps, facilitating the study of potential catastrophic earthquakes. In this research, a local ground motion model (GMM) for Turkey is developed using region-specific simulated records, thus constructing a homogeneous data set. The simulations employ the stochastic finite-fault approach and utilize validated input-model parameters in distinct regions, namely Afyon, Erzincan, Duzce, Istanbul and Van. To overcome the limitations of linear regression-based models, artificial neural network is used to establish the form of equations and coefficients. The predictive input parameters encompass fault mechanism (FM), focal depth (FD), moment magnitude (Mw), Joyner and Boore distance (RJB) and average shear wave velocity in the top 30 m (Vs30). The data set comprises 7359 records with Mw ranging between 5.0 and 7.5 and RJB ranging from 0 to 272 km. The results are presented in terms of spectral ordinates within the period range of 0.03–2.0 s, as well as peak ground acceleration and peak ground velocity. The quantification of the GMM uncertainty is achieved through the analysis of residuals, enabling insights into inter- and intra-event uncertainties. The simulation results and the effectiveness of the model are verified by comparing the predicted values of ground motion parameters with the observed values recorded during previous events in the region. The results demonstrate the efficacy of the proposed model in simulating physical phenomena.

Karimzadeh, Shaghayegh (ORCID:0000000337531676)↗
NASA NTRS1985

Qualitative properties of large buckled states of spherical shells

A system of 6th-order quasi-linear Ordinary Differential Equations is analyzed to show the global existence of axisymmetrically buckled states. A surprising nodal property is obtained which shows that everywhere along a branch of solutions that bifurcates from a simple eigenvalue of the linearized equation, the number of simultaneously vanishing points of both shear resultant and circumferential bending moment resultant remains invariant, provided that a certain auxiliary condition is satisfied.

Shih, K. G.↗
NASA NTRS1989

Differential Curing In Fiber/Resin Laminates

Modified layup schedule counteracts tendency toward delamination. Improved manufacturing process resembles conventional process, except prepregs partially cured laid on mold in sequence in degree of partial cure decreases from mold side to bag side. Degree of partial cure of each layer at time of layup selected by controlling storage and partial-curing temperatures of prepreg according to Arrhenius equation for rate of gel of resin as function of temperature and time from moment of mixing. Differential advancement of cure in layers made large enough to offset effect of advance bag-side heating in oven or autoclave. Technique helps prevent entrapment of volatile materials during manufacturing of fiber/resin laminates.

Webster, Charles N.↗
NASA NTRS1990

An EPR methodology for measuring the London penetration depth for the ceramic superconductors

The use is discussed of electron paramagnetic resonance (EPR) as a quick and easily accessible method for measuring the London penetration depth, lambda for the high T(sub c) superconductors. The method utilizes the broadening of the EPR signal, due to the emergence of the magnetic flux lattice, of a free radical adsorbed on the surface of the sample. The second moment, of the EPR signal below T(sub c) is fitted to the Brandt equation for a simple triangular lattice. The precision of this method compares quite favorably with those of the more standard methods such as micro sup(+)SR, Neutron scattering, and magnetic susceptibility.

Rakvin, B.↗
NASA NTRS1991

Microwave (EPR) measurements of the penetration depth measurements of high-Tc superconductors

The use is discussed of electron paramagnetic resonance (EPR) as a quick and easily accessible method for measuring the London penetration depth, lambda for the high T sub c superconductors. The method uses the broadening of the EPR signal, due to the emergence of the magnetic flux lattice, of a free radical adsorbed on the surface of the sample. The second moment, of the EPR signal below T sub c is fitted to the Brandt equation for a simple triangular lattice. The precision of this method compares quite favorably with those of the more standard methods such as micro sup(+)SR, neutron scattering, and magnetic susceptibility.

Dalal, N. S.↗
NASA NTRS1948

Damping in Pitch and Roll of Triangular Wings at Supersonic Speeds

The expressions for the damping derivatives in pitch and roll of triangular wings are derived by means of the linearized theory. In the method used, the wing is represented by an unknown distribution of doublets. An integral equation containing the unknown distribution is set up and solved by analogy with known incompressible flow relations. It is pointed out that the results may be used to obtain damping coefficients of a limited series of sweptback wings, the most interesting of which are the so-called "arrow wings."

EQUATIONS - MOTION↗
NASA NTRS1931

A study of curvilinear flight

When an airplane describes a curve it takes a certain time for it to turn from level to inclined position and then back to level again. In the following, we express the motion about the horizontal axis as "roll" or "bank" and the motion perpendicular to the vertical axis, i.e., the actual curve, as "turn." Equations and tables provide results on various aspects of turns, control settings, acceleration, inertia moments, and angular velocity.

Kruse, Helmuth↗
NASA NTRS1997

Simulation of a Controlled Airfoil with Jets

Numerical simulations of a two-dimensional airfoil, controlled by an applied moment in pitch and an airfoil controlled by jets, were investigated. These simulations couple the Reynolds-averaged Navier-Stokes equations and Euler's equations of rigid body motion, with an active control system. Controllers for both systems were designed to track altitude commands and were evaluated by simulating a closed-loop altitude step response using the coupled system. The airfoil controlled by a pitching moment used an optimal state feedback controller. A closed-loop simulation, of the airfoil with an applied moment, showed that the trajectories compared very well with quasi-steady aerodynamic theory, providing a measure of validation. The airfoil with jets used a controller designed by robust control methods. A linear plant model for this system was identified using open-loop data generated by the nonlinear coupled system. A closed-loop simulation of the airfoil with jets, showed good tracking of an altitude command. This simulation also showed oscillations in the control input as a result of dynamics not accounted for in the control design. This research work demonstrates how computational fluid dynamics, coupled with rigid body dynamics, and a control law can be used to prototype control systems in problematic nonlinear flight regimes.

Allan, Brian G.↗
NASA NTRS1960

Supersonic Aerodynamic Characteristics of a Low-Drag Aircraft Configuration having an Arrow Wing of Aspect Ratio 1.86 and a Body of Fineness Ratio 20

A free-flight rocket-propelled-model investigation was conducted at Mach numbers of 1.2 to 1.9 to determine the longitudinal and lateral aero-dynamic characteristics of a low-drag aircraft configuration. The model consisted of an aspect-ratio -1.86 arrow wing with 67.5 deg. leading-edge sweep and NACA 65A004 airfoil section and a triangular vertical tail with 60 deg. sweep and NACA 65A003 section in combination with a body of fineness ratio 20. Aerodynamic data in pitch, yaw, and roll were obtained from transient motions induced by small pulse rockets firing at intervals in the pitch and yaw directions. From the results of this brief aerodynamic investigation, it is observed that very slender body shapes can provide increased volumetric capacity with little or no increase in zero-lift drag and that body fineness ratios of the order of 20 should be considered in the design of long-range supersonic aircraft. The zero-lift drag and the drag-due-to-lift parameter of the test configuration varied linearly with Mach number. The maximum lift-drag ratio was 7.0 at a Mach number of 1.25 and decreased slightly to a value of 6.6 at a Mach number of 1.81. The optimum lift coefficient, normal-force-curve slope, lateral-force-curve slope, static stability in pitch and yaw, time to damp to one-half amplitude in pitch and yaw, the sum of the rotary damping derivatives in pitch and also in yaw, and the static rolling derivatives all decreased with an increase in Mach number. Values of certain rolling derivatives were obtained by application of the least-squares method to the differential equation of rolling motion. A comparison of the experimental and calculated total rolling-moment-coefficient variation during transient oscillations of the model indicated good agreement when the damping-in-roll contribution was included with the static rolling-moment terms.

Gillespie, Warren, Jr.↗
NASA NTRS2007

Estimating the Inertia Matrix of a Spacecraft

A paper presents a method of utilizing some flight data, aboard a spacecraft that includes reaction wheels for attitude control, to estimate the inertia matrix of the spacecraft. The required data are digitized samples of (1) the spacecraft attitude in an inertial reference frame as measured, for example, by use of a star tracker and (2) speeds of rotation of the reaction wheels, the moments of inertia of which are deemed to be known. Starting from the classical equations for conservation of angular momentum of a rigid body, the inertia-matrix-estimation problem is formulated as a constrained least-squares minimization problem with explicit bounds on the inertia matrix incorporated as linear matrix inequalities. The explicit bounds reflect physical bounds on the inertia matrix and reduce the volume of data that must be processed to obtain a solution. The resulting minimization problem is a semidefinite optimization problem that can be solved efficiently, with guaranteed convergence to the global optimum, by use of readily available algorithms. In a test case involving a model attitude platform rotating on an air bearing, it is shown that, relative to a prior method, the present method produces better estimates from few data.

Acikmese, Behcet↗
NASA NTRSDate not supplied

Explicit Discontinuous Galerkin Methods for Conservation Laws

The two explicit DG methods in this study are based on a ‘predictor-corrector’ formulation, the first introduced by Lörcher, Gassner, and Munz (2007, 2008) called space–time expansion discontinuous Galerkin or STE-DG scheme, and the second, introduced independently by the author (Huynh 2006, 2013) called the upwind moment scheme. The predictor step of the two methods is essentially identical using a Cauchy-Kovalevsky (CK) procedure, which involves no interaction of the data among neighboring cells. The corrector step also shares the same space-time integration formulation and is where interaction of the data among neighboring cells takes place; the difference, however, is in how the resulting space-time volume integral is estimated. As a consequence of the different estimates, for the case of advection in one spatial dimension (1D), the moment scheme has a CFL (Courant-Friedrichs-Lewy) condition of 1 for all p and is accurate to order 2p+1, i.e., it possesses the super accuracy property, whereas the STE-DG method has a more restrictive CFL condition and is accurate to the expected order of p+1. For 1D advection, compared with the CFL conditions of 1/(2p+1) of standard RK-DG (Runge-Kutta) scheme where space and time discretization are of the same order, the moment scheme allows a significantly larger time step size. It also turns out that the scheme yields a result identical to Van Leer’s scheme III (1977), which amounts to shifting the data a distance of advection corresponding to the time step and projecting the result onto the space of polynomial solutions. Contrary to Van Leer’s approach, however, the space-time ‘predictor-corrector’ formulation facilitates extensions to the case of systems of equations. Concerning 2D extensions, in the case of advection, when the flow is along the diagonal direction, the CFL conditions for the moment schemes become restrictive as will be shown by Fourier (Von Neumann) stability and accuracy analyses. Since the moment scheme employs the right Radau points as collocation points in time, the method is closely related to the implicit Radau IIA scheme, which is stable for any time step size. The role of Radau IIA in relieving stability restriction for these explicit DG schemes remains to be explored

Discontinuous Galerkin↗
NASA NTRS1992

Polarimetric signatures of a canopy of dielectric cylinders based on first and second order vector radiative transfer theory

Complete polarimetric signatures of a canopy of dielectric cylinders overlying a homogeneous half space are studied with the first and second order solutions of the vector radiative transfer theory. The vector radiative transfer equations contain a general nondiagonal extinction matrix and a phase matrix. The energy conservation issue is addressed by calculating the elements of the extinction matrix and the elements of the phase matrix in a manner that is consistent with energy conservation. Two methods are used. In the first method, the surface fields and the internal fields of the dielectric cylinder are calculated by using the fields of an infinite cylinder. The phase matrix is calculated and the extinction matrix is calculated by summing the absorption and scattering to ensure energy conservation. In the second method, the method of moments is used to calculate the elements of the extinction and phase matrices. The Mueller matrix based on the first order and second order multiple scattering solutions of the vector radiative transfer equation are calculated. Results from the two methods are compared. The vector radiative transfer equations, combined with the solution based on method of moments, obey both energy conservation and reciprocity. The polarimetric signatures, copolarized and depolarized return, degree of polarization, and phase differences are studied as a function of the orientation, sizes, and dielectric properties of the cylinders. It is shown that second order scattering is generally important for vegetation canopy at C band and can be important at L band for some cases.

Tsang, Leung↗
NASA NTRS1995

Computational analysis of forebody tangential slot blowing on the high alpha research vehicle

A numerical analysis of forebody tangential slot blowing as a means of generating side force and yawing moment is conducted using an aircraft geometry. The Reynolds-averaged, thin-layer, Navier-Stokes equations are solved using a partially flux-split, approximately-factored algorithm. An algebraic turbulence model is used to determine the turbulent eddy viscosity values. Solutions are obtained using both patched and overset grid systems. In the patched grid model, and actuator plane is used to introduce jet variables into the flow field. The overset grid model is used to model the physical slot geometry and facilitate modeling of the full aircraft configuration. A slot optimization study indicates that a short slot located close to the nose of the aircraft provided the most side force and yawing moment per unit blowing coefficient. Comparison of computed surface pressure with that obtained in full-scale wind tunnel tests produce good agreement, indicating the numerical method and grid system used in the study are valid. Full aircraft computations resolve the changes in vortex burst point due to blowing. A time-accurate full-aircraft solution shows the effect of blowing on the changes in the frequency of the aerodynamic loads over the vertical tails. A study of the effects of freestream Mach number and various jet parameters indicates blowing remains effective through the transonic Mach range. An investigation of the force onset time lag associated with forebody blowing shows the lag to be minimal. The knowledge obtained in this study may be applied to the design of a forebody tangential slot blowing system for use on flight aircraft.

Gee, Ken↗
NASA NTRS1944

On the general theory of thin airfoils for nonuniform motion

General thin-airfoil theory for a compressible fluid is formulated as boundary problem for the velocity potential, without recourse to the theory of vortex motion. On the basis of this formulation the integral equation of lifting-surface theory for an incompressible fluid is derived with the chordwise component of the fluid velocity at the airfoil as the function to be determined. It is shown how by integration by parts this integral equation can be transformed into the Biot-Savart theorem. A clarification is gained regarding the use of principal value definitions for the integral which occur. The integral equation of lifting-surface theory is used a s the starting point for the establishment of a theory for the nonstationary airfoil which is a generalization of lifting-line theory for the stationary airfoil and which might be called "lifting-strip" theory. Explicit expressions are given for section lift and section moment in terms of the circulation function, which for any given wing deflection is to be determined from an integral equation which is of the type of the equation of lifting-line theory. The results obtained are for airfoils of uniform chord. They can be extended to tapered airfoils. One of the main uses of the results should be that they furnish a practical means for the analysis of the aerodynamic span effect in the problem of wing flutter. The range of applicability of "lifting-strip" theory is the same as that of lifting-line theory so that its results may be applied to airfoils with aspect ratios as low as three.

BIBLIOGRAPHY-WINGS↗