Search NASA⌕ Search

SEARCH · Search NASA

Results for “Lagrangian function”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1982

Computer simulations of optimum boost and buck-boost converters

The development of mathematicl models suitable for minimum weight boost and buck-boost converter designs are presented. The facility of an augumented Lagrangian (ALAG) multiplier-based nonlinear programming technique is demonstrated for minimum weight design optimizations of boost and buck-boost power converters. ALAG-based computer simulation results for those two minimum weight designs are discussed. Certain important features of ALAG are presented in the framework of a comprehensive design example for boost and buck-boost power converter design optimization. The study provides refreshing design insight of power converters and presents such information as weight annd loss profiles of various semiconductor components and magnetics as a function of the switching frequency.

Rahman, S.↗
NASA NTRS1990

Dynamic Analysis and Control of Lightweight Manipulators with Flexible Parallel Link Mechanisms

The objective is the theoretical analysis and the experimental verification of dynamics and control of a two link flexible manipulator with a flexible parallel link mechanism. Nonlinear equations of motion of the lightweight manipulator are derived by the Lagrangian method in symbolic form to better understand the structure of the dynamic model. The resulting equation of motion have a structure which is useful to reduce the number of terms calculated, to check correctness, or to extend the model to higher order. A manipulator with a flexible parallel link mechanism is a constrained dynamic system whose equations are sensitive to numerical integration error. This constrained system is solved using singular value decomposition of the constraint Jacobian matrix. Elastic motion is expressed by the assumed mode method. Mode shape functions of each link are chosen using the load interfaced component mode synthesis. The discrepancies between the analytical model and the experiment are explained using a simplified and a detailed finite element model.

Lee, Jeh Won↗
NASA NTRS1993

Dynamics of the baryonic component in hierarchical clustering universes

I present self-consistent 3-D simulations of the formation of virialized systems containing both gas and dark matter in a flat universe. A fully Lagrangian code based on the Smoothed Particle Hydrodynamics technique and a tree data structure has been used to evolve regions of comoving radius 2-3 Mpc. Tidal effects are included by coarse-sampling the density of the outer regions up to a radius approx. 20 Mpc. Initial conditions are set at high redshift (z greater than 7) using a standard Cold Dark Matter perturbation spectrum and a baryon mass fraction of 10 percent (omega(sub b) = 0.1). Simulations in which the gas evolves either adiabatically or radiates energy at a rate determined locally by its cooling function were performed. This allows us to investigate with the same set of simulations the importance of radiative losses in the formation of galaxies and the equilibrium structure of virialized systems where cooling is very inefficient. In the absence of radiative losses, the simulations can be rescaled to the density and radius typical of galaxy clusters. A summary of the main results is presented.

Navarro, Julio↗
NASA NTRS1991

Numerical simulation of the reactive flow in advanced (HSR) combustors using KIVA-2

Recent work has been done with the goal of establishing ultralow emission aircraft gas turbine combustors. A significant portion of the effort is the development of three dimensional computational combustor models. The KIVA-II computer code which is based on the Implicit Continuous Eulerian Difference mesh Arbitrary Lagrangian Eulerian (ICED-ALE) numerical scheme is one of the codes selected by NASA to achieve these goals. This report involves a simulation of jet injection through slanted slots within the Rich burn/Quick quench/Lean burn (RQL) baseline experimental rig. The RQL combustor distinguishes three regions of combustion. This work specifically focuses on modeling the quick quench mixer region in which secondary injection air is introduced radially through 12 equally spaced slots around the mixer circumference. Steady state solutions are achieved with modifications to the KIVA-II program. Work currently underway will evaluate thermal mixing as a function of injection air velocity and angle of inclination of the slots.

Winowich, Nicholas S.↗
NASA NTRS2012

Progress Toward Affordable High Fidelity Combustion Simulations Using Filtered Density Functions for Hypersonic Flows in Complex Geometries

