Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Lagrangian modeling”

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 721 records · Page 40

Froggatt-Nielsen meets the SMEFT

We study the matching of Froggatt-Nielsen theories of flavour onto the Standard Model Effective Field Theory (SMEFT), upon integrating out a heavy Beyond-the-Standard-Model (BSM) scalar ‘flavon’ whose vacuum expectation value breaks an Abelian flavour symmetry at energies Λ FN well above the electroweak scale, Λ FN > Λ SM . We include matching contributions to the infrared d SM = 6 (Warsaw basis) SMEFT sourced from ultraviolet contact terms suppressed up to order 1/${\Lambda}_{\textrm{UV}}^2$ in the Froggatt-Nielsen Lagrangian, where Λ UV > Λ FN is an arbitrary ultraviolet scale where further unspecified BSM particles are dynamical. This includes tree-level (one-loop) ultraviolet diagrams with d FN = 6 (5) effective vertices. We first do so with a toy model, but then generalize our findings to arbitrary Frogatt-Nielsen charges. Our results indicate a rich and non-trivial signature of Froggatt-Nielsen theories on the (otherwise) model-independent operators of the SMEFT, and we briefly speculate on extending our analysis to broader classes of BSM flavour models, e.g. non-Abelian and/or gauged theories. We thus take an important step towards determining how to use rapidly developing theoretical and experimental SMEFT technologies to gain unambiguous insight into the SM’s longstanding fermion flavour puzzle.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A scalable planning framework of energy storage systems under frequency dynamics constraints

As the penetration of renewables increases in power systems, the declining system inertia can cause frequency stability issues. Battery energy storage systems (BESSs) respond fast and therefore can relieve the low inertia difficulty but need to be appropriately sized considering the associated cost. This paper presents a novel stochastic optimization model for economically planning BESS capacity while considering the spatial–temporal correlation of wind generation and generator outages under frequency stability constraints, which include the rate-of-change of frequency (RoCoF), frequency nadir (FN), and quasi-steady-state (QSS) frequency. A set of new FN constraints that can be easily linearized is developed. To account for renewable uncertainties, a realistic uncertainty modeling approach, Random Field, is adopted to generate wind generation scenarios by considering both spatial and temporal evolutions of wind speed profiles. The ESS sizing is formulated as a mixed-integer linear programming problem and solved by using a scalable decomposition-and-coordination approach, Surrogate Absolute Value Lagrangian Relaxation (SAVLR). To further improve the scalability and reduce computational burdens, a rolling-horizon-based update is developed and incorporated into SAVLR for providing a practical solution to the long-term planning of very large-scale power systems. Finally, a modified IEEE 118-bus system and the Polish system are used to validate the effectiveness and scalability of the model and solution methodology.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Quadrature Moments Method for the Simulation of Turbulent Reactive Flows

A sub-filter model for reactive flows, namely the DQMOM model, was formulated for Large Eddy Simulation (LES) using the filtered mass density function. Transport equations required to determine the location and size of the delta-peaks were then formulated for a 2-peak decomposition of the FDF. The DQMOM scheme was implemented in an existing structured-grid LES solver. Simulations of scalar shear layer using an experimental configuration showed that the first and second moments of both reactive and inert scalars are in good agreement with a conventional Lagrangian scheme that evolves the same FDF. Comparisons with LES simulations performed using laminar chemistry assumption for the reactive scalar show that the new method provides vast improvements at minimal computational cost. Currently, the DQMOM model is being implemented for use with the progress variable/mixture fraction model of Pierce. Comparisons with experimental results and LES simulations using a single-environment for the progress-variable are planned. Future studies will aim at understanding the effect of increase in environments on predictions.

Raman, Venkatramanan↗

Detailed chemistry modeling of rotating detonations with dilute n -heptane sprays and preheated air

