Search NASA⌕ Search

SEARCH · Search NASA

Results for “Discontinuous fields”

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 379 records · Page 21

Estuarine density fronts and their effect on oil slicks

The author has identified the following significant results. Estuarine fronts represent regions of extremely high gradient or discontinuity in various parameters of physical interest, the most important being the water velocity and density fields. Aircraft and boats were combined to study the behavior of different types of fronts in Delaware Bay and their effect on pollutants in order to provide a basis for improving an oil drift and spreading model. Imagery from the LANDSAT satellites provided the most effective means of determining the location and extent of frontal systems over all portions of the tidal cycle. This data is being used to modify the oil drift and spreading model.

Klemas, V.↗

Remote sensing of estuarine fronts and their effects on pollutants

Estuarine fronts represent regions of extremely high gradient or discontinuity in various parameters of physical interest, the most important being the water velocity and density fields. Such fronts strongly influence pollutant dispersion by capturing oil slicks and other pollutants concentrated in surface films and drawing them down into the water column. Aircraft and boats were combined to study the behavior of different types of fronts in Delaware Bay and their effect on pollutants in order to provide a basis for improving an oil drift and spreading model. Imagery from the Landsat satellites provided the most effective means of determining the location and extent of frontal systems over all portions of the tidal cycle. This information is being used to modify the oil drift and spreading model.

Klemas, V.↗

A study of density fronts and their effects on coastal pollutants

Density fronts represent regions of extremely high gradient or discontinuity in various parameters of physical interest, the most important being the water velocity and density fields. Such fronts strongly influence pollutant dispersion, by capturing oil slicks and other pollutants concentrated in surface films and drawing them down into the water column. Satellites, aircraft and boats were used to study the behavior of different types of fronts in Delaware Bay and their effect on pollutants in order to provide a basis for improving an oil drift and spreading model. Landsat satellites provided the most effective means of determining the location and extent of frontal systems over all portions of the tidal cycle. Satellite observations of flood-associated fronts on the New Jersey side of the Bay and ebb-associated fronts on the Delaware side agreed with boat measurements and model predictions.

Klemas, V.↗

Numerical simulation of viscous-inviscid interactions on indented nose tips

An implicit numerical algorithm to solve the unsteady thin-layer Navier-Stokes equations in a strong conservative form has been used to compute the viscous flow over indented nose tips placed in a supersonic free stream. Numerical solutions are presented for axisymmetric and three-dimensional indented configurations for laminar flow conditions. Results demonstrate the capability of the present numerical procedure to predict flow fields that contain strong viscous-inviscid interactions, including boundary-layer separation, reattachment, and embedded discontinuities. Comparisons with available experimental data for the shock shape and surface pressure are also given.

Rizk, Y. M.↗

Structure and length scales of rotational discontinuities

Results are presented of a series of 1D hybrid (kinetic ions and fluid electrons) simulations, using a nondynamic method to form rotational discontinuities (RDs) at an angle of 60 deg between the normal direction and the upstream magnetic field. Ion kinetic effects are found to be important with regard to the stability of the current layer, its thickness and scaling, and its dependence on initial conditions. A detailed examination of the downstream wave train and gradient scale size showed that RDs with alpha1 = +/-180 deg are stable within the framework of the simulation. With proper initialization, electron-sense rotations of -270 deg were found that do not disintegrate immediately but rather appear to consist of a superposition of a solitary wave with -360-deg rotation and a 90-deg ion-sense rotation.

Krauss-Varban, D.↗

Anomalous resistivity and the evolution of magnetic field topology

This paper explores the topological restructuring of a force-free magnetic field caused by the hypothetical sudden onset of a localized region of strong anomalous resistivity. It is shown that the topological complexity increases, with the primitive planar force-free field with straight field lines developing field lines that wrap half a turn around each other, evidently providing a surface of tangential discontinuity in the wraparound region. It is suggested that the topological restructuring contributes to the complexity of the geomagnetic substorm, the aurora, and perhaps some of the flare activity on the sun, or other star, and the Galactic halo.

