Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical continuation”

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 1,225 records · Page 68
All sourcesNASA NTRSDOE OSTI
NASA NTRS2001

Modelling of Trapped Radiation Near Jupiter

Energetic (62 to approx. 130 MeV) proton fluxes measured with the Galileo Probe inside Jupiter's main ring (radii 122,500 to 128,940 km) show a modest increase as Jupiter is approached. Solutions of a reduced (1-d) equatorial diffusion equation with-constant losses described by a lifetime tau match those data for tau greater than or equal to 10(exp 9) s for diffusion coefficient D(sub LL) = 10(exp -9) L(exp 4)/s, 10(exp -9) L(exp -3)/s, or 10(exp -10) L(exp 4)/s (this particle population may be undergoing nearly loss-free inward radial diffusion). Tau would increase as Jupiter is approached if its value were determined by pitch angle scattering due to EM wave-particle interactions. if the absolute amplitudes of the waves' magnetic fluctuations did not vary with radial distance. Exploration of numerical solutions of the diffusive transport equation for Jupiter's magnetospheric ions at die planetary ring and closer to Jupiter has been done. Explicit range-energy relationships for energy loss in ring matter (SiO2), modeled as a continuous disk, are incorporated. A spatial resolution of approx. 1000 km in the radial direction is used in the vicinity of the main ring; off-equatorial ion fluxes have been included (assuming the atmosphere is a perfect absorber). For protons the maximum effect of energy loss in the microscopic size range of ring particles is expected at energies approx. 0.1 Mev. Presumably there could be sufficiently large, unobserved bodies within the ring that would have to be modeled as discrete objects, as planetary satellites would be. Within limitations of the model and computational resources it's planned to characterize solutions for transport of magnetospheric ions past the planetary ring; the procedures that are used are fairly well known and may be applied generally to the simpler magnetospheric models.

Mihalov, John D.↗
NASA NTRS2002

Leveraging Improvements in Precipitation Measuring from GPM Mission to Achieve Prediction Improvements in Climate, Weather and Hydrometeorology

The main scientific goal of the GPM mission, currently planned for start in the 2007 time frame, is to investigate important scientific problems arising within the context of global and regional water cycles. These problems cut across a hierarchy of scales and include climate-water cycle interactions, techniques for improving weather and climate predictions, and better methods for combining observed precipitation with hydrometeorological prediction models for applications to hazardous flood-producing storms, seasonal flood/draught conditions, and fresh water resource assessments. The GPM mission will expand the scope of precipitation measurement through the use of a constellation of some 9 satellites, one of which will be an advanced TRMM-like "core" satellite carrying a dual-frequency Ku-Ka band precipitation radar and an advanced, multifrequency passive microwave radiometer with vertical-horizontal polarization discrimination. The other constellation members will include new dedicated satellites and co-existing Operational/research satellites carrying similar (but not identical) passive microwave radiometers. The goal of the constellation is to achieve approximately 3-hour sampling at any spot on the globe. The constellation's orbit architecture will consist of a mix of sun-synchronous and non-sun-synchronous satellites with the core satellite providing measurements of cloud-precipitation microphysical processes plus calibration-quality rainrate retrievals to be used with the other retrieval information to ensure bias-free constellation coverage. GPM is organized internationally, currently involving a partnership between NASA in the US and the National Space Development Agency in Japan. Additionally, the program is actively pursuing agreements with other international partners and domestic scientific agencies and institutions, as well as participation by individual scientists from academia, government, and the private sector to fulfill mission goals and to pave the way for what ultimately is expected to become an internationally-organized operational global precipitation observing system. Notably, the broad societal applications of GPM are reflected in the United Nation s identification of this mission as a foremost candidate for its Peaceful Uses of Space Program. In this presentation, an overview of the GPM mission design will be presented, followed by an explanation of its scientific agenda as an outgrowth of making improvements in rain retrieval accuracy, microphysics dexterity, sampling frequency, and global coverage. All of these improvements offer new means to observe variability in precipitation and water cycle fluxes and to achieve improved predictability of weather, climate, and hydrometeorology. Specifically, the scientific agenda of GPM has been designed to leverage the measurement improvements to improve prognostic model performance, particularly quantitative precipitation forecasting and its linked phenomena at short, intermediate, and extended time scales. The talk will address how GPM measurements will enable better detection of accelerations and decelerations in regional and global water cycle processes and their relationship to climate variability, better impacts of precipitation data assimilation on numerical weather prediction and global climate reanalysis, and better performance from basin scale hydrometeorological models for short and long term flood-drought forecasting and seasonal fresh water resource assessment. Improved hydrometeorological forecasting will be possible by using continuous global precipitation observations to obtain better closure in water budgets and to generate more realistic forcing of the models themselves to achieve more accurate estimates of interception, infiltration, evaporation/transpiration fluxes, storage, and runoff.