Utilization of liquid fuels is crucial to enabling commercialization of rotating detonation engines (RDE) in the near future. In this study, Eulerian-Lagrangian simulations are conducted for rotating detonative combustion with dilute n-heptane sprays and preheated air. Two-dimensional flattened configuration is used and a skeletal chemical mechanism with 44 species and 112 elementary reactions for n-heptane combustion is adopted. The flow structure, droplet distribution, and thermochemical parameters in the refill zone are first analyzed. It is shown that the mixture in the refill zone is heterogeneous, including evaporating droplets, vapor, and air. When the total temperature is below 950 K, the average equivalence ratio increases with the total temperature. When it is higher than 950 K, the average equivalence ratio is almost constant. Subsequently, the chemical explosive mode analysis is applied to identify the controlling reactions and dominant combustion modes in the fuel refill zone and reaction fronts. Results demonstrate that the initiation reaction (R104: n-C 7 H 16 + O 2 2-C 7 H 15 + HO 2 ) and low-temperature reaction (R107: RO 2 R'O 2 H) are dominant in the upstream and downstream of the refill zone, respectively. The intermediate species from low-temperature chemistry, R'O 2 H, is found to be important for the chemical explosive mode in the undetonated mixture. The influence of species diffusion and dispersed droplets is further analyzed. Results show that vapor autoignition facilitated by droplet evaporation occurs in the refill zone. Finally, the effects of the air total temperature on the detonation propagation speed and RDE propulsion performance are investigated. Further, it is found that the detonation propagation speed and specific impulse increase with air total temperature. The total pressure ratio first increases and then decreases as the air total temperature increases. Moreover, when the total temperature of the preheated air is above 1,300 K, the effects of low-temperature chemistry are negligible.

33 ADVANCED PROPULSION SYSTEMS↗

A process to verify numerical models for seismic fluid-structure interaction in advanced reactor vessels

Seismic design and qualification of a liquid-filled advanced nuclear reactor will have to account for fluid-structure interaction (FSI). Interaction between the tank, internal components, and contained liquid will rely on analysis of numerical models that must be verified and validated. Here this study demonstrates a verification process for models of a base-supported cylindrical tank by comparing numerical predictions and analytical solutions. The numerical models are consistent with the assumptions made to derive analytical solutions, namely, either a rigid or a linear elastic tank, ideal fluid, and small-amplitude, unidirectional, horizontal inputs. One software platform is used to illustrate the process. Seismic FSI analysis is performed using the Arbitrary Lagrangian-Eulerian (ALE) and Incompressible Computational Fluid Dynamics (ICFD) solvers in LS-DYNA. Reported responses are those used for design, including hydrodynamic pressures on the tank wall, shear forces and moments at the tank base, and wave heights of the contained liquid. The accuracy of the numerical results is discussed. The numerical models are verified for calculating the pressures on the tank wall and reactions at its base. Accurate simulation of wave action is challenging for both solvers. Recommendations for modeling, code development, and steps for verification are provided. Although focused on reactor vessels and one software platform, the verification process described herein is broadly applicable to liquid-filled vessels and other finite element codes.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Validity of standard model EFT studies of V H and V V production at NLO

The production of $W^±H, ZH, W^+W^-$, and $W^±Z$ pairs probes non-Standard-Model interactions of quarks, gauge bosons, and the Higgs boson. New effects can be parameterized in terms of an effective field theory (EFT) where the Lagrangian is expanded in terms of higher-dimension operators suppressed by increasing powers of a high scale Λ. We examine the importance of including next-to-leading-order QCD corrections in global fits to the coefficients of the EFT. The numerical implications on the fits due to different approaches to enforcing the validity of the EFT are quantified. We pay particular attention to the dependence of the fits on the ex-pansion in 1/Λ 2 since the differences between results calculated at $\mathcal{O}$(1/Λ 2 ) and $\mathcal{O}$(1/Λ 4 ) may give insight into the possible significance of dimension-8 effects.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

On the generation and evolution of internal gravity waves

The tidal generation and evolution of internal gravity waves is investigated experimentally and theoretically using a two-dimensional two-layer model. Time-dependent flow is created by moving a profile of maximum submerged depth 7.7 cm through a total stroke of 29 cm in water above a freon-kerosene mixture in an 8.6-m-long 30-cm-deep 20-cm-wide transparent channel, and the deformation of the fluid interface is recorded photographically. A theoretical model of the interface as a set of discrete vortices is constructed numerically; the rigid structures are represented by a source distribution; governing equations in Lagrangian form are obtained; and two integrodifferential equations relating baroclinic vorticity generation and source-density generation are derived. The experimental and computed results are shown in photographs and graphs, respectively, and found to be in good agreement at small Froude numbers. The reasons for small discrepancies in the position of the maximum interface displacement at large Froude numbers are examined.

Lansing, F. S.↗

Role of the Quasi-Biennial Oscillation in the Transport of Aerosols from the Tropical Stratospheric Reservoir to Midlatitudes

The temporal evolution of the stratospheric aerosol distribution in the tropical stratospheric reservoir after the eruption of Mount Pinatubo was observed from 1992 to 1995 by the HALOE instrument on the UARS satellite. Since the spatial gradient of aerosol loading is large at the boundaries of the tropical stratospheric reservoir due to the volcanic aerosols, the effect of the meridional circulation on the distribution is seen clearly. The mechanism for dispersal of aerosol in the lower stratosphere from the tropics into midlatitudes strongly depends on the phase of the equatorial zonal wind. The time-latitude cross sections of the normalized distribution of aerosol on isentropic surfaces are used to observe the equatorial variation as well as change in meridional dispersal during the quasi-biennial period. Observed tropical stratospheric winds are used with a simple analytical dynamical model to examine transport processes of tracers from the tropics during several phases of the quasi-biennial oscillation (QBO) from 1992 to 1995. The Lagrangian meridional circulation in the tropics is consistent with the vertical and meridional velocities correlated with the QBO in the zonal wind. We find that vertical motion plays a crucial role in vertical and subsequent meridional transport. The pattern of meridional divergence derived from the vertical velocity is closely related to the observed HALOE aerosol distributions and their temporal development in the equatorial region. The westerly (easterly) shear phase of the QBO is associated with sinking (rising) motions at the equator and subsequent poleward (equatorward) transport in the lower stratosphere.

Choi, Wookap↗

A combined Eulerian-volume of fraction-Lagrangian method for atomization simulation

The tracking of free surfaces between liquid and gas phases and analysis of the interfacial phenomena between the two during the atomization and breakup process of a liquid fuel jet is modeled. Numerical modeling of liquid-jet atomization requires the resolution of different conservation equations. Detailed formulation and validation are presented for the confined dam broken problem, the water surface problem, the single droplet problem, a jet breakup problem, and the liquid column instability problem.

Seung, S. P.↗

A Thermo‐Flow‐Mechanics‐Fracture Model Coupling a Phase‐Field Interface Approach and Thermo‐Fluid‐Structure Interaction

This work proposes a novel approach for coupling non-isothermal fluid dynamics with fracture mechanics to capture thermal effects within fluid-filled fractures accurately. This method addresses critical aspects of calculating fracture width in enhanced geothermal systems, where the temperature effects of fractures are crucial. The proposed algorithm features an iterative coupling between an interface-capturing phase-field fracture method and interface-tracking thermo-fluid-structure interaction using arbitrary Lagrangian–Eulerian coordinates. We use a phase-field approach to represent fractures and reconstruct the geometry to frame a thermo-fluid-structure interaction problem, resulting in pressure and temperature fields that drive fracture propagation. We developed a novel phase-field interface model accounting for thermal effects, enabling the coupling of quantities specific to the fluid-filled fracture with the phase-field model through the interface between the fracture and the intact solid domain. We provide several numerical examples to demonstrate the capabilities of the proposed algorithm. In particular, we analyze mesh convergence of our phase-field interface model, investigate the effects of temperature on crack width and volume in a static regime, and highlight the method's potential for modeling slowly propagating fractures.

fracture↗

Using XR for Improving Scientific Discovery With Numerical Weather Models

Earth science (ES) digital twins will help us understand the complex interactions and interrelationships that make up our Earth system and the impacts of earth science phenomena on it. Our work addresses two underdeveloped areas in current ES digital twin work: improving the understanding and interaction with ES model outputs by using Virtual and Mixed Reality (XR) tools and improving the non-intuitive mapping of continuous ES natural phenomena to gridded reference frames in current numerical models. Traditionally, scientists working on ES view and analyze the results of calculated or measured observables with static 1-dimensional (1D), 2D or 3D plots displayed on flat computer screens or paper. Using such limited mediums, it can be very difficult to identify, track and understand the evolution of key features due to poor viewing angles and the nature of flat computer screens. In addition, numerical models, such as the NASA Goddard Earth Observing System (GEOS) ES model, are almost exclusively formulated, visualized and analyzed in an Eulerian reference frame with fixed grid points in space and time. However, ES phenomena such as convective clouds, hurricanes and wildfire smoke plumes are visualized and analyzed in a Lagrangian reference frame: therefore it is often difficult and unnatural to understand these phenomena in relation to each other, visualized either in an Eulerian or Lagrangian context. In 3D visualizations, data generally takes one of three forms: gridded (e.g., voxelized) data, where space is divided into regions; point clouds, where data is represented as a set of points; and meshes, where objects are rendered as surfaces composed of small polygons (usually triangles). A gridded, Eulerian reference frame has been the default representation for the 2D visual analysis of atmospheric data in part because the numerical methods used to generate atmospheric model data in the first place use a gridded approach, with equations defining the relationships between the physical variables in each of a grid's cells across successive timesteps. In our work, we are particularly interested in data from GEOS. Another reason why gridded representations tend to be used for visualizing data from such models is because trajectories are difficult to interpret from representations on 2D surfaces, due to line-of-sight ambiguity. Instead of a fixed grid from GEOS, we embed a trajectory model to simulate particles' movement throughout a GEOS run. We then ingest these particle trajectories as animated point clouds with a NASA open source XR toolkit, the Mixed Reality Exploration Toolkit (MRET), and merge GEOS data with ES phenomena data onto one combined visualization that the user can intuitively interact with. Efficient rendering of arbitrarily large point clouds is an ongoing challenge being addressed by the computer science community, with the GPU-based optimizations and efficient GPU memory utilization a common theme of recent advances, especially for XR, where sustained high frame rate is mandatory to save the user from suffering due to simulation sickness. In this work, we describe and evaluate our progress in choosing and implementing appropriate methods for rendering arbitrarily large point clouds within MRET for XR. While tracking the XR headset enables the immersion of a user within a 3D scene of a data visualization, tracking of XR handheld controllers or user’s hands enables us to implement intuitive user interactions with the visualized datasets. Conventional tools require a user working with an ES visualization to conduct many interactions to commit their intended selections or manipulations with a visualized dataset; for example to specify a set of points in 3D space. Doing so in a 2D flat screen interface has traditionally required specifying a set of points in three distinct 2D coordinate systems (XY, XZ, and YZ), which is cumbersome. In other scientific domains, it has been shown that specifying or selecting a location or volume in XR using handheld controllers or tracked hands allows for greater speed and accuracy. We anticipate the same will hold true for atmospheric data, and we will share initial results of measuring the utility of such an interface. Notably, as the data being visualized is generated by GEOS as a prediction based on initial conditions, an intended application of our tool is to serve as part of an iterative feedback loop. Through XR, a scientist will review and manipulate a GEOS model run, modifying the conditions as needed to do subsequent runs of GEOS. Thereby, XR-based improvements to speed and accuracy of 3D tagging of points minimizes the effort required by both the scientist and the computer cluster conducting the necessary calculations.

Thomas Grubb↗

A dynamical basis for the parameterization of organized deep convection in large-scale numerical models

A hierarchy of steady, nonlinear, semianalytic models of different types of convection were produced. These provide a theoretical framework for determining cloud outflow fluxes of both dynamic and thermodynamic quantities, which can be used to formulate dynamical transports in parameterization schemes. This was achieved by exploiting certain Lagrangian conservation properties of steady flow, from which an equation for the vertical displacement of particles can be obtained and the outflow entropy, energy and momentum fluxes and the infow/outflow mass fluxes can be determined from solution to the equation. These fluxes are determined in terms of grid scale parameters such as convective available potential energy (CAPE), cloud layer shear, and horizontal pressure gradients. Five main types of system models are identified, respectively representing archtypes of convection in zero shear, large shear, midlatitude squall lines, tropical squall lines and cellular convection. The downdraught is an important aspect in the first four of these and the cloud scale transport of momentum is very distinctive.

Moncrieff, M. W.↗

Influence of the symmetry energy on the nuclear binding energies and the neutron drip line position

A clear connection can be established between properties of nuclear matter and finite-nuclei observables, such as the correlation between the slope of the symmetry energy and the dipole polarizability, or between compressibility and the isoscalar monopole giant resonance excitation energy. Establishing a connection between realistic atomic nuclei and an idealized infinite nuclear matter leads to a better understanding of underlying physical mechanisms that govern nuclear dynamics. In this work, we aim to study the dependence of the binding energies and related quantities (e.g., location of drip lines, the total number of bound even-even nuclei) on the symmetry energy S 2 (ρ). The properties of finite nuclei are calculated by employing the relativistic Hartree-Bogoliubov model, assuming even-even axial and reflection symmetric nuclei. Calculations are performed by employing two families of relativistic energy density functionals, based on different effective Lagrangians, constrained to a specific symmetry energy at the saturation density J within the interval of 30–36 MeV. Nuclear binding energies and related quantities of bound nuclei are calculated between 8 ≤ Z ≤ 104 from the two-proton to the two-neutron drip line. As the neutron drip line is approached, the interactions with stiffer J tend to predict more bound nuclei, resulting in a systematic shift of the two-neutron drip line towards more neutron-rich nuclei. Consequentially, a correlation between the number of bound nuclei N nucl and S 2 (ρ) is established for a set of functionals constrained using the similar optimization procedures. Furthermore the direction of the relationship between the number of bound nuclei and the symmetry energy highly depends on the density under consideration.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

User guide for the Geospatial Risk Analysis Tool (V.1.0.1)

Developed by the Pacific Northwest National Laboratory (PNNL), the Geospatial Risk Analysis Tool (GRAT) links to Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT), a complete system for computing simple air parcel trajectories, as well as complex transport, dispersion, chemical transformation, and deposition simulations. The HYSPLIT model was developed by the National Oceanic and Atmospheric Administration (NOAA) Air Resources Laboratory and the Australian Bureau of Meteorology Research Center in 1998. GRAT is used to batch process meteorological data into ARL format. While HYSPLIT GUI limits users to convert meteorological data up to 6 date time points per operation, GRAT enables users to convert a significantly large amount of data (For example, 10 years of meteorological data measured with an interval of 5 minutes) within one operation. Incorporated with the Potential source distribution function (PSDF), GRAT can also be used to determine the areas influenced by the emission of hazardous chemicals.

