Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Analytic gradients”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Determining Optimal Magnetometer Configuration on MAGIS-100

Long-baseline atom interferometers such as the Matter-wave Atomic Gradiometer Interferometric Sensor (MAGIS-100) require stringent control and continuous characterization of background magnetic fields and spatial gradients to prevent systemic phase shifts that mimic ultralight dark matter or gravitational wave signatures. Because direct sensor placement within the ultra-high vacuum beam pipe is infeasible, in-situ magnetic field monitoring relies on external sensor arrays situated in the surrounding annular region. This work demonstrates a field reconstruction framework for a 5.3-meter MAGIS-100 modular section using finite-element Opera simulations. Transverse magnetic fields are expanded using a cylindrical multipole framework as informed by Fermilab’s Muon g-2 experiment, with magnetometer array configurations optimized via Fisher information matrix D-optimality. Inverting external sensor readings through a Gauss-Newton scheme recovers interior tube fields across distinct axial positions. In the discontinuity-averse uniform region (slice pair P4), the model achieves sub-noise-floor performance with a cross-validated root-mean-square error (RMSE) of $6.7227 \times 10^{-4}\text{ A/m}$ ($0.845\times$ sensor noise floor) and an interior field coefficient of variation of $1.71\%$. An elbow criterion in the Fisher bounds establishes $n_{\text{max}} = 2$ as the optimal multipole truncation order to prevent noise amplification from over-parameterization, with $n_{\text{max}} = 3$ (sextupole) order chosen for analysis to demonstrate further complexity and cross-pair comparison. Furthermore, analytical differentiation of the fitted multipole coefficients yields dense spatial maps of the transverse Jacobian gradient matrix $\nabla \mathbf{H}$ along with propagated $1\sigma$ uncertainty bounds across the beam region ($r \le 2.75\text{ in}$). This operational framework confirms that external magnetometer arrays can reliably monitor magnetic field uniformity and spatial gradients along the 100-meter flight path given appropriate sampling for any complexity order.