Smith, Eric A.↗
NASA NTRS2008

DMAC and NMP as Electrolyte Additives for Li-Ion Cells

Dimethyl acetamide (DMAC) and N-methyl pyrrolidinone (NMP) have been found to be useful as high-temperature-resilience-enhancing additives to a baseline electrolyte used in rechargeable lithium-ion electrochemical cells. The baseline electrolyte, which was previously formulated to improve low-temperature performance, comprises LiPF6 dissolved at a concentration of 1.0 M in a mixture comprising equal volume proportions of ethylene carbonate, diethyl carbonate, and dimethyl carbonate. This and other electrolytes comprising lithium salts dissolved in mixtures of esters (including alkyl carbonates) have been studied in continuing research directed toward extending the lower limits of operating temperatures and, more recently, enhancing the high-temperature resilience of such cells. This research at earlier stages, and the underlying physical and chemical principles, were reported in numerous previous NASA Tech Briefs articles. Although these electrolytes provide excellent performance at low temperatures (typically as low as -40 C), when the affected Li-ion cells are subjected to high temperatures during storage and cycling, there occur irreversible losses of capacity accompanied by power fade and deterioration of low-temperature performance. The term "high-temperature resilience" signifies, loosely, the ability of a cell to resist such deterioration, retaining as much as possible of its initial charge/discharge capacity during operation or during storage in the fully charged condition at high temperature. For the purposes of the present development, a temperature is considered to be high if it equals or exceeds the upper limit (typically, 30 C) of the operating-temperature range for which the cells in question are generally designed.

Smart, Marshall↗
NASA NTRS2010

Development and Breakdown of Goertler Vortices in High Speed Boundary Layers

The nonlinear development of G rtler instability over a concave surface gives rise to a highly distorted stationary flow in the boundary layer that has strong velocity gradients in both spanwise and wall-normal directions. This distorted flow is susceptible to strong, high frequency secondary instability that leads to the onset of transition. For high Mach number flows, the boundary layer is also subject to the second mode instability. The nonlinear development of G rtler vortices and the ensuing growth and breakdown of secondary instability, the G rtler vortex interactions with second mode instabilities as well as oblique second mode interactions are examined in the context of both internal and external hypersonic configurations using nonlinear parabolized stability equations, 2-D eigenvalue analysis and direct numerical simulation. For G rtler vortex development inside the Purdue Mach 6 Ludwieg tube wind tunnel, multiple families of unstable secondary eigenmodes are identified and their linear and nonlinear evolution is examined. The computation of secondary instability is continued past the onset of transition to elucidate the physical mechanisms underlying the laminar breakdown process. Nonlinear breakdown scenarios associated with transition over a Mach 6 compression cone configuration are also explored.

Li, Fei↗
NASA NTRS2011

Coupling Between Doppler Radar Signatures and Tornado Damage Tracks