Significant progress has been made in the development of subgrid scale (SGS) closures based on a filtered density function (FDF) for large eddy simulations (LES) of turbulent reacting flows. The FDF is the counterpart of the probability density function (PDF) method, which has proven effective in Reynolds averaged simulations (RAS). However, while systematic progress is being made advancing the FDF models for relatively simple flows and lab-scale flames, the application of these methods in complex geometries and high speed, wall-bounded flows with shocks remains a challenge. The key difficulties are the significant computational cost associated with solving the FDF transport equation and numerically stiff finite rate chemistry. For LES/FDF methods to make a more significant impact in practical applications a pragmatic approach must be taken that significantly reduces the computational cost while maintaining high modeling fidelity. An example of one such ongoing effort is at the NASA Langley Research Center, where the first generation FDF models, namely the scalar filtered mass density function (SFMDF) are being implemented into VULCAN, a production-quality RAS and LES solver widely used for design of high speed propulsion flowpaths. This effort leverages internal and external collaborations to reduce the overall computational cost of high fidelity simulations in VULCAN by: implementing high order methods that allow reduction in the total number of computational cells without loss in accuracy; implementing first generation of high fidelity scalar PDF/FDF models applicable to high-speed compressible flows; coupling RAS/PDF and LES/FDF into a hybrid framework to efficiently and accurately model the effects of combustion in the vicinity of the walls; developing efficient Lagrangian particle tracking algorithms to support robust solutions of the FDF equations for high speed flows; and utilizing finite rate chemistry parametrization, such as flamelet models, to reduce the number of transported reactive species and remove numerical stiffness. This paper briefly introduces the SFMDF model (highlighting key benefits and challenges), and discusses particle tracking for flows with shocks, the hybrid coupled RAS/PDF and LES/FDF model, flamelet generated manifolds (FGM) model, and the Irregularly Portioned Lagrangian Monte Carlo Finite Difference (IPLMCFD) methodology for scalable simulation of high-speed reacting compressible flows.

Drozda, Tomasz G.↗
NASA NTRS1992

Numerical modeling for dilute and dense sprays

We have successfully implemented a numerical model for spray-combustion calculations. In this model, the governing gas-phase equations in Eulerian coordinate are solved by a time-marching multiple pressure correction procedure based on the operator-splitting technique. The droplet-phase equations in Lagrangian coordinate are solved by a stochastic discrete particle technique. In order to simplify the calculation procedure for the circulating droplets, the effective conductivity model is utilized. The k-epsilon models are utilized to characterize the time and length scales of the gas phase in conjunction with turbulent modulation by droplets and droplet dispersion by turbulence. This method entails random sampling of instantaneous gas flow properties and the stochastic process requires a large number of computational parcels to produce the satisfactory dispersion distributions even for rather dilute sprays. Two major improvements in spray combustion modelings were made. Firstly, we have developed a probability density function approach in multidimensional space to represent a specific computational particle. Secondly, we incorporate the Taylor Analogy Breakup (TAB) model for handling the dense spray effects. This breakup model is based on the reasonable assumption that atomization and drop breakup are indistinguishable processes within a dense spray near the nozzle exit. Accordingly, atomization is prescribed by injecting drops which have a characteristic size equal to the nozzle exit diameter. Example problems include the nearly homogeneous and inhomogeneous turbulent particle dispersion, and the non-evaporating, evaporating, and burning dense sprays. Comparison with experimental data will be discussed in detail.

Chen, C. P.↗
DOE OSTI2025

CoCoMET v1.0: a unified open-source toolkit for atmospheric object tracking and analysis

Advances in performance and analysis capabilities have accelerated the development of object tracking algorithms for atmospheric research. This has resulted in a growing number of studies using Lagrangian tracking techniques to analyze the evolution of atmospheric phenomena and the underlying processes. However, the increasing complexity and variety of tracking algorithms present a steep learning curve for new users and make it difficult for existing users to compare algorithm performance. We introduce CoCoMET (Community Cloud Model Evaluation Toolkit), an open-source toolkit that addresses these issues. CoCoMET simplifies the process of running multiple tracking algorithms simultaneously and analyzing objects in both model and observational datasets by specifying parameters in a single configuration file. It standardizes input data from different sources into a consistent format and unifies the tracking output across algorithms. CoCoMET enhances the functionality of existing tracking methods by calculating additional properties such as cell growth and dissipation rates, perimeter, surface area, convexity, and irregularity. In addition, CoCoMET includes a novel method for identifying mergers and splits in 2D and 3D tracks and supports the integration of Eulerian/stationary datasets external to the tracking data for process studies. Its potential utility is demonstrated through examples of model intercomparison, model evaluation against observations, and comparisons between tracking algorithms. Designed for open-source environments, CoCoMET will continue to expand with future releases, incorporating more input data types and tracking algorithms.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1988