Parker, E. N.↗

Sound diffraction at wall impedance discontinuities in a circular cylinder - Investigated using Wiener-Hopf technique

The results of ground observations as well as high resolution rocket electric field and particle observations during a breakup event of an intense magnetospheric substorm over northern Scandinavia are discussed. In particular, the characteristics of the substorm-associated electric field, ionospheric currents, and power dissipation during a time period about 15 minutes after substorm onset are addressed. A comparison of the observations with those of a pre-breakup event earlier in the day (Marklund et al., 1982) showed that the ionospheric substorm-related electric field could be split up into two parts: (1) an ambient LT-dependent field, probably of magnetospheric origin; and (2) a small-scale electric field associated with the bright auroral structures, which is superimposed on the LT-dependent field. The consequences for the location of the ionospheric currents and the Joule energy dissipation relative to the auroral forms are discussed. Previously announced in STAR as N83-23117

Cho, Y.-C.↗

Correlated particle and magnetic field observations of a large-scale magnetic loop structure behind an interplanetary shock

From a survey of observations on ISEE-3, an example of correlated particle and magnetic field observations of a large-scale magnetic loop structure is presented. Bidirectional proton fluxes were observed for a period of 40 hours in the energy range 35-1600 keV approximately 12 hours after the passage of the interplanetary shock of December 11, 1980, and directly after the passage of a discontinuity. For each of the eight logarithmically spaced energy channels, a three-dimensional anisotropy analysis reveals streaming along both directions of the magnetic field. The magnetic field rotated slowly but steadily through approximately 180 deg during this same 40-hour period; this is consistent with the existence of a large-scale loop with extent greater than 0.5 AU. The observations suggest that the particles are being injected into the loop sunward of the spacecraft; they appear as bidirectional fluxes in the outermost regions of the loop arising from a combination of focusing and near scatter-free transport.

Sanderson, T. R.↗

Electromagnetic scattering from a class of open-ended waveguide discontinuities

A relatively simple high frequency analysis of electromagnetic scattering from a class of open-ended waveguide discontinuites was developed. The waveguides are composed of perfectly-conducting sections in which the electromagnetic field can be written as the sum of waveguide modes. Junctions are formed at the open end and also within interior regions where different sections are joined. The reflection and transmission properties of each junction are described in terms of a scattering matrix which is determined by combining the modal ray picture with high frequency techniques such as the Geometrical Theory of Diffraction (GTD), the Equivalent Current Method (ECM), and modifications of the Physical Theory of Diffraction (PTD). A new set of equivalent circuits are employed in this ECM analysis which leads to a simple treatment of many types of junction discontinuities. Also, a new procedure is presented to improve the efficiency of the aperture integration at the open end which is required in the PTD procedure for finding the fields radiated from (or coupled to) the open end. Once the scattering matrices are determined, they are then combined using a self-consistent multiple scattering method to obtain the total scattered fields.

Altintas, A.↗

Angular-spatial hp -adaptivity for radiative transfer with discontinuous Galerkin spectral element methods

Radiative transfer is important for many science and engineering applications, and numerical simulations of radiative transfer can be challenging. For instance, the radiation field is seven-dimensional – three spatial, two angular, one wavelength, and one temporal – and often features steep gradients. Therefore, memory usage is a key issue. To reduce memory, some past work has investigated the use of adaptive mesh refinement (AMR), typically for either the spatial or angular coordinate, and typically for only h -adaptivity. Here, we propose the use of AMR for the spatial and angular coordinates together, and the use of h - and p -adaptivity together as hp -AMR for the potential for further memory savings. We implemented the proposed method for several test cases in two spatial and one angular dimension, with the discontinuous Galerkin spectral element method. These test cases featured highly anisotropic angular radiation, with or without steep spatial gradients. Our primary findings from these test cases were: (1) Angular hp -adaptivity can deliver the radiation solution with the same accuracy as, and with much less computational memory than, uniform angular h - or p -refinements, or angular h -adaptivity alone. This is most obvious when the incoming radiation is highly anisotropic, in which case the savings can be orders of magnitude. (2) Full spatial-angular hp -adaptivity is more efficient in solution representation, compared to solely spatial or solely angular -adaptivity. This is most evident when steep gradients are present in both the spatial and angular distribution. These results suggest that adaptive spatial- hp angular-refinement may perform well in large-scale seven-dimensional applications.