Appleby, Darwin [William Rainey Harper Coll.] (ORC↗
DOE OSTI2026

Determining Optimal Magnetometer Configuration on MAGIS-100

Long-baseline atom interferometers such as the Matter-wave Atomic Gradiometer Interferometric Sensor (MAGIS-100) require stringent control and continuous characterization of background magnetic fields and spatial gradients to prevent systemic phase shifts that mimic ultralight dark matter or gravitational wave signatures. Because direct sensor placement within the ultra-high vacuum beam pipe is infeasible, in-situ magnetic field monitoring relies on external sensor arrays situated in the surrounding annular region. This work demonstrates a field reconstruction framework for a 5.3-meter MAGIS-100 modular section using finite-element Opera simulations. Transverse magnetic fields are expanded using a cylindrical multipole framework as informed by Fermilab’s Muon g-2 experiment, with magnetometer array configurations optimized via Fisher information matrix D-optimality. Inverting external sensor readings through a Gauss-Newton scheme recovers interior tube fields across distinct axial positions. In the discontinuity-averse uniform region (slice pair P4), the model achieves sub-noise-floor performance with a cross-validated root-mean-square error (RMSE) of $6.7227 \times 10^{-4}\text{ A/m}$ ($0.845\times$ sensor noise floor) and an interior field coefficient of variation of $1.71\%$. An elbow criterion in the Fisher bounds establishes $n_{\text{max}} = 2$ as the optimal multipole truncation order to prevent noise amplification from over-parameterization, with $n_{\text{max}} = 3$ (sextupole) order chosen for analysis to demonstrate further complexity and cross-pair comparison. Furthermore, analytical differentiation of the fitted multipole coefficients yields dense spatial maps of the transverse Jacobian gradient matrix $\nabla \mathbf{H}$ along with propagated $1\sigma$ uncertainty bounds across the beam region ($r \le 2.75\text{ in}$). This operational framework confirms that external magnetometer arrays can reliably monitor magnetic field uniformity and spatial gradients along the 100-meter flight path given appropriate sampling for any complexity order.

Appleby, Darwin [William Rainey Harper Coll.] (ORC↗
DOE OSTI2026

In-Situ Magnetic Field Reconstruction in the MAGIS-100 Experiment

Long-baseline atom interferometers such as the Matter-wave Atomic Gradiometer Interferometric Sensor (MAGIS-100) require stringent control and continuous characterization of background magnetic fields and spatial gradients to prevent systemic phase shifts that mimic ultralight dark matter or gravitational wave signatures. Because direct sensor placement within the ultra-high vacuum beam pipe is infeasible, in-situ magnetic field monitoring relies on external sensor arrays situated in the surrounding annular region. This work demonstrates a field reconstruction framework for a 5.3-meter MAGIS-100 modular section using finite-element Opera simulations. Transverse magnetic fields are expanded using a cylindrical multipole framework as informed by Fermilab’s Muon g-2 experiment, with magnetometer array configurations optimized via Fisher information matrix D-optimality. Inverting external sensor readings through a Gauss-Newton scheme recovers interior tube fields across distinct axial positions. In the discontinuity-averse uniform region (slice pair P4), the model achieves sub-noise-floor performance with a cross-validated root-mean-square error (RMSE) of $6.7227 \times 10^{-4}\text{ A/m}$ ($0.845\times$ sensor noise floor) and an interior field coefficient of variation of $1.71\%$. An elbow criterion in the Fisher bounds establishes $n_{\text{max}} = 2$ as the optimal multipole truncation order to prevent noise amplification from over-parameterization, with $n_{\text{max}} = 3$ (sextupole) order chosen for analysis to demonstrate further complexity and cross-pair comparison. Furthermore, analytical differentiation of the fitted multipole coefficients yields dense spatial maps of the transverse Jacobian gradient matrix $\nabla \mathbf{H}$ along with propagated $1\sigma$ uncertainty bounds across the beam region ($r \le 2.75\text{ in}$). This operational framework confirms that external magnetometer arrays can reliably monitor magnetic field uniformity and spatial gradients along the 100-meter flight path given appropriate sampling for any complexity order.

Appleby, Darwin [William Rainey Harper Coll.; Ferm↗
NASA NTRS2023

Finite Element Model Tuning using Analytical Sensitivity Values

This paper presents an efficient approach for tuning finite element models to match the measured ground vibrationtest and proof test data. Frequencies, mode shapes, total weight, location of the center of gravity, and staticdeformation computed from the finite element model are matched to the measured data. The model tuningprocedure used in this work is based on solving an optimization problem in which the errors for the consideredmetrics, between the finite element prediction and the measured data, are minimized. Analytical sensitivity values ofperformance indices are computed using the NASTRAN-generated sensitivity values together with the in-housecomputer programs, which allow for faster computational time and the use of gradient-based optimizers. The methodis applied to the Aerostructures Test Wing 4 model. The study shows that the military standard and the NASAstandard for comparing analytical and experimental modal data are all satisfied. The final finite element modelcorrelates well with the test data. The flutter speed decreases by 8.91% after model tuning compared with the original Aerostructures Test Wing 4 design.

Chan-gi Pak↗
DOE OSTI2021

Simple, general, realistic, robust, analytic tokamak equilibria. Part 2. Pedestals and flow

Part 1 described a wide range of analytic tokamak equilibria modelling smooth limiter surfaces, double- and single-null divertor surfaces, arbitrary aspect ratio, elongation, triangularity and beta. Part 2 generalizes the analysis to further include edge pedestals and toroidal flow. Specifically, edge pedestals are allowed in the pressure, pressure gradient and toroidal current density. Also, an edge-localized contribution to the bootstrap current is treated. In terms of flow, analytic solutions are obtained for two cases: a $\gamma = 2$ adiabatic and a $\gamma = \infty $ incompressible energy conservation relation.

Physics↗
DOE OSTI2024

Write-up on the analytical formalism for DTL Synchrotron and Betatron frequencies estimates

Appendix A contains tables of the DTL tanks geometry, generated using PARMILA code with input parameters of initial and final energy, accelerating gradient, synchronouse phase, and using DTLfish code (part of the SUPERFISH/POISSON package of codes). Details of the design and generation od the DTL geometry will be presented in a separate report. After generation of the DTL geometry and initial beam dynamics simulations, the analytic estimates of the synchrotron and betatron oscillations, as well as presence of the parametric resonance of first order were done. In present report the analytic background of those estimates, as well as the results of these are presented.

43 PARTICLE ACCELERATORS↗
DOE OSTI2026

Molecular concentration field design using closed-form steady-state solutions

Control over spatial concentration fields represents a fundamental challenge in designing synthetic biological systems and programmable soft materials. While nature creates morphogen gradients that orchestrate complex developmental processes, synthetic approaches have largely relied on empirical optimization and computationally intensive simulations. Here, we present an analytical framework for steady-state concentration fields generated by finite-sized localized sources in diffusion–degradation systems and derive closed-form solutions for one-, two-, and three-dimensional geometries. By expressing these solutions in dimensionless form, we show that gradient steepness and spatial structure are organized by the Thiele modulus, which captures the competition between diffusion and degradation length scales. The analysis reveals distinct design regimes: in degradation-dominated systems, gradient shape is governed by exponential decay and becomes dimension-independent, whereas in diffusion-dominated systems, gradient magnitude and extent follow dimension-dependent power-law scaling. Building on these results, we introduce a quantitative design strategy that uses threshold-based criteria to program concentration ranges by tuning physically accessible parameters, most directly the production rate, while holding transport and degradation properties fixed. Comparisons with numerical solutions and reported experimental systems demonstrate consistency with the predicted scaling behavior. Together, this work provides a generalizable and physically transparent framework for designing steady-state concentration fields in synthetic biological and soft matter systems, enabling predictive control of gradient-mediated organization without reliance on extensive numerical optimization.

Kim, Dong Woo [Johns Hopkins University, Baltimore↗
NASA NTRS2003

Analysis of Radial Segregation in Directionally Solidified Hg(0.89)Mn(0.11)Te

Bridgman growth experiments were performed on Hg(0.89)Mn(0.11)Te (MMT) to determine the extent of radial Manganese segregation during directional solidification. MMT crystals were directionally solidified at rates of 0.09 and 0.18 p d s and in axial thermal gradients of 83 and 68"C/cm. Wavelength Dispersive Spectroscopy (WDS) and Fourier Transform Infra-Red (FTIR) analytical techniques were used to determine the radial homogeneity in all boules and the deflection of the solid-liquid interface (SLI) in two boules that were rapidly quenched after 5 to 6 cm of directional solidification. For all growth runs, the measured radial coinpositional variations were on the order of 0.01 molar percent MnTe in the steady state region of growth. Comparison of the measured radial compositional results of the crystals to predicted values in the diffusion-limited regime indicate a strong influence of convection near the solid-liquid interface. This conclusion is supported by the weak influence of the translation rates and axial thermal gradients utilized in this study upon radial compositional homogeneity.

Price, M. W.↗
NASA NTRS1992

Partially coherent scattering in stellar chromospheres. II - The first-order escape probability method. III - A second-order escape probability method

Approximate analytic expressions are derived for resonance-line wing diagnostics, accounting for frequency redistribution effects, for homogeneous slabs, and slabs with a constant Planck function gradient. Resonance-line emission profiles from a simplified conceptual standpoint are described in order to elucidate the basic physical parameters of the line-forming layers prior to the performance of detailed numerical calculations. An approximate analytic expression is derived for the dependence on stellar surface gravity of the location of the Ca II and Mg II resonance-line profile peaks. An approximate radiative transfer equation using generalized second-order escape probabilities, applicable even in the presence of nearly coherent scattering in the damping wings of resonance lines, is derived. Approximate analytic solutions that can be applied in special regimes and achieve good agreement with accurate numerical results are found.

Gayley, K. G.↗
DOE OSTI2025

A New Framework for the Attribution of Air‐Sea CO 2 Exchange

Abstract The air‐sea transfer of carbon dioxide can be viewed as a dynamic system through which atmospheric and oceanic processes push surface waters away from thermodynamic equilibrium, while diffusive gas transfer pulls them back toward local equilibrium. These push/pull processes drive significant sub‐seasonal, seasonal, and interannual variability in air‐sea carbon fluxes, the quantification of which is critical both for diagnosing the ocean response to fossil fuel emissions and for attempts to mitigate anthropogenic climate disruption through intentional modification of surface ocean biogeochemistry. In this study, we present a new approach for attributing air‐sea carbon fluxes to specific mechanisms. The new framework is first applied to a two‐box ocean nutrient and carbon cycle model as an illustrative example. Next, outputs from a regional eddy‐resolving model of the Southern Ocean are analyzed. The roles of multiple physical and biogeochemical processes are identified. The decomposition of the seasonal air‐sea carbon flux shows the dominant role of biological carbon pumps that are partially compensated by the transport convergence. Finally, the framework is used to diagnose the response to mesoscale iron and alkalinity release, explicitly quantifying transport feedback and eventual impacts on net air‐sea carbon flux. Ocean carbon transport has divergent influences between iron and alkalinity release, due to opposing near‐surface gradients of dissolved inorganic carbon. We suggest that our attribution framework may be a useful analytical technique for monitoring natural ocean carbon fluxes and quantifying the impacts of human intervention on the ocean carbon cycle.

Ito, Takamitsu [School of Earth and Atmospheric Sc↗
NASA NTRS1978

The accuracy of the narrow seal approximation in analyzing radial face seals

The accuracy of the narrow seal approximation which enables closed form analytical solutions for radial face seals is examined in this paper. Both hydrostatic and hydrodynamic effects in a flat, misaligned seal are considered. Analytical results obtained from a simplified Reynolds equation, based on the neglect of circumferential pressure gradient and seal curvature, are compared with accurate results from numerical solution of the complete Reynolds equation. The agreement between the approximate and accurate solutions is quite reasonable over a wide range of a seal's inner-to-outer radius ratio. For radius ratios greater than 0.8, the accuracy of the narrow seal approximation is better than 1 percent over most of the range of angular misalignment.

Etsion, I.↗
NASA NTRS1985

A space-marching method for incompressible Navier-Stokes equations

This paper deals with the development of a space-marching method for incompressible flows. The method solves the continuity and momentum equations as a coupled system at each streamwise station. The character of the system of equations has been changed from elliptic to hyperbolic/parabolic in order to enable the equations to be marched in space. The present method has many advantages compared to the existing parabolic or space-marching methods for incompressible flow: (1) it avoids the solution of Poisson equations, (2) it conserves the mass flow with no additional computation, (3) it does not require the specification of an assumed pressure field when used in the prediction of duct flows. The present method can capture strong secondary velocities and strong transverse pressure gradients. Predictions of the flow through straight and curved ducts are in good agreement with analytical and experimental results.

Pouagare, M.↗
NASA NTRS1988

Cooling flows and the stability of radio jets

This paper considers what happens to radio jets which, propagating out from the nucleus of an elliptical galaxy, encounters a point at which the ambient gas pressure increases very rapidly. The pressure variation within the sonic radius of a cooling flow is calculated, and a semianalytical calculation of the stability of a jet which encounters a rapidly increasing pressure gradient is given. The results of numerical hydrodynamic simulations which test the results of the analytic calculations are shown. The case where the jet encounters the sonic region at an oblique angle is considered, and it is suggested that this may provide an explanation for jets which are bent on small scales. These results are applied to radio sources associated with central dominant cluster galaxies and with more normal elliptical galaxies.

Soker, Noam↗
NASA NTRS1995

A deterministic particle method for one-dimensional reaction-diffusion equations

We derive a deterministic particle method for the solution of nonlinear reaction-diffusion equations in one spatial dimension. This deterministic method is an analog of a Monte Carlo method for the solution of these problems that has been previously investigated by the author. The deterministic method leads to the consideration of a system of ordinary differential equations for the positions of suitably defined particles. We then consider the time explicit and implicit methods for this system of ordinary differential equations and we study a Picard and Newton iteration for the solution of the implicit system. Next we solve numerically this system and study the discretization error both analytically and numerically. Numerical computation shows that this deterministic method is automatically adaptive to large gradients in the solution.

Mascagni, Michael↗
NASA NTRS2012

Probable Rotation States of Rocket Bodies in Low Earth Orbit

In order for Active Debris Removal to be accomplished, it is critically important to understand the probable rotation states of orbiting, spent rocket bodies. As compared to the question of characterizing small unresolved debris, in this problem there are several advantages: (1) objects are of known size, mass, shape and color, (2) they have typically been in orbit for a known period of time, (3) they are large enough that resolved images may be obtainable for verification of predicted orientation, and (4) the dynamical problem is simplified to first order by largely cylindrical symmetry. It is also nearly certain for realistic rocket bodies that internal friction is appreciable in the case where residual liquid or, to a lesser degree, unconsolidated solid fuels exist. Equations of motion have been developed for this problem in which internal friction as well as torques due to solar radiation, magnetic induction, and gravitational gradient are included. In the case of pure cylindrical symmetry, the results are compared to analytical predictions patterned after the standard approach for analysis of symmetrical tops. This is possible because solar radiation and gravitational torques may be treated as conservative. Agreement between results of both methods ensures their mutual validity. For monotone symmetric cylinders, solar radiation torque vanishes if the center of mass resides at the geometric center of the object. Results indicate that in the absence of solar radiation effects, rotation states tend toward an equilibrium configuration in which rotation is about the axis of maximum inertia, with the axis of minimum inertia directed toward the center of the earth. Solar radiation torque introduces a modification to this orientation. The equilibrium state is asymptotically approached within a characteristic timescale given by a simple ratio of relevant characterizing parameters for the body in question. Light curves are simulated for the expected asymptotic final rotation states of model objects, and these are compared to data derived from physical models of the same objects, tested in the Optical Measurements Center at JSC. Comparison to relevant light curves from actual orbiting rocket bodies are also performed, and diagnostic features of such curves are examined.

Ojakangas, Gregory W.↗
DOE OSTI2023

Stabilizing Ni-rich layered cathode for high-voltage operation through hierarchically heterogeneous doping with concentration gradient

High-nickel LiNi x Mn y Co 1-x-y O 2 (NMC) cathodes have demonstrated superior energy density, yet their stability is compromised under high voltage conditions. To address this, here we propose a strategy of heterogeneous doping with a concentration gradient, specifically through Sr–Zr co-modification. We synthesized Ni-rich NMC particles featuring several micron-sized secondary particles composed of micron-sized primary grains. This design aims to harness the structural robustness of single-crystalline grains and the favorable diffusion kinetics of polycrystalline secondary particles. Systematic characterization using a combination of electrochemical measurements and synchrotron analytics reveals an intriguing pattern of hierarchically heterogeneous Sr–Zr co-doping. It demonstrates a depth-dependent concentration gradient at the secondary particle level and competing dopant segregation over the buried grain boundaries. This unique characteristic creates opportunities for enhancing battery performance, particularly by optimizing precursors and implementing advanced modulation techniques. We also investigate the dissolution and precipitation of the cathode's transition metal cations upon high-voltage cycling. These insights suggest that a tailored compositional variation can be a viable approach to effectively design the next-generation high-Ni NMC cathode materials for high-voltage lithium batteries.

36 MATERIALS SCIENCE↗