Computation of two-phase shear-layer flow using an Eulerian-Lagrangian analysis

A combined Eulerian-Lagrangian analysis which combines a linearized block implicit Navier-Stokes analysis for the continuous phase with a Lagrangian analysis for the discrete phase has been utilized to simulate a two-phase shear-layer flow in a two-dimensional channel flow configuration. The motion of the particles is tracked in computational coordinate space, resulting in improved computational efficiency, and interphase coupling terms for the Eulerian analysis are computed from the instantaneous distribution of the particles rather than the trajectory information. The use of the implicit Navier-Stokes analysis for the continuous phase has made it possible to efficiently use a highly stretched mesh. A low Reynolds number form of the k-epsilon turbulence model was used with near-wall resolution, thus eliminating the need to use the 'wall-function' approach. The streamwise velocity profiles for the continuous and the discrete phases have been compared with experimental data for two test cases.

Sabnis, Jayant S.↗
NASA NTRS2012

A Regional CO2 Observing System Simulation Experiment Using ASCENDS Observations and WRF-STILT Footprints

Knowledge of the spatiotemporal variations in emissions and uptake of CO2 is hampered by sparse measurements. The recent advent of satellite measurements of CO2 concentrations is increasing the density of measurements, and the future mission ASCENDS (Active Sensing of CO2 Emissions over Nights, Days and Seasons) will provide even greater coverage and precision. Lagrangian atmospheric transport models run backward in time can quantify surface influences ("footprints") of diverse measurement platforms and are particularly well suited for inverse estimation of regional surface CO2 fluxes at high resolution based on satellite observations. We utilize the STILT Lagrangian particle dispersion model, driven by WRF meteorological fields at 40-km resolution, in a Bayesian synthesis inversion approach to quantify the ability of ASCENDS column CO2 observations to constrain fluxes at high resolution. This study focuses on land-based biospheric fluxes, whose uncertainties are especially large, in a domain encompassing North America. We present results based on realistic input fields for 2007. Pseudo-observation random errors are estimated from backscatter and optical depth measured by the CALIPSO satellite. We estimate a priori flux uncertainties based on output from the CASA-GFED (v.3) biosphere model and make simple assumptions about spatial and temporal error correlations. WRF-STILT footprints are convolved with candidate vertical weighting functions for ASCENDS. We find that at a horizontal flux resolution of 1 degree x 1 degree, ASCENDS observations are potentially able to reduce average weekly flux uncertainties by 0-8% in July, and 0-0.5% in January (assuming an error of 0.5 ppm at the Railroad Valley reference site). Aggregated to coarser resolutions, e.g. 5 degrees x 5 degrees, the uncertainty reductions are larger and more similar to those estimated in previous satellite data observing system simulation experiments.

Wang, James S.↗
NASA NTRS2002

Molten-Metal Droplet Deposition on a Moving Substrate in Microgravity: Aiding the Development of Novel Technologies for Microelectronic Assembly