Adaptive refinement↗

A phase-field diffraction model for thermo-hydro-mechanical propagating fractures

This paper introduces a novel diffraction based thermo-hydraulic–mechanical (THM) model for fracture propagation using a phase-field fracture (PFF) approach. The key innovation of the THM-PFF model lies in its integrated treatment of four solution variables—displacements, phase-field, pressure, and temperature—each governed by a combination of conservation of momentum (mechanics problem), a variational inequality (constrained minimization problem), mass conservation (pressure problem), and energy conservation (temperature problem). This leads to a new formulation of a coupled variational inequality system. A major advancement is the development of an extended fixed-stress algorithm, where displacements, phase-field, pressures, and temperatures are solved in a staggered sequence. An important aspect of this work is the global coupling of pressures and temperatures across the domain using diffraction systems, with diffraction coefficients defined by material parameters weighted by the diffusive phase-field variable. To ensure robust local mass conservation, we employ enriched Galerkin finite elements (EG) for both pressure and temperature diffraction equations. By enriching the continuous Galerkin basis functions with discontinuous piecewise constants, EG accurately represents solution and parameter discontinuities while preserving local mass and energy conservation—crucial aspects for THM problems and realistic behavior. Moreover, the use of a predictor–corrector local mesh adaptivity scheme is employed, allowing the model to handle small phase-field length-scale parameters while maintaining high numerical accuracy and reasonable computational cost. Furthermore, these new model and algorithmic developments represent significant advances in the field and have been substantiated through rigorous numerical tests.

Diffraction systems↗

Humate Amendment Injection Viability Testing – Supplementary Batch Mixture and Soil Column Studies

Groundwater in the Lost Lake Aquifer Zone (LLAZ) in the Southern Sector of the M-Area Hazardous Waste management Facility (HWMF) is contaminated with chlorinated ethenes, including trichloroethylene (TCE) and tetrachloroethylene (PCE). Treatment of contaminated groundwater with humic acid is being evaluated as a potential corrective action for these volatile organic compounds (VOCs) in the LLAZ in Southern Sector (SRNS, 2019). Pilot scale injections of Huma-K brand humic acid for groundwater treatment were previously performed in M-Area between 2017 and 2020 (Amidon, 2023). Groundwater was extracted from the Lost Lake Aquifer Zone (LLAZ), then a solution of humate was mixed with recovered groundwater and subsequently re-injected into the same aquifer unit. Operations of the historical humate pilot testing were discontinued in 2020 due to COVID-19 restrictions, low injection rates, and recirculation fluid spillage. No further field scale or laboratory research to investigate humate amended groundwater injection was performed at the SRS until 2024. At the request of Area Completion Projects (ACP), the Savannah River National Laboratory (SRNL) performed supplementary testing to help determine the viability of further injection of humate as a remedial option for VOCs within the LLAZ in the Southern Sector of the M-Area HWMF.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

On large poleward-directed electric fields at sub-auroral latitudes

Intense poleward-directed electric fields were observed by OGO-6 in the premagnetic midnight F-region, close to the Harang discontinuity and at latitudes of less than 60 deg invariant. These results support the recent measurements of Smiddy et al. (1977). Both sets of data are consistent with a compression of equipotentials on the low latitude edge of the premidnight convection cell, caused by an impediment of the substorm-related southwestward expansion of the convection pattern.

