Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Computing continuum”

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 703 records · Page 39

𝐹 𝐾 ⁡/𝐹 𝜋 from Möbius domain-wall fermions solved on gradient-flowed HISQ ensembles

We report the results of a lattice quantum chromodynamics calculation of 𝐹 𝐾 ⁡/𝐹 𝜋 using Möbius domain-wall fermions computed on gradient-flowed 𝑁 𝑓 =2 +1 +1 highly improved staggered quark (HISQ) ensembles. The calculation is performed with five values of the pion mass ranging from 130 ≲ 𝑚 𝜋 ≲ 400 MeV , four lattice spacings of 𝑎 ∼ 0.15, 0.12, 0.09 and 0.06 fm and multiple values of the lattice volume. The interpolation/extrapolation to the physical pion and kaon mass point, the continuum, and infinite volume limits are performed with a variety of different extrapolation functions utilizing both the relevant mixed-action effective field theory expressions as well as discretization-enhanced continuum chiral perturbation theory formulas. We find that the 𝑎 ∼ 0.06 fm ensemble is helpful, but not necessary to achieve a subpercent determination of 𝐹 𝐾 ⁡/𝐹 𝜋 . We also include an estimate of the strong isospin breaking corrections and arrive at a final result of 𝐹 $\hat{K}$ + ⁡ /𝐹 $\hat{𝜋}$ + = 1.1942⁢(45) with all sources of statistical and systematic uncertainty included. This is consistent with the Flavour Lattice Averaging Group average value, providing an important benchmark for our lattice action. Combining our result with experimental measurements of the pion and kaon leptonic decays leads to a determination of |𝑉 𝑢⁢𝑠 |/|𝑉 𝑢⁢𝑑 | = 0.2311⁢(10).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Rarefield-Flow Shuttle Aerodynamics Flight Model

A model of the Shuttle Orbiter rarefied-flow aerodynamic force coefficients has been derived from the ratio of flight acceleration measurements. The in-situ, low-frequency (less than 1Hz), low-level (approximately 1 x 10(exp -6) g) acceleration measurements are made during atmospheric re-entry. The experiment equipment designed and used for this task is the High Resolution Accelerometer Package (HiRAP), one of the sensor packages in the Orbiter Experiments Program. To date, 12 HiRAP re-entry mission data sets spanning a period of about 10 years have been processed. The HiRAP-derived aerodynamics model is described in detail. The model includes normal and axial hypersonic continuum coefficient equations as function of angle of attack, body-flap deflection, and elevon deflection. Normal and axial free molecule flow coefficient equations as a function of angle of attack are also presented, along with flight-derived rarefied-flow transition bridging formulae. Comparisons are made between the aerodynamics model, data from the latest Orbiter Operational Aerodynamic Design Data Book, applicable computer simulations, and wind-tunnel data.

Blanchard, Robert C.↗

Scientific machine learning for modeling and simulating complex fluids

The formulation of rheological constitutive equations—models that relate internal stresses and deformations in complex fluids—is a critical step in the engineering of systems involving soft materials. While data-driven models provide accessible alternatives to expensive first-principles models and less accurate empirical models in many engineering disciplines, the development of similar models for complex fluids has lagged. The diversity of techniques for characterizing non-Newtonian fluid dynamics creates a challenge for classical machine learning approaches, which require uniformly structured training data. Consequently, early machine-learning based constitutive equations have not been portable between different deformation protocols or mechanical observables. Here, we present a data-driven framework that resolves such issues, allowing rheologists to construct learnable models that incorporate essential physical information, while remaining agnostic to details regarding particular experimental protocols or flow kinematics. These scientific machine learning models incorporate a universal approximator within a materially objective tensorial constitutive framework. By construction, these models respect physical constraints, such as frame-invariance and tensor symmetry, required by continuum mechanics. We demonstrate that this framework facilitates the rapid discovery of accurate constitutive equations from limited data and that the learned models may be used to describe more kinematically complex flows. This inherent flexibility admits the application of these “digital fluid twins” to a range of material systems and engineering problems. We illustrate this flexibility by deploying a trained model within a multidimensional computational fluid dynamics simulation—a task that is not achievable using any previously developed data-driven rheological equation of state.

Science & Technology - Other Topics↗

Toward a Ginsparg-Wilson lattice Hamiltonian