Driven by advancements in microelectronics manufacturing, this research investigates the oblique (non-axisymmetric) impact of liquid-metal droplets on flat substrates. The problem of interest is relevant to the development of the novel technology of on-demand dispension (printing) of microscopic solder deposits for the surface mounting of microelectronic devices. The technology, known as solder jetting, features on-demand deposition of miniature solder droplets (30 to 120 microns in diameter) in very fine, very accurate patterns using techniques analogous to those developed for the ink-jet printing industry. Despite its promise, severe limitations exist currently with regards to the throughput rates of the technology; some of these limitations are largely due to the lack of the capability for reliable prediction of solder bump positioning and shapes, especially under ballistic deposition conditions where the droplet impact phenomena are inherently three-dimensional. The study consists of a theoretical and an experimental component. The theoretical work uses a finite element formulation to simulate numerically the non-axisymmetric (3-D) fluid mechanics and heat transfer phenomena of a liquid solder droplet impacting at an angle alpha on a flat substrate. The work focuses on the pre-solidification regime. The modeling of the most challenging fluid mechanics part of the process has been completed successfully. It is based upon the full laminar Navier-Stokes equations employing a Lagrangian frame of reference. Due to the large droplet deformation, the surface (skin) as well as the volumetric mesh have to be regenerated during the calculations in order to maintain the high accuracy of the numerical scheme. The pressure and velocity fields are then interpolated on the newly created mesh. The numerical predictions are being tested against experiments, for cases where wetting phenomena are not important. For the impact parameters used in the example shown (We = 2.38, Fr = 16300, Re = 157), the droplet rolls along the substrate, but its shape remains practically axisymmetric for all impact angles within the range from 0 to 60 deg. Interestingly, the substrate/droplet contact area during the recoiling phase of the impact is not a monotonically decreasing function of time. The experimental component of the research tests the numerical predictions and provides necessary input data (contact angles) for the theoretical model. The experiments are performed in microgravity (2.2s drop tower of the NASA GRC) in order to allow for the use of mm-size solder droplets, which make feasible the performance of accurate measurements, while maintaining similitude of the relevant fluid dynamic groups (Re, Fr, We, Ste). Preliminary oblique impact experiments have been performed using water droplets in normal gravity.

Megaridis, C. M.↗
NASA NTRS2004

Model of Mixing Layer With Multicomponent Evaporating Drops

A mathematical model of a three-dimensional mixing layer laden with evaporating fuel drops composed of many chemical species has been derived. The study is motivated by the fact that typical real petroleum fuels contain hundreds of chemical species. Previously, for the sake of computational efficiency, spray studies were performed using either models based on a single representative species or models based on surrogate fuels of at most 15 species. The present multicomponent model makes it possible to perform more realistic simulations by accounting for hundreds of chemical species in a computationally efficient manner. The model is used to perform Direct Numerical Simulations in continuing studies directed toward understanding the behavior of liquid petroleum fuel sprays. The model includes governing equations formulated in an Eulerian and a Lagrangian reference frame for the gas and the drops, respectively. This representation is consistent with the expected volumetrically small loading of the drops in gas (of the order of 10 3), although the mass loading can be substantial because of the high ratio (of the order of 103) between the densities of liquid and gas. The drops are treated as point sources of mass, momentum, and energy; this representation is consistent with the drop size being smaller than the Kolmogorov scale. Unsteady drag, added-mass effects, Basset history forces, and collisions between the drops are neglected, and the gas is assumed calorically perfect. The model incorporates the concept of continuous thermodynamics, according to which the chemical composition of a fuel is described probabilistically, by use of a distribution function. Distribution functions generally depend on many parameters. However, for mixtures of homologous species, the distribution can be approximated with acceptable accuracy as a sole function of the molecular weight. The mixing layer is initially laden with drops in its lower stream, and the drops are colder than the gas. Drop evaporation leads to a change in the gas-phase composition, which, like the composition of the drops, is described in a probabilistic manner

Bellan, Josette↗
DOE OSTI2025

Off-shell vertices in heavy particle effective theories and B → Dπℓν

We study the modifications to decay amplitudes in heavy to heavy semileptonic decays with multiple hadrons in the final state due to intermediate heavy hadrons being off-shell or having a finite width. Combining Heavy Hadron Chiral Perturbation Theory (HHχPT) with a BCFW on-shell factorization formula, we show that these effects induce O (1/M) corrections to the standard results computed in the narrow-width approximation and therefore are important in extracting form factors from data. A combination of perturbative unitarity, analyticity, and reparameterization invariance fully determine these corrections in terms of known Isgur-Wise functions without the need to introduce new form factors. In doing so, we develop a novel technique to compute the boundary term at complex infinity in the BCFW formula for theories with derivatively coupled scalars. While we have used the $\overline{B}$ → Dπℓν decay as an example, these techniques can generally be applied to effective field theories with (multiple) distinct reference vectors. Article PDF