On April 27, 2011, the southeastern United States was raked with several episodes of severe weather. Numerous tornadoes caused extensive damage, and tragically, the deaths of over 300 people. In Alabama alone, there were 61 confirmed tornados, 4 of them produced EF5 damage, and several were on the ground an hour or more with continuous damage tracks exceeding 80km. The use of Doppler radars covering the region provided reflectivity and velocity signatures that allowed forecasters to monitors the severe storms from beginning to end issuing hundreds of severe weather warnings throughout the day. Meteorologists from the the NWS performed extensive surveys to assess the intensity, duration, and ground track of tornadoes reported during the event. Survey activities included site visits to the affected locations, analysis of radar and satellite data, aerial surveys, and interviews with eyewitnesses. Satellite data from NASA's MODIS and ASTER instruments played a helpful role in determining the location of tornado damage paths and in the assessment. High resolution multispectral and temporal composites helped forecasters corroborate their damage assessments, determine starting and ending points for tornado touchdowns, and helped to provide forecasters with a better big-picture view of the damage region. The imagery also helped to separate damage from the April 27th tornados from severe weather that occurred earlier that month. In a post analysis of the outbreak, tornado damage path signatures observed in the NASA satellite data have been correlated to "debris ball" signatures in the NWS Doppler radars and a special ARMOR dual-polarization radar operated by the University of Alabama Huntsville during the event. The Doppler radar data indicates a circular enhanced reflectivity signal and rotational couplet in the radial velocity likely associated with the tornado that is spatially correlated with the damage tracks in the observed satellite data. An algorithm to detect and isolate the "debris ball" from precipitation signatures in the dual polarization radar data has been developed and verified using the NASA damage track data.

Jedlovec, Gary J.↗
DOE OSTI2024

MODELING, IMPLEMENTATION AND CONTROL OF A CVT BASED PTO FOR A SMALL SCALE MHK-TURBINE IN LOW FLOW SPEED OPERATION

Modeling, implementation, field testing and control of a power takeoff (PTO) device equipped with a ball-type continuously variable transmission (B-CVT) for a small marine hydrokinetic (MHK) turbine deployed from a floating unmanned autonomous mobile catamaran platform is described. The turbine is a partially submerged multi-blade undershot waterwheel (USWW). A validated numerical torque model for the MHK turbine has been derived and a speed controller has been developed, implemented and tested in the field. The dependance of the power generated as a function of number and submergence level of turbine blades has been investigated and the number of blades that maximizes power production is determined. Bench and field testing in support of characterizing the power conversion capabilities of MHK turbine and PTO are described. Detailed results of the final torque and power coefficient models, the controls architecture, and the MHK turbine performance with varying numbers of blades are provided.

Pimentel, Hugo↗
DOE OSTI2023

Blackout Diffusion: Generative Diffusion Models in Discrete-State Spaces

Typical generative diffusion models rely on a Gaussian diffusion process for training the backward transformations, which can then be used to generate samples from Gaussian noise. However, real world data often takes place in discrete-state spaces, including many scientific applications. Here, we develop a theoretical formulation for arbitrary discrete-state Markov processes in the forward diffusion process using exact (as opposed to variational) analysis. We relate the theory to the existing continuous-state Gaussian diffusion as well as other approaches to discrete diffusion, and identify the corresponding reverse-time stochastic process and score function in the continuous-time setting, and the reverse-time mapping in the discrete-time setting. As an example of this framework, we introduce “Blackout Diffusion”, which learns to produce samples from an empty image instead of from noise. Numerical experiments on the CIFAR-10, Binarized MNIST, and CelebA datasets confirm the feasibility of our approach. Generalizing from specific (Gaussian) forward processes to discrete-state processes without a variational approximation sheds light on how to interpret diffusion models, which we discuss.

Santos, Javier E.↗
NASA NTRS2024

Direct Numerical Simulation of Smooth-Body Flow Separation Around a Ramp

Spanwise-periodic computation of a turbulent flow past a two-dimensional smooth ramp geometry is performed in the form of a direct numerical simulation. The Reynolds number based on the ramp height is about 147,000. A straight section that precedes the smooth ramp allows the incoming turbulent boundary layer to grow under a weak favorable pressure gradient. The boundary layer introduced at the domain inlet has a momentum-thickness based Reynolds number of 2000. The turbulent boundary layer nearing the ramp first interacts with a relatively stronger favorable pressure gradient, followed by a strong adverse pressure gradient. Consequently, the boundary layer experiences a modest acceleration before decelerating and separating. Analysis of the data over this region hints at the formation of an internal layer beneath the accelerated boundary layer. The analysis also reveals that this internal layer forms the origin of the free shear layer that emerges in the deceleration region and separates. The streamwise extent of the separated region is comparable to the ramp length, while the viscous layer thickness near reattachment is about the same as the ramp height; hence, the boundary layer undergoing separation and subsequent reattachment in the present configuration experiences its thickness being amplified by about tenfold. The reattached flow continues to develop further under a diminishing pressure gradient in the recovery region in a similar fashion to a zero pressure gradient turbulent boundary layer.