54 ENVIRONMENTAL SCIENCES↗

Quantifying vertical wind shear effects in shallow cumulus clouds over Amazonia

Abstract. This study analyses and quantifies the effects of vertical wind shear (VWS) on the properties of shallow cumulus cloud fields over Central Amazonia. We perform idealised simulations with high resolution (50 m horizontally and 20 m vertically) using the Dutch Atmospheric Large-Eddy Simulation (DALES) model, changing the initial conditions and large-scale forcing of VWS. The resulting cloud field is analysed by applying a cloud tracking algorithm to generate Lagrangian datasets of the life cycle of individual clouds as well as their time-varying core and margin dimensions. The reference run has no wind speed or directional shear and represents a typical day in the local dry season. Numerical experiments with moderate and high wind speed shear are simulated by adding linear increases in the wind speed of 1.2 and 2.4 m s−1 km−1, respectively. Three additional runs are made by adding 90∘ of wind rotation between the surface and the top of the domain (5 km) on top of the three wind speed shear conditions. We find that clouds developing in a sheared environment have horizontal equivalent diameter increased by up to 100 m on average, but the cloud depth is reduced. Our quantification shows that VWS tends to increase the size of the cloud cores but reduces their relative area, volume, and mass fractions compared to the overall cloud dimensions. The addition of 2.4 m s−1 km−1 of VWS decreases the relative core area by about 0.03 (about 10 % of the overall average) and its volume and mass ratios by about 0.05 (10 %–25 % in relative terms). Relevant for the cloud transport properties is that the updraught speed and the liquid water content are lower within the cores, and consequently so is the upward mass flux. All quantifications of mean cloud properties point to the inhibition of convective strength by VWS, therefore hampering the shallow-to-deep transition. However, open questions still remain given that the individually deepest clouds were simulated under high environmental shear, even though they occur in small numbers. This could indicate other indirect effects of VWS that have opposite effects on cloud development if found to be significant in the future.