Maynard, N. C.↗

Numerical solution of two-dimensional turbulent blunt body flows with an impinging shock

An implicit finite-difference method has been developed to compute two-dimensional, turbulent, blunt body flows with an impinging shock wave. The full time-averaged Navier-Stokes equations are solved with algebraic eddy viscosity and turbulent Prandtl number models employed for shear stress and heat flux. The irregular-shaped bow shock is treated as a discontinuity across which the Rankine-Hugoniot equations are applied. A Type III turbulent shock interference flow field has been computed and the numerical results compare favorably with existing experimental data. In addition, comparisons are made between the present implicit code and a previous explicit code.

Tannehill, J. C.↗

The optical analogy for vector fields

This paper develops the optical analogy for a general vector field. The optical analogy allows the examination of certain aspects of a vector field that are not otherwise readily accessible. In particular, in the cases of a stationary Eulerian flow v of an ideal fluid and a magnetostatic field B, the vectors v and B have surface loci in common with their curls. The intrinsic discontinuities around local maxima in absolute values of v and B take the form of vortex sheets and current sheets, respectively, the former playing a fundamental role in the development of hydrodyamic turbulence and the latter playing a major role in heating the X-ray coronas of stars and galaxies.

Parker, E. N.↗

The Analysis and Design of Low Boom Configurations Using CFD and Numerical Optimization Techniques

The use of computational fluid dynamics (CFD) for the analysis of sonic booms generated by aircraft has been shown to increase the accuracy and reliability of predictions. CFD takes into account important three-dimensional and nonlinear effects that are generally neglected by modified linear theory (MLT) methods. Up to the present time, CFD methods have been primarily used for analysis or prediction. Some investigators have used CFD to impact the design of low boom configurations using trial and error methods. One investigator developed a hybrid design method using a combination of Modified Linear Theory (e.g. F-functions) and CFD to provide equivalent area due to lift driven by a numerical optimizer to redesign or modify an existing configuration to achieve a shaped sonic boom signature. A three-dimensional design methodology has not yet been developed that completely uses nonlinear methods or CFD. Constrained numerical optimization techniques have existed for some time. Many of these methods use gradients to search for the minimum of a specified objective function subject to a variety of design variable bounds, linear and nonlinear constraints. Gradient based design optimization methods require the determination of the objective function gradients with respect to each of the design variables. These optimization methods are efficient and work well if the gradients can be obtained analytically. If analytical gradients are not available, the objective gradients or derivatives with respect to the design variables must be obtained numerically. To obtain numerical gradients, say, for 10 design variables, might require anywhere from 10 to 20 objective function evaluations. Typically, 5-10 global iterations of the optimizer are required to minimize the objective function. In terms of using CFD as a design optimization tool, the numerical evaluation of gradients can require anywhere from 100 to 200 CFD computations per design for only 10 design variables. If one CFD computation requires an hour of computational time on a Cray computer, one can see that the use of constrained numerical optimization quickly becomes impractical.Hence, in order to practically couple a numerical design optimization technique with a CFD method, the CFD method must be extremely efficient with running times on the order of only minutes. The CFD Euler code developed under NASA sponsorship and referred to as MIM3D-SB for the most part fulfills these efficiency requirements. Analysis of wing- body configurations can be computed in a matter of a few minutes. The present study will concentrate on the feasibility of the use of this CFD code in conjunction with a numerical design optimization technique for the sonic boom reduction of candidate HSCT configurations. A preliminary supersonic aircraft design system has been established that utilizes the numerical design optimization code NPSOL developed at Stanford University coupled with the supersonic NUM3D-SB CFD code. Many questions still need to be answered in regard to using CFD and numerical optimizers as design tools. There are difficulties related to both the CFD codes and the numerical optimizers. Numerical optimizers can converge to a local minima rather than a global minima. This behavior is largely a function of the initial guess in the design space. The optimizer also is searching for a minimum of the function in terms of its derivative without any regard to the actual function value. Numerically (i.e. CFD) determined gradients can also generate spurious numerical local minima. In addition, for the sonic boom problem, grid fineness will also determine the accuracy of the final design solution. Design optimization methods work well on problems defined by continuous objective functions. The sonic boom signature design problem is not necessarily defined by a continuous objective function. The signature can have a variety of shapes; i.e. from N-wave to multiple shocks. The far-field or ground signature may not transition continuously from one shape to another and hence, may exhibit discontinuous behavior. This is also a source of difficulty in using design optimization methods.In the following sections, several low boom and one reference aircraft configuration will be analyzed to predict their sonic boom signature characteristics. Modifications to some of these designs will also be presented to demonstrate the feasibility of using CFD as a design tool and to demonstrate the feasibility of designing shaped sonic boom signatures. Design modifications to some configurations will be presented to demonstrate the feasibility of achieving shaped signatures with reduced levels and not necessarily to represent realistic or aerodynamically efficient design modifications. Fuselage volume or camber are used as design variables in order to have a minimal effect on the primary wing aerodynamics. The paper will also seek to demonstrate whether a hybrid or ramped signature is feasible to achieve. For the low-boom configurations, the CFD predicted signatures will be compared qualitatively to their MLT design signatures.