Chiral Lagrangian↗
NASA NTRS1999

Detailed Studies on the Structure and Dynamics of Reacting Dusty Flows at Normal and Microgravity

Dusty reacting flows are of particular interest for a wide range of applications. Inert particles can alter the flammability and extinction limits of a combustible mixture. Reacting particles can release substantial amount of heat and can be used either for power generation or propulsion. Accumulation of combustible particles in air can result in explosions which, for example, can occur in grain elevators, during lumber milling and in mine galleries. Furthermore, inert particles are used as flow velocity markers in reacting flows, and their velocity is measured by non-intrusive laser diagnostic techniques. Despite their importance, dusty reacting flows have been less studied and understood compared to gas phase as well as sprays. The addition of solid particles in a flowing gas stream can lead to strong couplings between the two phases, which can be of dynamic, thermal, and chemical nature. The dynamic coupling between the two phases is caused by the inertia that causes the phases to move with different velocities. Furthermore, gravitational, thermophoretic, photophoretic, electrophoretic, diffusiophoretic, centrifugal, and magnetic forces can be exerted on the particles. In general, magnetic, electrophoretic, centrifugal, photophoretic, and diffusiophoretic can be neglected. On the other hand, thermophoretic forces, caused by steep temperature gradients, can be important. The gravitational forces are almost always present and can affect the dynamic response of large particles. Understanding and quantifying the chemical coupling between two phases is a challenging task. However, all reacting particles begin this process as inert particles, and they must be heated before they participate in the combustion process. Thus, one must first understand the interactions of inert particles in a combustion environment. The in-detail understanding of the dynamics and structure of dusty flows can be only advanced by considering simple flow geometries such as the opposed-jet, stagnation-type. In such configurations the imposed strain rate is well characterized, and the in-depth understanding of the details of the physico-chemical processes can be systematically obtained. A number of computational and experimental studies on spray and particle flows have been conducted in stagnation-type configurations. Numerically, the need for a hybrid Eulerian-Lagrangian approach has been identified by Continillo and Sirignano, and the use of such approach has allowed for the prediction of the phenomenon of droplet flow reversal. Gomez and Rosner have conducted a detailed study on the particle response in the opposed-jet configuration, and the particle thermophoretic diffusivities were determined experimentally. Sung, Law and co-workers have conducted numerical studies on the effect of strain rate and temperature gradients on the dynamics of inert particles, as a way of understanding potential errors in experimental LDV data that may arise from thermophoretic forces. This investigation is a combined experimental and numerical study on the details of reacting dusty flows. The specific tasks are: (1) Experimental determination of laminar flame speeds, and extinction strain rates of dusty flows at normal- and micro-gravity as functions of the particle type, particle initial diameter, particle initial number density, and gas phase chemical composition; (2) Detailed numerical simulation of the experiments. Results are compared with experiments and the adequacy of theoretical models is assessed; and (3) Provision of enhanced insight into the thermo-chemical coupling between the two phases.

Andac, M. Gurhan↗
NASA NTRS1984

Structure and Growth of the Marine Boundary Layer

LANDSAT visible imagery and a one-dimensional Lagrangian boundary layer model were used to hypothesize the nature and the development of the marine boundary layer during a winter episode of strong seaward cold air advection. Over-water heating and moistening of the cold, dry continental air is estimable from linear relations involving horizontal gradients of the near-surface air temperature and humidity. A line of enhanced convection paralleling the Atlantic U.S. coast from south of New York Bay to the vicinity of Virginia Beach, VA was attributed to stronger convergence at low levels. This feature was characterized as a mesoscale front. With the assistance of a three-dimensional mesoscale boundary layer model, initialized with data obtained from the MASEX, the marine boundary layer can be mapped over the entire Atlantic coastal domain and the evolution of the boundary layer can be studied as a function of different characteristics of important surface level forcings. The effects on boundary layer growth due to the magnitude and pattern of sea surface temperature, to the shape of the coastline, and to atmospheric conditions, such as the orientation of the prevailing wind are examined.