To address quantum computation of quantities in quantum chromodynamics (QCD) for which chiral symmetry is important, it would be useful to have the Hamiltonian for a fermion satisfying the Ginsparg-Wilson (GW) equation. I work with a solution to the GW equation which is fractional linear in time derivatives. The resulting Hamiltonian is nonlocal and has ghosts, but is free of doublers and has the correct continuum limit. This construction works in general odd spatial dimensions, and I provide an explicit expression for the Hamiltonian in one spatial dimension. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

Rarefied-flow shuttle aerodynamics flight model

A model of the Shuttle Orbiter rarefied-flow aerodynamic force coefficients has been derived from the ratio of flight acceleration measurements. The in-situ, low frequency, low level acceleration measurements are made during atmospheric reentry. The experiment equipment designed and used for this task is the High Resolution Accelerometer Package (HiRAP), one of the sensor packages in the Orbiter Experiments Program. To date, 12 HiRAP reentry mission data sets spanning a period of about 10 years have been processed. The HiRAP derived aerodynamics model is described in detail. The model includes normal and axial hypersonic continuum coefficient equations as functions of angle-of-attack, body flap deflection, and elevon deflection. Normal and axial free molecule flow coefficient equations as a function of angle-of-attack are also presented, along with flight derived rarefied-flow transition bridging formulas. Comparisons are made between the aerodynamics model and data from the Orbiter Operational Aerodynamic Design Data Book, and applicable computer simulations and wind-tunnel data.

Blanchard, Robert C.↗

Logarithmic potential with super-super-exponential kink profiles and tails

In this work, we consider a novel one dimensional model of a logarithmic potential which has super-super-exponential kink profiles as well as kink tails. We provide analytic kink solutions of the model—it has 3 kinks, 3 mirror kinks and the corresponding antikinks. While some of the kink tails are super-super-exponential, some others are super-exponential whereas the remaining ones are exponential. The linear stability analysis reveals that there is a gap between the zero mode and the onset of continuum. Finally, we compare this potential and its kink solutions with those of very high order field theories harboring seven degenerate minima and their attendant kink solutions, specifically φ 14 , φ 16 and φ 18 .

97 MATHEMATICS AND COMPUTING↗

Laboratory molecular spectroscopy

The precision required in making spectroscopic measurements is discussed. Remarks are directed specifically to vibration-rotation spectra rather than continuum absorptions. The ultimate precision that is required for line positions is related to the width of the lines which may be no narrower than the Doppler width. The spectroscopic methods considered are those which are of the most general value to the astronomers; those which acquire and can handle large volumes of spectra in digital form, or in a form which is compatible with computer analysis, and in a form which is at least internally consistent. The use of dye laser, grating instruments, and the most versatile instrument for laboratory spectroscopy, the Fourier transform spectrometer is discussed.

Margolis, J.↗

Molecular Hydrogen in Diffuse Interstellar Clouds of Arbitrary Three-Dimensional Geometry

We have constructed three-dimensional models for the equilibrium abundance of molecular hydrogen in diffuse interstellar clouds of arbitrary geometry that are illuminated by ultraviolet radiation. The position-dependent photodissociation rate of H2 in such clouds was computed with a 26 ray approximation to model the attenuation of the incident ultraviolet radiation field by dust and by H2 line absorption. We have applied our modeling technique to the isolated diffuse cloud G236+39, assuming that the cloud has a constant density and that the thickness of the cloud along the line of sight is at every point proportional to the 100 micron continuum intensity measured by IRAS. We find that our model can successfully account for observed variations in the ratio of 100 micron continuum intensity to H I column density, with larger values of that ratio occurring along lines of sight in which the molecular hydrogen fraction is expected to be the largest. Using a standard chi square analysis to assess the goodness of fit of our models, we find (at the 60 a level) that a three-dimensional model is more successful in matching the observational data than a one-dimensional model in which the geometrical extent of the cloud along the line of sight is assumed to be much smaller than its extent in the plane of the sky. If D is the distance to G236 + 39, and given standard assumptions about the rate of grain-catalyzed H2 formation, we find that the cloud has an extent along the line of sight that is 0.9 +/- 0.1 times its mean extent projected onto the plane of the sky and a gas density of 53 +/- 8 (100 pc/D) H nuclei/cc and is illuminated by a radiation field of 1.1 +/- 0.2 (100 pc/D) times the mean interstellar radiation field. The derived 100 micron emissivity per nucleon is 1.13 +/- 0.06 x 10(exp -20) MJy/sr sq cm.

Spaans, Marco↗

A new continuum model for suspensions of gyrotactic micro-organisms

