Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation Cases”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Thermal mechanical assessment of a SiC-SiC-composite clad fuel pin concept in a light water reactor environment

Accident Tolerant Fuels (ATFs) are designed to increase coping time following an accident scenario while preserving or improving current steady state reactor operational performance. A potential ATF concept is SiC-SiC composite claddings. Fuel performance simulations were conducted on a SiC-SiC based cladding concept utilizing a multilayered approach for improved performance. This cladding concept referred to as the Duplex concept is a duplex structure composed of a monolithic SiC layer placed on the outside of a SiC-SiC composite. A liquid metal is added to fuel-cladding gap for improved heat dissipation from the fuel. The monolithic SiC layer is used to improve the coolant corrosion characteristics and protect the SiC-SiC composite layer from exposure to the coolant. The fuel performance code BISON was used to conduct fuel performance simulations on the cladding concepts. Comparisons are made with a current prototypic fuel rod design (UO 2 fuel enclosed in Zircaloy-4 cladding). Representative steady-state cases were considered for normal power and two cycle power histories. Additionally, a PCI ramp case was simulated to analyze potential anticipated operational occurrences. Transient response during a Loss of Coolant Accident and a Reactivity Initiated Accident were also simulated. This computational study demonstrated that for normal operating conditions, the SiC concept cladding performed as well as the baseline for the standard power cases evaluated. The ramping evaluations indicate potential fracturing of the SiC-SiC composite of the composite cladding compared to the Zircaloy-4 cladding due to the temperature gradient and the subsequent differential thermal conductivity degradation and swelling across the composite thickness. In conclusion, the rod fails early at low enthalpy for RIA but survives a LOCA with minimal material loss due to high temperature steam corrosion.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1981

Why use case studies rather than simulation-gaming techniques or library research?

Method which present a student with a more challenging and true to life situation of needing to conduct research in a problem solving context--and not thinking about organization of format until research and thinking are complete are investigated. Simulation-gaming techniques which attempt to teach initiative and creativity that library research are used for this purpose. However, it is shown case studies provide the greatest opportunities to engage the students in problem solving situations in which they develop skills as researchers and writers.

Mcdonald, S. W.↗
NASA NTRS2003

Simulation of the Onset of the Southeast Asian Monsoon During 1997 and 1998: The Impact of Surface Processes

The onset of the southeast Asian monsoon during 1997 and 1998 was simulated with a coupled mesoscale atmospheric model (MM5) and a detailed land surface model. The rainfall results from the simulations were compared with observed satellite data fiom the TRMM (Tropical Rainfall Measuring Mission) TMI (TRMM Microwave Imager) and GPCP (Global Precipitation Climatology Project). The simulation with the land surface model captured basic signatures of the monsoon onset processes and associated rainfall statistics. The sensitivity tests indicated that land surface processes had a greater impact on the simulated rainfall results than that of a small sea surface temperature change during the onset period. In both the 1997 and 1998 cases, the simulations were significantly improved by including the land surface processes. The results indicated that land surface processes played an important role in modifying the low-level wind field over two major branches of the circulation; the southwest low-level flow over the Indo- China peninsula and the northern cold front intrusion from southern China. The surface sensible and latent heat exchange between the land and atmosphere modified the lowlevel temperature distribution and gradient, and therefore the low-level. The more realistic forcing of the sensible and latent heat from the detailed land surface model improved the monsoon rainfall and associated wind simulation.

Wang, Yansen↗
NASA NTRS2004

Simulation of the Onset of the Southeast Asian Monsoon During 1997 and 1998: The Impact of Surface Processes