54 ENVIRONMENTAL SCIENCES↗

Effect of steady state coning angle and damping on whirl flutter stability

The main object of this investigation is to find the effect of the steady state coning angle and the damping at the flapping hinge of the blades on the whirl flutter stability boundary and thus to determine the role they can play in narrowing down the gap between theory and experiment. The governing equations of motion, with these two parameters included are derived by the classical Lagrangian approach using quasi-steady blade element theory for aerodynamic forces. A linearized analysis of these equations is applied to two of the wind tunnel models of the previous investigations. The results indicate that these parameters have a marked effect on stability boundary and they may even change the mode of flutter from backward whirl to forward whirl.

Kaza, K. R. V.↗

Lagrangian methods in nonlinear plasma wave interaction

Analysis of nonlinear plasma wave interactions is usually very complicated, and simplifying mathematical approaches are highly desirable. The application of averaged-Lagrangian methods offers a considerable reduction in effort, with improved insight into synchronism and conservation (Manley-Rowe) relations. This chapter indicates how suitable Lagrangian densities have been defined, expanded, and manipulated to describe nonlinear wave-wave and wave-particle interactions in the microscopic, macroscopic and cold plasma models. Recently, further simplifications have been introduced by the use of techniques derived from Lie algebra. These and likely future developments are reviewed briefly.