Turbulence↗
DOE OSTI2024

Continuous-variable quantum computation of the O(3) model in 1+1 dimensions

We formulate the $O(3)$ non-linear sigma model in $1+1$ dimensions as a limit of a three-component scalar field theory restricted to the unit sphere in the large squeezing limit. This allows us to describe the model in terms of the continuous variable (CV) approach to quantum computing. Here we construct the ground state and excited states using the coupled cluster ansatz and find excellent agreement with the exact diagonalization results for a small number of lattice sites. We then present the simulation protocol for the time evolution of the model using CV gates, estimate the discretization error, and present numerical results obtained from a photonic quantum simulator. We expect that the methods developed in this work will be useful for exploring interesting dynamics for a wide class of sigma models and gauge theories, as well as for simulating scattering events on quantum hardware in the coming decade.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Continuous-time look-ahead flexible ramp scheduling in real-time operation

This paper develops a novel continuous-time look-ahead optimization model for co-optimizing balancing energy and flexible ramp products in real-time power systems operation. In addition, a continuous-time day-ahead operation model is developed to commit and schedule an adequate subset of generating units for providing balancing energy and ramping requirements in real-time operation. The up and down flexible ramp products are defined as continuous-time decision trajectories that continuously supply the real-time up and down ramping requirements of the system. The deliverability of flexible ramp products is secured through reserving the respective power requirements in generation capacity constraints of generating units. A function space solution method is proposed to reduce the dimensionality of day-ahead and real-time models through projecting them into finite-dimensional function spaces spanned by Bernstein polynomials, where the real-time operation decisions are modeled with a higher accuracy level than day-ahead decisions. Simulations are conducted on the IEEE-RTS system using CAISO load data, and the numerical results showcase the effectiveness of the proposed model in supplying the energy and flexible ramp requirements and avoiding ramping scarcity in the system, while reducing the operation costs.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS2021

The Near and Intermediate Wake of a Thick Flat Plate with a Circular Trailing Edge

The near and intermediate wake regions of a thick flat plate with a circular trailing edge (TE) are investigated with a direct numerical simulation (DNS). The focus is on the region x/D < 43.0(D is the diameter of the TE). The upper and lower separating boundary layers are both turbulent and are statistically identical; the resulting wake is symmetric in the mean. The present study is a continuation of an earlier investigation (Case AX). The earlier research dealt with shed vortex structure, in particular the persistence of the coherence of these vortices with increasing streamwise distance (except for infrequent, random disruption/rupture). In addition, across-centerline movement (ACM)of shed vortices, that is prevalent past x/D = 20.0 in a cylinder case at a comparable Reynolds number, was found to be essentially absent in Case AX. In the present investigation, an underlying cause for the occasional disruption/rupture of shed vortices is explored. Strong interactions between subsequent shed vortices in the trailing edge region is found to be an important contributor to this phenomenon. The distributions of phase-averaged streamwise and cross-stream components of intensity and the corresponding budget terms in the near wake are compared with those obtained in the intermediate wake here. In addition, the main features of the distribution of the phase-averaged velocity pressure-gradient term are compared with those of the production term, for the two intensity components. Except for a change in sign, the important features in the two budget terms are qualitatively similar for both components.

Direct Numerical Simulations of Flat Plate Wake↗
DOE OSTI2026

Characteristics of Fluid‐Solid Interaction Constitutive Models Within Poroelastodynamics at Higher Strain‐Rates and Large Deformations Implemented in 1D