The onset of the southeast Asian monsoon during 1997 and 1998 was simulated with a coupled mesoscale atmospheric model (MM5) and a detailed land surface model. The rainfall results from the simulations were compared with observed satellite data from the TRMM (Tropical Rainfall Measuring Mission) TMI (TRMM Microwave Imager) and GPCP (Global Precipitation Climatology Project). The simulation with the land surface model captured basic signatures of the monsoon onset processes and associated rainfall statistics. The sensitivity tests indicated that land surface processes had a greater impact on the simulated rainfall results than that of a small sea surface temperature change during the onset period. In both the 1997 and 1998 cases, the simulations were significantly improved by including the land surface processes. The results indicated that land surface processes played an important role in modifying the low-level wind field over two major branches of the circulation; the southwest low-level flow over the Indo-China peninsula and the northern cold front intrusion from southern China. The surface sensible and latent heat exchange between the land and atmosphere modified the low-level temperature distribution and gradient, and therefore the low-level. The more realistic forcing of the sensible and latent heat from the detailed land surface model improved the monsoon rainfall and associated wind simulation. The model results will be compared to the simulation of the 6-7 May 2000 Missouri flash flood event. In addition, the impact of model initialization and land surface treatment on timing, intensity, and location of extreme precipitation will be examined.

Tao, W.-K.↗
NASA NTRS2004

Simulation of the Onset of the Southeast Asian Monsoon during 1997 and 1998: The Impact of Surface Processes

The onset of the southeast Asian monsoon during 1997 and 1998 was simulated with a coupled mesoscale atmospheric model (MM5) and a detailed land surface model. The rainfall results from the simulations were compared with observed satellite data from the TRMM (Tropical Rainfall Measuring Mission) TMI (TRMM Microwave Imager) and GPCP (Global Precipitation Climatology Project). The simulation with the land surface model captured basic signatures of the monsoon onset processes and associated rainfall statistics. The sensitivity tests indicated that land surface processes had a greater impact on the simulated rainfall results than that of a small sea surface temperature change during the onset period. In both the 1997 and 1998 cases, the simulations were significantly improved by including the land surface processes. The results indicated that land surface processes played an important role in modifying the low-level wind field over two major branches of the circulation; the southwest low-level flow over the Indo-China peninsula and the northern cold front intrusion from southern China. The surface sensible and latent heat exchange between the land and atmosphere modified the low-level temperature distribution and gradient, and therefore the low-level. The more realistic forcing of the sensible and latent heat from the detailed land surface model improved the monsoon rainfall and associated wind simulation. The model results will be compared to the simulation of the 6-7 May 2000 Missouri flash flood event. In addition, the impact of model initialization and land surface treatment on timing, intensity, and location of extreme precipitation will be examined.

Tao, W.-K.↗
NASA NTRS2007

Effects of Deep Convection on Atmospheric Chemistry

This presentation will trace the important research developments of the last 20+ years in defining the roles of deep convection in tropospheric chemistry. The role of deep convection in vertically redistributing trace gases was first verified through field experiments conducted in 1985. The consequences of deep convection have been noted in many other field programs conducted in subsequent years. Modeling efforts predicted that deep convection occurring over polluted continental regions would cause downstream enhancements in photochemical ozone production in the middle and upper troposphere due to the vertical redistribution of ozone precursors. Particularly large post-convective enhancements of ozone production were estimated for convection occurring over regions of pollution from biomass burning and urban areas. These estimates were verified by measurements taken downstream of biomass burning regions of South America. Models also indicate that convective transport of pristine marine boundary layer air causes decreases in ozone production rates in the upper troposphere and that convective downdrafts bring ozone into the boundary layer where it can be destroyed more rapidly. Additional consequences of deep convection are perturbation of photolysis rates, effective wet scavenging of soluble species, nucleation of new particles in convective outflow, and the potential fix stratosphere-troposphere exchange in thunderstorm anvils. The remainder of the talk will focus on production of NO by lightning, its subsequent transport within convective clouds . and its effects on downwind ozone production. Recent applications of cloud/chemistry model simulations combined with anvil NO and lightning flash observations in estimating NO Introduction per flash will be described. These cloud-resolving case-study simulations of convective transport and lightning NO production in different environments have yielded results which are directly applicable to the design of lightning parameterizations for global chemical transport models. The range of mean values (factor of 3) of NO production per flash (or per meter of lightning channel length) that have been deduced from the model will be shown and compared with values of production in the literature that have been deduced using other methods, Results show that on a per flash basis, IC flashes are nearly as productive of NO as CG flashes. When combined with the global flash rate of 44 flashes per second from NASA's Optical Transient Detector (OTD) measurements, these estimates and the results from other techniques yield global NO production rates of 2-9 TgN/year. Vertical profiles of lightning NOx mass at the end of the 3-D storm simulations have been summarized to yield suggested profiles for use in global models. Simulations of the photochemistry over the 24 hours following a storm have been performed to determine the additional ozone production which can be attributed to lightning NO.