Siclari, Michael J.↗

Simulations of Sea-Ice Dynamics Using the Material-Point Method

In recent years, the availability of large volumes of recorded ice motion derived from high-resolution SAR data has provided an amazingly detailed look at the deformation of the ice cover. The deformation is dominated by the appearance of linear kinematic features that have been associated with the presence of leads. These remarkable data put us in a position to begin detailed evaluation of current coupled mechanical and thermodynamic models of sea ice. This presentation will describe the material point method (MPM) for solving these model equations. MPM is a numerical method for continuum mechanics that combines the best aspects of Lagrangian and Eulerian discretizations. The material points provide a Lagrangian description of the ice that models convection naturally. Thus, properties such as ice thickness and compactness are computed in a Lagrangian frame and do not suffer from errors associated with Eulerian advection schemes, such as artificial diffusion, dispersion, or oscillations near discontinuities. This desirable property is illustrated by solving transport of ice in uniform, rotational and convergent velocity fields. Moreover, the ice geometry is represented by unconnected material points rather than a grid. This representation facilitates modeling the large deformations observed in the Arctic, as well as localized deformation along leads, and admits a sharp representation of the ice edge. MPM also easily allows the use of any ice constitutive model. The versatility of MPM is demonstrated by using two constitutive models for simulations of wind-driven ice. The first model is a standard viscous-plastic model with two thickness categories. The MPM solution to the viscous-plastic model agrees with previously published results using finite elements. The second model is a new elastic-decohesive model that explicitly represents leads. The model includes a mechanism to initiate leads, and to predict their orientation and width. The elastic-decohesion model can provide similar overall deformation as the viscous-plastic model; however, explicit regions of opening and shear are predicted. Furthermore, the efficiency of MPM with the elastic-decohesive model is competitive with the current best methods for sea ice dynamics. Simulations will also be presented for an area of the Beaufort Sea, where predictions can be validated against satellite observations of the Arctic.

Sulsky, D.↗

Don't go with the Flow: An Invitation to Magnetosheath and Foreshock Studies

This talk reviews the predictions of gasdynamic, magnetohydrodynamic, and kinetic models for the magnetosheath and foreshock and compares these predictions with observations by the recent Cluster and THEMIS missions. Topics of interest include: the depletion layer, dawn/dusk asymmetries, the transmission of solar wind discontinuities, the formation of hot flow anomalies and cavities in the foreshock, and flows accelerated by field-line tension. We conclude by discussing opportunities for magnetosheath imaging.

Sibeck, D. G.↗