The large deformation, mixed formulation, finite element (FE) modeling approach presented in Irwin et al. 2024 is extended herein to include improved constitutive models for representing dynamic solid-fluid interactions at higher strain rates (𝒪⁢(1⁢0 2 −1⁢0 3 )⁢s −1 ) and larger overpressure magnitudes (𝒪⁡(1⁢0 2 )⁢kPa) within a biphasic soft porous material using Theory of Porous Media (TPM) at finite strain. Specifically, these constitutive modeling improvements are the following: (i) a more physically robust constitutive model for pore fluid seepage velocity via inclusion of pore fluid viscous stress, and (ii) a modified deformation-dependent-permeability model and updated hyperelastic constitutive model better suited for handling larger volumetric compressions and extensions. The novelty of the present work is mainly the contribution (i): inclusion of pore fluid viscous stress at higher strain-rate and large deformations, which requires 𝐶 1 continuity in the weak formulation, accomplished by employing Hermite cubic interpolation functions within a mixed nonlinear poromechanical finite element formulation. In (ii), the model is updated to weakly enforce solid phase incompressibility, such that this assumption is not violated numerically, which provides improved numerical stability for achieving larger overpressure magnitudes on 𝒪⁡(1⁢0 2 ) kPa, which were not achievable with the previous Kozeny–Carman model in Irwin et al. 2024. Also in (ii), the volumetric part of the solid skeleton free energy function is modified to ensure proper bounds on the solid skeleton Jacobian of deformation 𝐽 s related to incompressibility of the solid phase. Uniaxial strain, unidirectional flow examples at higher strain rates (𝒪⁢(1⁢0 2 −1⁢0 3 )⁢s −1 ) and larger deformations (up to 0.2 (or 20%) nominal axial strain) demonstrate the improved physical representation—and numerical stability—of these constitutive model improvements.

42 ENGINEERING↗
DOE OSTI2020

A random unitary circuit model for black hole evaporation

Inspired by the Hayden-Preskill protocol for black hole evaporation, we consider the dynamics of a quantum many-body qudit system coupled to an external environment, where the time evolution is driven by the continuous limit of certain 2-local random unitary circuits. We study both cases where the unitaries are chosen with and without a conserved U(1) charge and focus on two aspects of the dynamics. First, we study analytically and numerically the growth of the entanglement entropy of the system, showing that two different time scales appear: one is intrinsic to the internal dynamics (the scrambling time), while the other depends on the system-environment coupling. In the presence of a U(1) conserved charge, we show that the entanglement follows a Page-like behavior in time: it begins to decrease in the middle stage of the “evaporation”, and decreases monotonically afterwards. Second, we study the time needed to retrieve information initially injected in the system from measurements on the environment qudits. Based on explicit numerical computations, we characterize such time both when the retriever has control over the initial configuration or not, showing that different scales appear in the two cases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

A hereditary integral, transient network approach to modeling permanent set and viscoelastic response in polymers

An efficient numerical framework is presented for modeling viscoelasticity and permanent set of polymers. It is based on the hereditary integral form of transient network theory, in which polymer chains belong to distinct networks each with different natural equilibrium states. Chains continually detach from previously formed networks and reattach to new networks in a state of zero stress. The free energy of these networks is given in terms of the deformation gradient relative to the configuration at which the network was born . A decomposition of the kernel for various free energies allows for a recurrence relationship to be established, bypassing the need to integrate over all time history. The technique is established for both highly compressible and nearly incompressible materials through the use of neo-Hookean, Blatz–Ko, Yeoh, and Ogden-Hill material models. Multiple examples are presented showing the ability to handle rate-dependent response and residual strains under complex loading histories.

Finite element method↗
DOE OSTI2021

On the implementation of dislocation reactions in continuum dislocation dynamics modeling of mesoscale plasticity

The continuum dislocation dynamics framework for mesoscale plasticity is intended to capture the dislocation density evolution and the deformation of crystals when subjected to mechanical loading. It does so by solving a set of transport equations for dislocations concurrently with crystal mechanics equations, with the latter being cast in the form of an eigenstrain problem. Incorporating dislocation reactions in the dislocation transport equations is essential for making such continuum dislocation dynamics predictive. A formulation is proposed to incorporate dislocation reactions in the transport equations of the vector density-based continuum dislocation dynamics. This formulation aims to rigorously enforce dislocation line continuity using the concept of virtual dislocations that close all dislocation loops involved in cross slip, annihilation, and glissile and sessile junction reactions. The addition of virtual dislocations enables us to accurately enforce the divergence free condition upon the numerical solution of the dislocation transport equations for all slip systems individually. A set of tests were performed to illustrate the accuracy of the formulation and the solution of the transport equations within the vector density-based continuum dislocation dynamics. Comparing the results from these tests with an earlier approach in which the divergence free constraint was enforced on the total dislocation density tensor or the sum of two densities when only cross slip is considered shows that the new approach yields highly accurate results. Bulk simulations were performed for a face centered cubic crystal based on the new formulation and the results were compared with discrete dislocation dynamics predictions of the same. The microstructural features obtained from continuum dislocation dynamics were also analyzed with reference to relevant experimental observations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Prospects of Nanoscience with Nanocrystals: 2025 Edition