Pickering, Kenneth E.↗
DOE OSTI2025

Surrogate-driven design optimization with uncertainty constraints in Monte Carlo simulations

In multi-objective design tasks, the computational cost increases rapidly when high-fidelity simulations are used to evaluate objective functions. Surrogate models help mitigate this cost by approximating the simulation output, simplifying the design process. However, under high uncertainty, surrogate models trained on noisy data can produce inaccurate predictions, as their performance depends heavily on the quality of training data. This study investigates the impact of data uncertainty on two multi-objective design problems modelled using Monte Carlo transport simulations: a neutron moderator and an ion-to-neutron converter. For each, a grid search was performed using five different tally uncertainty levels to generate training data for neural network surrogate models. These models were then optimized using NSGA-III. The recovered Pareto-fronts were analyzed across uncertainty levels: in the moderator problem, normalized hypervolume dropped from 0.886 at 1.0% uncertainty to 0.748 at 10% uncertainty, while in the converter problem it remained near 0.50 for all cases. Average simulation times were also compared to evaluate the trade-off between accuracy and computational cost. Results show that the influence of simulation uncertainty is strongly problem-dependent. In the neutron moderator case, higher uncertainties led to exaggerated objective sensitivities and distorted Pareto-fronts, reducing normalized hypervolume. In contrast, the ion-to-neutron converter task was less affected—low-fidelity simulations produced results similar to those from high-fidelity data. These findings suggest that a fixed-fidelity approach is not optimal. Surrogate models can recover the Pareto-front under noisy conditions, and multi-fidelity studies help identify suitable uncertainty levels for each problem to balance efficiency and accuracy.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2022

Specific Ion and Electric Field Controlled Diverse Ion Distribution and Electroosmotic Transport in a Polyelectrolyte Brush Grafted Nanochannel

Controlling ion distribution inside a charged nanochannel is central to using such channels in diverse applications. Here, we show the possibility of using a charged polyelectrolyte (PE) brush grafted nanochannel for triggering diverse nanoscopic ion distribution and nanofluidic electroosmotic transport by controlling the valence and size of the counterions (that screen the charges of the PE brushes) and the strength of an externally applied axial electric field. We atomistically simulate separate cases of fully charged Polyacrylic acid (PAA) brush functionalized nanochannels with Na + , Cs + , Ca 2+ , Ba 2+ , and Y 3+ counterions screening the PE charges. Four key findings emerge from our simulations. First, we find that the counterions with a greater valence and a smaller size prefer to remain localized inside the brush layer. Second, for the case where there is an added chloride salt with the same cation (as the screening counterions), there are more coions (Cl - ions) in the brush-free bulk than counterions (for counterions Na + , Ca 2+ , Ba 2+ , Y 3+ ): this is a manifestation of the overscreening (OS) of the PE brush layer. Contrastingly, the number of Cs + ions remain higher than the Cl - ions inside the brush-free bulk, ensuring that there is no OS effect for this case. Third, large applied electric field enables a few Na + , Cs + , and Ba 2+ counterions to leave the brush layer and to go to the bulk: this makes the OS of the PE brush layer disappear for the cases of PE brushes being screened by the Na+ and Ba 2+ ions. On the other hand, no such electric-field-mediated disappearance of OS is observed for the cases of Ca 2+ and Y 3+ screening counterions; we attribute this to the firm attachment of these counterions to the negatively charged monomers. Free energy associated to a counterion binding to a PE chain corroborates this diversity in the counterion-specific response to the applied electric field. Lastly, we demonstrate that such diverse ion distributions, along with specific electric-field-strength-dependent ion properties, lead to (1) EOS transport in nanochannels grafted with PAA brushes screened with Cs + ions to be always counterion dominated, (2) EOS transport in nanochannels grafted with PAA brushes screened with Ca 2+ and Y 3+ ions to be always coion dominated, and (3) EOS transport in nanochannels grafted with PAA brushes screened with Na + and Ba 2+ ions to be coion dominated for smaller electric fields and counterion dominated for larger electric fields.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