Mccumber, M.↗
NASA NTRS1984

Center of mass estimation in closed vortices - A verification in principle and practice

The problem of tracking closed mesoscale vortices using center of mass estimation techniques is studied. Three estimators are evaluated using data from a warm core Gulf Stream ring. The comparisons show that a method based on the intersection of perpendicular bisectors and one using a least-squares fit of a conic section perform comparably. The perpendicular bisector algorithm is used in conjunction with a Gaussian ring model and a star-shaped survey pattern to produce an expected error curve as a function of vortex translation, survey speed and vortex size. For typical ring parameters, center estimation is usually possible to within + or - 5 km. The feasibility of using differing data sets to construct a history of ring motion based on a coordinate system moving with the ring is also investigated. In this way, the validity of using satellite-derived data and drifter trajectories to estimate the center of mass of a mesoscale feature is assessed. The results of the analysis demonstrate that the location of the deeper structure of the ring and the surface expression are sufficiently well correlated to permit dynamically relevant calculations based on surface measurements. It is shown that satellite-derived data can be used to approximate the center of mass trajectory to within the error in the individual center estimates for the period analyzed. The Lagrangian-drifter-derived centers are offset from the center of mass trajectory in a manner consistent with kinematic arguments.

Hooker, S. B.↗
NASA NTRS1998

LSPRAY: Lagrangian Spray Solver for Applications With Parallel Computing and Unstructured Gas-Phase Flow Solvers

Sprays occur in a wide variety of industrial and power applications and in the processing of materials. A liquid spray is a phase flow with a gas as the continuous phase and a liquid as the dispersed phase (in the form of droplets or ligaments). Interactions between the two phases, which are coupled through exchanges of mass, momentum, and energy, can occur in different ways at different times and locations involving various thermal, mass, and fluid dynamic factors. An understanding of the flow, combustion, and thermal properties of a rapidly vaporizing spray requires careful modeling of the rate-controlling processes associated with the spray's turbulent transport, mixing, chemical kinetics, evaporation, and spreading rates, as well as other phenomena. In an attempt to advance the state-of-the-art in multidimensional numerical methods, we at the NASA Lewis Research Center extended our previous work on sprays to unstructured grids and parallel computing. LSPRAY, which was developed by M.S. Raju of Nyma, Inc., is designed to be massively parallel and could easily be coupled with any existing gas-phase flow and/or Monte Carlo probability density function (PDF) solver. The LSPRAY solver accommodates the use of an unstructured mesh with mixed triangular, quadrilateral, and/or tetrahedral elements in the gas-phase solvers. It is used specifically for fuel sprays within gas turbine combustors, but it has many other uses. The spray model used in LSPRAY provided favorable results when applied to stratified-charge rotary combustion (Wankel) engines and several other confined and unconfined spray flames. The source code will be available with the National Combustion Code (NCC) as a complete package.

Raju, Manthena S.↗
NASA NTRS2016

Pressure Fluctuations Induced by a Hypersonic Turbulent Boundary Layer

Direct numerical simulations (DNS) are used to examine the pressure fluctuations generated by a spatially-developed Mach 5.86 turbulent boundary layer. The unsteady pressure field is analyzed at multiple wall-normal locations, including those at the wall, within the boundary layer (including inner layer, the log layer, and the outer layer), and in the free stream. The statistical and structural variations of pressure fluctuations as a function of wall-normal distance are highlighted. Computational predictions for mean velocity pro les and surface pressure spectrum are in good agreement with experimental measurements, providing a first ever comparison of this type at hypersonic Mach numbers. The simulation shows that the dominant frequency of boundary-layer-induced pressure fluctuations shifts to lower frequencies as the location of interest moves away from the wall. The pressure wave propagates with a speed nearly equal to the local mean velocity within the boundary layer (except in the immediate vicinity of the wall) while the propagation speed deviates from the Taylor's hypothesis in the free stream. Compared with the surface pressure fluctuations, which are primarily vortical, the acoustic pressure fluctuations in the free stream exhibit a significantly lower dominant frequency, a greater spatial extent, and a smaller bulk propagation speed. The freestream pressure structures are found to have similar Lagrangian time and spatial scales as the acoustic sources near the wall. As the Mach number increases, the freestream acoustic fluctuations exhibit increased radiation intensity, enhanced energy content at high frequencies, shallower orientation of wave fronts with respect to the flow direction, and larger propagation velocity.