Crawford, F. W.↗

Sensitivity of Age-of-Air Calculations to the Choice of Advection Scheme

The age of air has recently emerged as a diagnostic of atmospheric transport unaffected by chemical parameterizations, and the features in the age distributions computed in models have been interpreted in terms of the models' large-scale circulation field. This study shows, however, that in addition to the simulated large-scale circulation, three-dimensional age calculations can also be affected by the choice of advection scheme employed in solving the tracer continuity equation, Specifically, using the 3.0deg latitude X 3.6deg longitude and 40 vertical level version of the Geophysical Fluid Dynamics Laboratory SKYHI GCM and six online transport schemes ranging from Eulerian through semi-Lagrangian to fully Lagrangian, it will be demonstrated that the oldest ages are obtained using the nondiffusive centered-difference schemes while the youngest ages are computed with a semi-Lagrangian transport (SLT) scheme. The centered- difference schemes are capable of producing ages older than 10 years in the mesosphere, thus eliminating the "young bias" found in previous age-of-air calculations. At this stage, only limited intuitive explanations can be advanced for this sensitivity of age-of-air calculations to the choice of advection scheme, In particular, age distributions computed online with the National Center for Atmospheric Research Community Climate Model (MACCM3) using different varieties of the SLT scheme are substantially older than the SKYHI SLT distribution. The different varieties, including a noninterpolating-in-the-vertical version (which is essentially centered-difference in the vertical), also produce a narrower range of age distributions than the suite of advection schemes employed in the SKYHI model. While additional MACCM3 experiments with a wider range of schemes would be necessary to provide more definitive insights, the older and less variable MACCM3 age distributions can plausibly be interpreted as being due to the semi-implicit semi-Lagrangian dynamics employed in the MACCM3. This type of dynamical core (employed with a 60-min time step) is likely to reduce SLT's interpolation errors that are compounded by the short-term variability characteristic of the explicit centered-difference dynamics employed in the SKYHI model (time step of 3 min). In the extreme case of a very slowly varying circulation, the choice of advection scheme has no effect on two-dimensional (latitude-height) age-of-air calculations, owing to the smooth nature of the transport circulation in 2D models. These results suggest that nondiffusive schemes may be the preferred choice for multiyear simulations of tracers not overly sensitive to the requirement of monotonicity (this category includes many greenhouse gases). At the same time, age-of-air calculations offer a simple quantitative diagnostic of a scheme's long-term diffusive properties and may help in the evaluation of dynamical cores in multiyear integrations. On the other hand, the sensitivity of the computed ages to the model numerics calls for caution in using age of air as a diagnostic of a GCM's large-scale circulation field.

Eluszkiewicz, Janusz↗