The Importance of Being Adaptable: An Exploration of the Power and Limitations of Domain Adaptation for Simulation-Based Inference with Galaxy Clusters

The application of deep machine learning methods in astronomy has exploded in the last decade, with new models showing remarkably improved performance on benchmark tasks. Not nearly enough attention is given to understanding the models' robustness, especially when the test data are systematically different from the training data, or "out of domain." Domain shift poses a significant challenge for simulation-based inference, where models are trained on simulated data but applied to real observational data. In this paper, we explore domain shift and test domain adaptation methods for a specific scientific case: simulation-based inference for estimating galaxy cluster masses from X-ray profiles. We build datasets to mimic simulation-based inference: a training set from the Magneticum simulation, a scatter-augmented training set to capture uncertainties in scaling relations, and a test set derived from the IllustrisTNG simulation. We demonstrate that the Test Set is out of domain in subtle ways that would be difficult to detect without careful analysis. We apply three deep learning methods: a standard neural network (NN), a neural network trained on the scatter-augmented input catalogs, and a Deep Reconstruction-Regression Network (DRRN), a semi-supervised deep model engineered to address domain shift. Although the NN improves results by 17% in the Training Data, it performs 40% worse on the out-of-domain Test Set. Surprisingly, the Scatter-Augmented Neural Network (SANN) performs similarly. While the DRRN is successful in mapping the training and Test Data onto the same latent space, it consistently underperforms compared to a straightforward Yx scaling relation. These results serve as a warning that simulation-based inference must be handled with extreme care, as subtle differences between training simulations and observational data can lead to unforeseen biases creeping into the results.