A new continuum model is formulated for dilute suspensions of swimming micro-organisms with asymmetric mass distributions. Account is taken of randomness in a cell's swimming direction, p, by postulating that the probability density function for p satisfies a Fokker-Planck equation analogous to that obtained for colloid suspensions in the presence of rotational Brownian motion. The deterministic torques on a cell, viscous and gravitational, are balanced by diffusion, represented by an isotropic rotary diffusivity Dr, which is unknown a priori, but presumably reflects stochastic influences on the cell's internal workings. When the Fokker-Planck equation is solved, macroscopic quantities such as the average cell velocity Vc, the particle diffusivity tensor D and the effective stress tensor sigma can be computed; Vc and D are required in the cell conservation equation, and sigma in the momentum equation. The Fokker-Planck equation contains two dimensionless parameters, lambda and epsilon; lambda is the ratio of the rotary diffusion time Dr-1 to the torque relaxation time B (balancing gravitational and viscous torques), while epsilon is a scale for the local vorticity or strain rate made dimensionless with B. In this paper we solve the Fokker-Planck equation exactly for epsilon = 0 (lambda arbitrary) and also obtain the first-order solution for small epsilon. Using experimental data on Vc and D obtained with the swimming alga, Chlamydomonas nivalis, in the absence of bulk flow, the epsilon = 0 results can be used to estimate the value of lambda for that species (lambda approximately 2.2; Dr approximately 0.13 s-1). The continuum model for small epsilon is then used to reanalyse the instability of a uniform suspension, previously investigated by Pedley, Hill & Kessler (1988). The only qualitatively different result is that there no longer seem to be circumstances in which disturbances with a non-zero vertical wavenumber are more unstable than purely horizontal disturbances. On the way, it is demonstrated that the only significant contribution to sigma, other than the basic Newtonian stress, is that derived from the stresslets associated with the cells' intrinsic swimming motions.

NASA Discipline Number 28-20↗

Low-Velocity Impact Analyses of Composite Panels Using Cohesive Zone Modeling and Continuum Damage Mechanics

A parametric study was conducted to gage the sensitivity of simulations of low-velocity impacts on laminated composite panels to various finite element model and material model features. The material model used in this parametric study was previously verified and validated for conducting progressive damage analyses of post-buckled composite structures subjected to quasi-static loading conditions. Based on the results of this parametric study, updates have been made to a progressive damage analysis material model and the best practices for its application. Experimental results are compared to computational results obtained with and without the updates to the models. Both thermoset and thermoplastic fiber-reinforced materials are considered. The resulting enhancements to the material model are presented in detail.

Frank A Leone↗

Low-Velocity Impact Analyses of Composite Panels Using Cohesive Zone Modeling and Continuum Damage Mechanics

A parametric study was conducted to gage the sensitivity of simulations of low-velocity impacts on laminated composite panels to various finite element model and material model features. The material model used in this parametric study was previously verified and validated for conducting progressive damage analyses of post-buckled composite structures subjected to quasi-static loading conditions. Based on the results of this parametric study, updates have been made to a progressive damage analysis material model and the best practices for its application. Experimental results are compared to computational results obtained with and without the updates to the models. Both thermoset and thermoplastic fiber-reinforced materials are considered. The resulting enhancements to the material model are presented in detail.

Frank Leone↗

Optimal 1D Ly α forest power spectrum estimation – I. DESI-lite spectra

ABSTRACT The 1D Ly α forest flux power spectrum P1D is sensitive to scales smaller than a typical galaxy survey, and hence ties to the intergalactic medium’s thermal state, suppression from neutrino masses, and new dark matter models. It has emerged as a competitive framework to study new physics, but also has come with various challenges and systematic errors in analysis. In this work, we revisit the optimal quadratic estimator for P1D, which is robust against the relevant problems such as pixel masking, time evolution within spectrum, and quasar continuum errors. We further improve the estimator by introducing a fiducial power spectrum, which enables us to extract more information by alleviating the discreteness of band powers. We meticulously apply our method to synthetic Dark Energy Spectroscopic Instrument (DESI) spectra and demonstrate how the estimator overcomes each challenge. We further apply an optimization scheme that approximates the Fisher matrix to three elements per row and reduces computation time by 60 per cent. We show that we can achieve per cent precision in P1D with 5-yr DESI data in the absence of systematics and provide forecasts for different spectral qualities.

79 ASTRONOMY AND ASTROPHYSICS↗

Continuum Model Development for Flow and Transport in Electrochemical Systems