Lian Duan↗
NASA NTRS2012

Extension of the Time-Spectral Approach to Overset Solvers for Arbitrary Motion

Forced periodic flows arise in a broad range of aerodynamic applications such as rotorcraft, turbomachinery, and flapping wing configurations. Standard practice involves solving the unsteady flow equations forward in time until the initial transient exits the domain and a statistically stationary flow is achieved. It is often required to simulate through several periods to remove the initial transient making unsteady design optimization prohibitively expensive for most realistic problems. An effort to reduce the computational cost of these calculations led to the development of the Harmonic Balance method [1, 2] which capitalizes on the periodic nature of the solution. The approach exploits the fact that forced temporally periodic flow, while varying in the time domain, is invariant in the frequency domain. Expanding the temporal variation at each spatial node into a Fourier series transforms the unsteady governing equations into a steady set of equations in integer harmonics that can be tackled with the acceleration techniques afforded to steady-state flow solvers. Other similar approaches, such as the Nonlinear Frequency Domain [3,4,5], Reduced Frequency [6] and Time-Spectral [7, 8, 9] methods, were developed shortly thereafter. Additionally, adjoint-based optimization techniques can be applied [10, 11] as well as frequency-adaptive methods [12, 13, 14] to provide even more flexibility to the method. The Fourier temporal basis functions imply spectral convergence as the number of harmonic modes, and correspondingly number of time samples, N, is increased. Some elect to solve the equations in the frequency domain directly, while others choose to transform the equations back into the time domain to simplify the process of adding this capability to existing solvers, but each harnesses the underlying steady solution in the frequency domain. These temporal projection methods will herein be collectively referred to as Time-Spectral methods. Time-Spectral methods have demonstrated marked success in reducing the computational costs associated with simulating periodic forced flows, but have yet to be fully applied to overset or Cartesian solvers for arbitrary motion with dynamic hole-cutting. Overset and Cartesian grid methodologies are versatile techniques capable of handling complex geometry configurations in practical engineering applications, and the combination of the Time-Spectral approach with this general capability potentially provides an enabling new design and analysis tool. In an arbitrary moving-body scenario for these approaches, a Lagrangian body moves through a fixed Eulerian mesh and mesh points in the Eulerian mesh interior to the solid body are removed (cut or blanked), leaving a hole in the Eulerian mesh. During the dynamic motion some gridpoints in the domain are blanked and do not have a complete set of time-samples preventing a direct implementation of the Time-Spectral method. Murman[6] demonstrated the Time-Spectral approach for a Cartesian solver with a rigid domain motion, wherein the hole cutting remains constant. Similarly, Custer et al. [15, 16] used the NASA overset OVERFLOW solver and limited the amount of relative motion to ensure static hole-cutting and interpolation. Recently, Mavriplis and Mundis[17] demonstrated a qualitative method for applying the Time-Spectral approach to an unstructured overset solver for arbitrary motion. The goal of the current work is to develop a robust and general method for handling arbitrary motion with the Time-Spectral approach within an overset or Cartesian mesh method, while still approaching the spectral convergence rate of the original Time-Spectral approach. The viscous OVERFLOW solver will be augmented with the new Time-Spectral algorithm and the capability of the method for benchmark problems in rotorcraft and turbomachinery will be demonstrated. This abstract begins with a brief synopsis of the Time-Spectral approach for overset grids and provides details of e current approach to allow for arbitrary motion. Model problem results in one and two dimensions are included to demonstrate the viability of the method and the convergence properties. Section IV briefly outlines the implementation into the OVERFLOW solver, and the abstract closes with a description of the benchmark test cases which will be included in the final paper.

Leffell, Joshua Isaac↗