Ntampaka, Michelle [Baltimore, Space Telescope Sci↗
DOE OSTI2022

Modelica-based modeling and simulation of district cooling systems: A case study

While equation-based object-oriented modeling language Modelica can evaluate practical energy improvements for district cooling systems, few have adopted Modelica for this type of large-scale thermo-fluid system. Further, to our best knowledge, district cooling modeling studies have yet to include hydraulics in piping networks alongside plant models featuring realistic mechanical systems and controls. These are critical details to include when looking to make energy and control improvements in many physical system installations. To fill these gaps, this study released new open-source district cooling models at the Modelica Buildings Library and leveraged these models for a real-world case study at the University of Colorado Boulder. Here, the site includes six buildings connected to a central chiller plant featuring a waterside economizer. Several energy saving strategies are pursued based on the validated model, including control setpoint optimization, equipment modification, and pump setpoint adjustments. Results indicate that a combination of the studied measures can save the campus annually 84.6 MWh of energy, 8.9% of electricity costs, 58.0 metric tons of carbon dioxide emissions, while the waterside economizer cuts down chillers’ run times by 201 days/year, reducing maintenance costs and extending chiller life.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

FIDES-II/P2M Simulation Exercise on AN3 and AN10 bump tests: main results

This paper presents simulations of bump tests conducted as part of the P2M (Power to Melt and Maneuverability) project within the OECD/NEA FIDES-II framework (Framework for Irradiation Experiments). To prepare fuel performance codes for analyzing the P2M experiments, an international Simulation Exercise (SE) was initiated under the P2M project. The objective of this exercise was to model the AN3 and AN10 bump tests carried out in the 1980s at the DR3 reactor (3rd RISØ Fission Gas Project), during which instrumented fuel rodlets (equipped with thermocouples and pressure sensors) were subjected to prolonged power transients and power dips. The ongoing simulation exercise involves seven organizations from industry, research, and regulatory bodies, each using different fuel performance codes. This paper provides detailed descriptions of the test cases and simulation results from six fuel performance codes (ALCYONE, BISON, FAST, FEMAXI, FRAPCON, and TRANSURANUS), with a particular focus on fission gas release kinetics.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS1997

Molecular Dynamics Simulation of a Multi-Walled Carbon Nanotube Based Gear

We used molecular dynamics to investigate the properties of a multi-walled carbon nanotube based gear. Previous work computationally suggested that molecular gears fashioned from (14,0) single-walled carbon nanotubes operate well at 50-100 gigahertz. The gears were formed from nanotubes with teeth added via a benzyne reaction known to occur with C60. A modified, parallelized version of Brenner's potential was used to model interatomic forces within each molecule. A Leonard-Jones 6-12 potential was used for forces between molecules. The gear in this study was based on the smallest multi-walled nanotube supported by some experimental evidence. Each gear was a (52,0) nanotube surrounding a (37,10) nanotube with approximate 20.4 and 16,8 A radii respectively. These sizes were chosen to be consistent with inter-tube spacing observed by and were slightly larger than graphite inter-layer spacings. The benzyne teeth were attached via 2+4 cycloaddition to exterior of the (52,0) tube. 2+4 bonds were used rather than the 2+2 bonds observed by Hoke since 2+4 bonds are preferred by naphthalene and quantum calculations by Jaffe suggest that 2+4 bonds are preferred on carbon nanotubes of sufficient diameter. One gear was 'powered' by forcing the atoms near the end of the outside buckytube to rotate to simulate a motor. A second gear was allowed to rotate by keeping the atoms near the end of its outside buckytube on a cylinder. The ends of both gears were constrained to stay in an approximately constant position relative to each other, simulating a casing, to insure that the gear teeth meshed. The stiff meshing aromatic gear teeth transferred angular momentum from the powered gear to the driven gear. The simulation was performed in a vacuum and with a software thermostat. Preliminary results suggest that the powered gear had trouble turning the driven gear without slip. The larger radius and greater mass of these gears relative to the (14,0) gears previously studied requires a smaller rotation rate and multiple rows of teeth to avoid excessive force on the gear teeth resulting, in slip and failure of the driven gear to turn. We hope that studies such as these will eventually lead to synthesis of components that can be assembled into atomically precise fullerene machines. These machines, in turn, may someday be used in machine-phase fullerene materials with remarkable properties.

Han, Jie↗
NASA NTRS2005

Effect of Free Jet on Refraction and Noise

This article investigates the role of a free jet on the sound radiated from a jet. In particular, the role of an infinite wind tunnel, which simulates the forward flight condition, is compared to that of a finite wind tunnel. The second configuration is usually used in experiments, where the microphones are located in a static ambient medium far outside the free jet. To study the effect of the free jet on noise, both propagation and source strength need to be addressed. In this work, the exact Green's function in a locally parallel flow is derived for a simulated flight case. Numerical examples are presented that show a reduction in the magnitude of the Green's function in the aft arc and an increase in the forward arc for the simulated flight condition. The effect of finite wind tunnel on refraction is sensitive to the source location and is most pronounced in the aft arc. A Reynolds-averaged Navier-Stokes solution (RANS) yields the required mean flow and turbulence scales that are used in the jet mixing noise spectrum calculations. In addition to the sound/flow interaction, the separate effect of source strength and elongation of the noise-generating region of the jet in a forward flight is studied. Comparisons are made with experiments for the static and finite tunnel cases. Finally, the standard free-jet shear corrections that convert the finite wind tunnel measurements to an ideal wind tunnel arrangement are evaluated.

Khavaran, Abbas↗
DOE OSTI2026

Detector alignment for X-ray crystallography using Millepede-II

I describe a method for accurately refining the geometrical parameters of segmented X-ray area detectors on the basis of serial crystallography data, using 'Millepede' – an algorithm created for a very similar problem in high-energy physics. The Millepede method for serial crystallography builds on the approach of Brewster et al. [Acta Cryst. (2018), D74, 877–894], in which the detector parameters are refined simultaneously with the parameters for each individual crystal. This accounts for the mutual dependency between the parameters and thereby avoids the bias and slow convergence problems that have afflicted older approaches in which the deviations between observed and calculated Bragg peak positions were taken directly as the updates for the detector panel positions. The Millepede method uses the special structure of the least-squares normal equations to reduce them to a much smaller form that can be solved very quickly, even compared with the sparse matrix methods used previously. This makes it practical to refine the detector geometry frequently and thereby maintain accurate calibration without specialized alignment campaigns. Tilts of detector panels out of the plane can be reliably refined, as can the overall distance of the detector in the beam direction. With a simulated test case, the new method produced panel shifts within 7% of the correct values with only one iteration, and produced almost exactly correct shifts after a second iteration. A simulated out-of-plane panel rotation was correctly determined to within 0.001°. Applied to experimental data from an X-ray free-electron laser, the method increased the indexable fraction of frames from 30% to 91% in a single iteration, and to 96% after two further iterations. Computing the geometry updates on the basis of 2060 crystals took only 0.819 s on desktop computing hardware, including the time taken to read the required data from disk. The scaling was found to be very close to linear for up to 100 980 sets of crystal parameters, which took only 78.2 s to process under the same conditions. The method has been applied as part of a real-time feedback system at a synchrotron radiation beamline, in which an out-of-plane detector tilt of 0.04° was detected and corrected. Possible further applications are also described here.

Millepede-II↗
NASA NTRS1992

Large transient fault current test of an electrical roll ring

The space station uses precision rotary gimbals to provide for sun tracking of its photoelectric arrays. Electrical power, command signals and data are transferred across the gimbals by roll rings. Roll rings have been shown to be capable of highly efficient electrical transmission and long life, through tests conducted at the NASA Lewis Research Center and Honeywell's Satellite and Space Systems Division in Phoenix, AZ. Large potential fault currents inherent to the power system's DC distribution architecture, have brought about the need to evaluate the effects of large transient fault currents on roll rings. A test recently conducted at Lewis subjected a roll ring to a simulated worst case space station electrical fault. The system model used to obtain the fault profile is described, along with details of the reduced order circuit that was used to simulate the fault. Test results comparing roll ring performance before and after the fault are also presented.

Yenni, Edward J.↗
NASA NTRS1991

Large transient fault current test of an electrical roll ring

The Space Station Freedom uses precision rotary gimbals to provide for sun tracking of its photoelectric arrays. Electrical power, command signals, and data are transferred across the gimbals by roll rings. Roll rings have been shown to be capable of highly efficient electrical transmission and long life, through tests conducted at the NASA Lewis Research Center and Honeywell's Satellite and Space Systems Division in Phoenix, AZ. Large potential fault currents inherent to the power system's DC distribution architecture have brought about the need to evaluate the effects of large transient fault currents on roll rings. A test recently conducted at Lewis subjected a roll ring to a simulated worst case space station electrical fault. The system model used to obtain the fault profile is described, along with details of the reduced order circuit that was used to simulate the fault. Test results comparing roll ring performance before and after the fault are also presented.

Yenni, Edward J.↗
DOE OSTI2023

Unsaturated alluvium disposal modelling with improved geological realism

This report describes specific activities in the Fiscal Year (FY) 2023 associated with the Geologic Disposal Safety Assessment (GDSA) Repository Systems Analysis (RSA) work package funded by the Spent Fuel and Waste Science and Technology (SFWST) Campaign of the U.S. Department of Energy Office of Nuclear Energy (DOE-NE), Office of Spent Fuel and Waste Disposition (SFWD). The overall objective of the GDSA RSA work package is to develop generic deep geologic repository concepts and system performance assessment (PA) models in several potential host-rock environments, and to simulate and analyze these generic repository concepts and models using the GDSA Framework toolkit and other tools as needed. The specific GDSA goal addressed in this report is reference case development, simulation, and analysis for the unsaturated alluvium (UZ), one of the four potential host-rocks considered by the GDSA. Further, we aim to exercise the simulation tools and methodologies under development by GDSA for PA modelling.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1993

Numerical simulation of the subsolar magnetopause current layer in the sun-earth meridian plane

The formation and stability of the magnetopause current layer near the subsolar point in the sun-earth meridian plane are examined using a 2D electromagnetic particle simulation. For the case of zero IMF, the simulation results show that the current layer remains stable and is essentially the same as in the 1D simulation. The width of the current layer is given by the electron-ion hybrid gyroradius which is much smaller than the ion gyroradius. The current layer is found to remain stable for the northward IMF as well. As in the 1D simulation, the jump of the magnetic field at the current layer for the northward IMF remains small. For the southward IMF, collisionless magnetic reconnection is found to develop, leading to the formation of magnetic islands and density peaking within the current layer.

Okuda, H.↗