Stanford University (Subcontractor) shall provide the services of qualified multiscale modelers to perform tasks that contribute to reaching the objectives of the LDRD project entitled Automated and Accelerated Continuum Model Development for Electrochemical Systems (Tracking number 24-ERD-051). These tasks relate to the development, deployment, and validation of multiscale models relevant to flow and transport in electrochemical systems.

42 ENGINEERING↗

Dissociative recombination of O2(+), NO(+) and N2(+)

A new L(2) approach for the calculation of the threshold molecular capture width needed for the determination of DR cross sections was developed. The widths are calculated with Fermi's golden rule by substituting Rydberg orbitals for the free electron continuum coulomb orbital. It is shown that the calculated width converges exponentially as the effective principal quantum number of the Rydberg orbital increases. The threshold capture width is then easily obtained. Since atmospheric recombination involves very low energy electrons, the threshold capture widths are essential to the calculation of DR cross sections for the atmospheric species studied here. The approach described makes use of bound state computer codes already in use. A program that collects width matrix elements over CI wavefunctions for the initial and final states is described.

Guberman, S. L.↗

Shape design sensitivity analysis and optimal design of structural systems

The material derivative concept of continuum mechanics and an adjoint variable method of design sensitivity analysis are used to relate variations in structural shape to measures of structural performance. A domain method of shape design sensitivity analysis is used to best utilize the basic character of the finite element method that gives accurate information not on the boundary but in the domain. Implementation of shape design sensitivty analysis using finite element computer codes is discussed. Recent numerical results are used to demonstrate the accuracy obtainable using the method. Result of design sensitivity analysis is used to carry out design optimization of a built-up structure.

Choi, Kyung K.↗

Water vapor absorption coefficients in the 8-13-micron spectral region - A critical review

Measurements of water vapor absorption coefficients in the thermal IR atmospheric window (8-13 microns) during the past 20 years obtained by a variety of techniques are reviewed for consistency and compared with computed values based on the AFGL spectral data tapes. The methods of data collection considered were atmospheric long path absorption with a CO2 laser or a broadband source and filters, a White cell and a CO2 laser or a broadband source and a spectrometer, and a spectrophone with a CO2 laser. Advantages and disadvantages of each measurement approach are given as a guide to further research. Continuum absorption has apparently been measured accurately to about the 5-10 percent level in five of the measurements reported.

Grant, William B.↗

Statistical discrete particle simulation studies

A particle simulation code was developed that is suited to Multiple Instruction Multiple Data multiprocessor computers. The resulting code is designed to be portable between a large class of multiprocessor computer architectures, including hypercubes (e.g., Intel iPSC/860), shared memory machines (e.g., Cray, SGI, DASH), and simple uniprocessors (e.g., SUN, VAX, IBM PC). This code provides a powerful general framework that can be adapted to specific applications, such as thermomechanical modeling, general 3-D geometry support, or visualization support via a codeveloped CPlot data reduction and visualization system. Chemistry models were implemented and validated with comparison to continuum solutions of thermochemically relaxing gas mixtures and with experimental results from the high speed flow about a circular cylinder. Furthermore, the exchange of energy between translational and internal modes was modeled. New models were developed, combining greater efficiency than earlier phenomenological models. They offer greater physical detail by addressing vibrational energy as a discretely distributed quantity.

Mcdonald, Jeffrey D.↗

Spinon excitations and spin correlations in the one-dimensional quantum magnet β -VOSO 4 probed by Raman spectroscopy

Fractionalized excitations such as spinons and anyons have emerged as a central theme in condensed matter physics with broad implications for superconductivity, quantum statistics, and quantum computation. The nearly ideal one-dimensional S = 1/2 system β-VOSO 4 without long-range order down to 85 mK provides a promising platform to experimentally explore such fractionalized excitations. Here, we employ Raman spectroscopy to probe magnetic excitations and the evolution of spin correlations in β-VOSO 4 . Spinon signatures are found along the chain direction, evidenced by a broad, gapless scattering continuum at low temperatures. The temperature dependence of the spinon spectral weight aligns considerably with numerical density matrix renormalization group calculations. By comparing the experimental spinon spectral weight with calculated results and evaluating the associated quantum Fisher information (QFI) therefrom, we observe a steep increase in QFI upon cooling, indicating rapidly growing correlation lengths. Our study showcases QFI as a probe of spin correlations in quantum magnets.

Wulferding, Dirk [Sejong Univ., Seoul (Korea, Repu↗