Abstract Nanocrystals (NCs) of various compositions have made important contributions to science and technology, with their impact recognized by the 2023 Nobel Prize in Chemistry for the discovery and synthesis of semiconductor quantum dots (QDs). Over four decades of research into NCs has led to numerous advancements in diverse fields, such as optoelectronics, catalysis, energy, medicine, and recently, quantum information and computing. The last 10 years since the predecessor perspective “Prospect of Nanoscience with Nanocrystals” was published in ACS Nano have seen NC research continuously evolve, yielding critical advances in fundamental understanding and practical applications. Mechanistic insights into NC formation have translated into precision control over NC size, shape, and composition. Emerging synthesis techniques have broadened the landscape of compounds obtainable in colloidal NC form. Sophistication in surface chemistry, jointly bolstered by theoretical models and experimental findings, has facilitated refined control over NC properties and represents a trusted gateway to enhanced NC stability and processability. The assembly of NCs into superlattices, along with two-dimensional (2D) photolithography and three-dimensional (3D) printing, has expanded their utility in creating materials with tailored properties. Applications of NCs are also flourishing, consolidating progress in fields targeted early on, such as optoelectronics and catalysis, and extending into areas ranging from quantum technology to phase-change memories. In this perspective, we review the extensive progress in research on NCs over the past decade and highlight key areas where future research may bring further breakthroughs.

Chemistry↗
DOE OSTI2025

Operator dynamics in Floquet many-body systems

We study operator dynamics in many-body quantum systems, focusing on generic features of systems that are ergodic, spatially extended, and lack conserved densities. Quantum circuits of various types provide simple models for such systems. We focus on Floquet quantum circuits, comparing their behavior with what has been found previously for circuits that are random in time. Floquet circuits, which have discrete time-translation symmetry, represent an intermediate case between circuits that are random in time and lack any symmetry, and systems with a time-independent Hamiltonian and continuous time-translation invariance. By making this comparison, one of our aims is to identify signatures of time-translation symmetry in Floquet operator dynamics. To characterize behavior we examine a variety of quantities in solvable models and numerically: operator autocorrelation functions; the partial spectral form factor; the out-of-time-order correlator (OTOC); and the paths in operator space that make the dominant contributions to the ensemble-averaged autocorrelation functions. Our most striking result is that ensemble-averaged autocorrelation functions show behavior that is distinctively different in Floquet systems compared to systems in which successive time-steps are independent. Specifically, while average autocorrelation functions decay on a microscopic timescale for circuits that are random in time, in Floquet systems they have a late-time tail with a duration that grows parametrically with the size of the operator support. In the simplest models this tail is separated from the initial decay by a minimum, so that the average autocorrelation function has an intermediate-time peak. The existence of these tails provides a way to understand deviations of the spectral form factor from random matrix behavior at times shorter than the Thouless time. In contrast to this feature in autocorrelation functions, we find no new aspects to the behavior of OTOCs for Floquet models compared to random-in-time circuits. We show that this difference between averaged autocorrelation functions and OTOCs can be understood in terms of the paths in operator space that contribute to the two quantities: paths for the former retain a limited support at late times, while paths for the latter are dominated by operator spreading. Published by the American Physical Society 2025

Yoshimura, Takato (ORCID:0000000309159846)↗
DOE OSTI2022

Generation of polarization entanglement via the quantum Zeno effect

The quantum Zeno effect reveals that continuous observation of a quantum system can significantly alter its evolution. Here, we present a method for establishing polarization entanglement between two initially unentangled photons in coupled waveguides via the quantum Zeno effect. We support our analytical investigation with numerical simulations of the underlying Schrodinger equation describing the system. Further, we extend our technique to three coupled waveguides in a planar configuration and determine the parameters required to generate three-qubit W-states. In contrast to existing schemes based on a vacuum and single-photon encoding, the polarization encoding in our approach is compatible with quantum information protocols that remove photon loss through post-selection. Our findings offer a powerful quantum state engineering approach for photonic quantum information technologies.

Nodurft, Ian C.↗