Fact-checked by Grok 2 weeks ago

Boundary conditions in fluid dynamics

In , boundary conditions are the mathematical constraints imposed on the governing partial differential equations, such as the Navier-Stokes equations, at the edges of the computational or physical to yield a unique and physically meaningful solution for variables like , , and . These conditions reflect real-world interactions at interfaces, including solid surfaces, fluid-fluid boundaries, inlets, and outlets, and are essential for closing the that describe fluid motion. Boundary conditions play a pivotal role in both theoretical analysis and (CFD), where they directly influence , , and the accuracy of simulations by enforcing physical realism and preventing ill-posed problems. Their specification is particularly critical in viscous flows, where small perturbations at boundaries can amplify through mechanisms like transition, affecting , , and overall flow behavior by orders of magnitude. In CFD applications, improper boundary conditions can introduce significant uncertainties, especially in complex scenarios involving multiphase flows or high-speed , underscoring the need for precise measurement and modeling of real-world inputs like or inlet disturbances. Common types of boundary conditions in fluid dynamics include Dirichlet conditions, which prescribe fixed values for variables (e.g., at an inlet), Neumann conditions, which specify gradients or fluxes (e.g., zero normal stress at a ), and mixed or Robin conditions combining both. For solid walls, the is fundamental in viscous flows, requiring the fluid to match the wall (typically zero for stationary walls), arising from viscous adhesion effects validated by macroscopic experiments like capillary flow. In contrast, slip conditions apply to inviscid or rarefied flows, allowing tangential discontinuity, while conditions enforce zero gradients across planes of to reduce computational domain size. At fluid-fluid interfaces, continuity conditions ensure matching of and tangential stress, often incorporating jumps via the normal stress balance, as in droplet or bubble dynamics. For inlets and outlets, conditions like specified uniform or extrapolated gradients maintain mass conservation and allow to propagate without reflection. These boundary conditions must be tailored to the flow regime—viscous versus inviscid, compressible versus incompressible—and integrated into numerical methods like finite volume or finite element schemes to capture phenomena such as boundary layers, where velocity gradients near walls dominate and separation. Advances in CFD techniques, including exact solutions for inflow/outflow and supersonic boundaries, continue to refine their implementation for high-fidelity predictions in , environmental, and biomedical applications.

Fundamentals

Definition and Role

Boundary conditions in fluid dynamics refer to the mathematical specifications of flow variables, such as and , at the edges of a physical or computational domain, enabling the unique solution of governing partial differential equations like the Navier-Stokes equations. These conditions dictate how the fluid interacts with its surroundings, ensuring that the mathematical model aligns with physical reality by constraining the otherwise of equations. The origins of boundary conditions trace back to foundational problems in partial differential equations, including the , which addresses initial value specifications for hyperbolic equations, and the , which imposes fixed values on boundaries for elliptic equations; these concepts were adapted to during the by key figures such as and George Gabriel Stokes. Navier's 1822 work on viscous flow introduced early boundary assumptions, while Stokes's 1845 derivation of the incorporated stress-based constraints using and symmetry considerations, laying the groundwork for modern fluid boundary formulations. Their role is indispensable because, absent boundary conditions, solutions to the Navier-Stokes equations remain non-unique or ill-posed, as the equations describe infinite possible states without domain-specific constraints; these conditions physically represent interfaces like solid walls or open inlets, translating real-world constraints into solvable problems. Boundaries in fluid dynamics are distinguished as physical, which correspond to tangible interfaces such as vessel walls, or artificial, which are imposed in computational simulations to approximate infinite domains by truncating the solution space. A representative example is steady in a , where prescribing a uniform or parabolic profile at the ensures a determinate to the Navier-Stokes equations, averting non-physical infinite or arbitrary fields that could otherwise emerge. The exemplifies a classic physical at walls, where fluid equals the stationary surface .

Classification Schemes

Boundary conditions in fluid dynamics are systematically classified into mathematical categories based on how they constrain the solution variables at the domain boundaries. The primary mathematical classifications are Dirichlet, Neumann, and Robin (or mixed) conditions, which arise from the need to solve partial differential equations governing fluid flow. Dirichlet boundary conditions specify the exact value of the dependent variable on the boundary Γ. For a scalar field u, this takes the general form u = g on Γ, where g is a prescribed function. In fluid dynamics, this often applies to velocity components at solid walls, such as the no-slip condition where velocity is fixed to zero relative to the wall. Neumann boundary conditions specify the normal derivative (or ) of the variable on the boundary, given by ∂u/∂n = h on Γ, where n is the outward normal and h is prescribed. This type is common for outflow boundaries or when specifying or , as it enforces a rather than an . Robin boundary conditions combine Dirichlet and Neumann types through a , typically αu + β ∂u/∂n = f on Γ, where α and β are constants. These are useful for convective at walls or coupled fluid-structure interactions, balancing value and flux. Physically, boundary conditions are classified according to the type of or the . Boundaries are categorized as solid walls, free surfaces, inlets, or outlets, each imposing conditions that reflect the interaction with the surrounding environment. For instance, solid walls typically enforce no-penetration and tangential constraints, while inlets specify incoming properties and outlets allow free exit. Another physical classification distinguishes conditions based on : viscous flows require detailed enforcement of effects, such as no-slip at walls, whereas inviscid flows use slip conditions that neglect tangential . This distinction is crucial for high-Reynolds-number approximations in Euler equations versus full Navier-Stokes solutions. In numerical implementations, particularly for incompressible flows, boundary conditions are tailored to the pressure-velocity coupling. Dirichlet conditions are commonly applied to fields at walls and inlets to enforce specified values, ensuring mass conservation. Neumann conditions are typically used for , such as zero-gradient at outlets or walls, to avoid over-constraining the for pressure correction in methods. This combination prevents singularities in the pressure solver while satisfying the divergence-free constraint. A common pitfall in applying these classifications is over-specification, where too many conditions (e.g., both and at an ) are imposed, leading to inconsistencies, ill-posed problems, or solver in simulations. Careful selection based on physical principles and is essential to avoid such issues.

Wall Boundary Conditions

No-Slip Condition

The no-slip condition arises from the physical of viscous molecules to a , primarily due to intermolecular forces and , which prevent relative tangential motion at the . This creates a sharp gradient near the wall, where particles in direct contact with the stationary exhibit zero relative to the surface. In most macroscopic flows of liquids and dense gases, this condition holds as an empirical observation supported by extensive experiments, contrasting with inviscid flows where slip may occur. Mathematically, the no-slip condition is expressed as the fluid velocity vector \mathbf{u} equaling the wall velocity \mathbf{U}_w at the boundary, or \mathbf{u} = \mathbf{U}_w for all components, with the normal component ensuring impermeability (\mathbf{u} \cdot \mathbf{n} = \mathbf{U}_w \cdot \mathbf{n}) and the tangential components enforcing zero slip. For stationary walls, this simplifies to \mathbf{u} = 0 at the surface. This boundary condition is imposed in the solution of the Navier-Stokes equations, as the equations themselves are second-order in spatial derivatives and require such supplemental constraints; it does not directly follow from the equations but is justified by the continuity of the stress tensor at the fluid-solid , where viscous stresses balance to accommodate the velocity gradient without discontinuity. Historically, the was first proposed by in 1823 as part of his derivation of the equations governing viscous fluid motion, allowing for a general slip model that includes the no-slip limit. It was later rigorously justified and validated experimentally by George Gabriel Stokes in the 1840s through studies of oscillatory fluid motion and shear flows, confirming the condition's applicability in viscous regimes. In applications, the is fundamental to modeling laminar and turbulent s over solid surfaces, where it enables accurate prediction of and momentum transfer. For instance, Ludwig Prandtl's 1904 boundary layer theory relies on no-slip to resolve the transition from zero wall velocity to free-stream conditions, facilitating calculations in for airfoils and vehicle bodies. This has high impact in , as it underpins viscous estimation in high-Reynolds-number flows, contributing to designs that minimize consumption in . The fails in rarefied gas flows, such as those in microchannels or high-altitude aerodynamics, where the Kn > 0.01 indicates significant relative to system scale, leading to velocity slip at due to insufficient molecular collisions. In such cases, slip boundary conditions, like the Navier slip model, must replace no-slip to capture partial . This contrasts with the full assumed in viscous flows, though slip is also relevant in high-speed or superhydrophobic surface scenarios.

Slip Condition

The slip boundary in applies to situations where the fluid does not adhere completely to a solid wall, permitting tangential while enforcing no penetration. This arises physically from the absence of viscous shear stress at the wall, which is characteristic of inviscid flows described by the Euler equations, where molecular interactions are negligible compared to inertial forces. It also holds in rarefied gas flows, such as low-density microflows with high Knudsen numbers (Kn > 0.01), where the mean free path of gas molecules is significant relative to the scale, leading to incomplete momentum accommodation through a combination of diffuse and specular reflections at the surface. Mathematically, the basic slip condition for inviscid flows specifies zero and zero tangential : \mathbf{u} \cdot \mathbf{n} = 0 \boldsymbol{\tau} \cdot \mathbf{t} = 0 where \mathbf{u} is the fluid , \mathbf{n} the unit to the wall, \mathbf{t} the unit tangential vector, and \boldsymbol{\tau} the . For transitional slip regimes in rarefied gases, the slip provides a correction, relating the tangential slip to the wall- : u_t = \lambda \frac{\partial u_t}{\partial n} with slip length \lambda = \frac{2 - \sigma}{\sigma} \ell, where \sigma is the tangential momentum accommodation coefficient (0 < \sigma ≤ 1) and \ell the molecular mean free path. This boundary condition finds applications in high-Reynolds-number approximations, such as hypersonic flows around re-entry vehicles, where rarefaction effects (Kn ≈ 0.001–0.1) necessitate slip to capture aerodynamic heating and drag accurately without resolving full viscous layers. In microfluidics, it models gas or liquid transport in microchannels and nanopores, enabling predictions of enhanced flow rates and drag reduction on engineered surfaces like superhydrophobic coatings, where slip lengths up to 50 nm can boost efficiency by factors of up to 20. Adaptations to the Navier-Stokes equations involve directly applying the slip formulation at walls for rarefied or multiphase simulations, or approximating inviscid limits by setting molecular viscosity to zero, which recovers Euler-like behavior. In turbulent flows, wall functions in models like k-ε integrate the slip condition implicitly by bridging the viscous sublayer to the logarithmic region, allowing coarser meshes while estimating near-wall slip velocities. Despite its utility, the slip condition is limited in accuracy for flows where viscous diffusion dominates near walls, such as low-Reynolds-number regimes, as it neglects the detailed momentum transfer that the captures more effectively.

Flow Entry and Exit Conditions

Inlet Boundary Conditions

Inlet boundary conditions in fluid dynamics are applied at the entry points of the computational domain to prescribe the incoming flow properties, ensuring accurate simulation of fluid behavior as it enters the system. These conditions are essential in (CFD) to define the upstream flow state, typically specifying , , or while allowing other variables like or scalars to be extrapolated from the interior solution. The primary types of inlet boundary conditions include inlet, mass flow inlet, and inlet. In a inlet condition, the flow components are directly specified at the boundary, often assuming a uniform profile where the component is set to a prescribed value U_0 directed inward (i.e., \mathbf{u} = -U_0 \mathbf{n}, with \mathbf{n} as the outward unit ), and tangential components are zero for simplicity; alternatively, a parabolic profile may be used for fully developed pipe flows to represent realistic distribution. A mass flow inlet condition fixes the \rho \mathbf{u} \cdot \mathbf{n} = \dot{m}/A, where \rho is , \dot{m} is the , and A is the inlet area, allowing to adjust based on local for compressible flows. The inlet condition specifies a fixed total or p = p_{\text{in}}, with and other properties extrapolated from the interior, making it suitable when inlet is unknown but is controlled, such as in external flows. For turbulent flows, inlet boundary conditions must also specify turbulence quantities to initialize the turbulent structures accurately. Common parameters include turbulence intensity I, defined as the ratio of the root-mean-square of velocity fluctuations to the mean velocity (typically 5-10% for fully turbulent flows), and the integral length scale L_t, which represents the size of the largest energy-containing eddies (often estimated as 0.07 times the for pipe flows). These are used to compute turbulent kinetic energy k = \frac{3}{2} (I U_0)^2 and dissipation rate \epsilon = C_\mu^{3/4} \frac{k^{3/2}}{L_t} (with C_\mu \approx 0.09 for standard k-\epsilon models), ensuring downstream development matches physical expectations. Inlet boundary conditions find applications in simulating duct flows, where or inlets enforce uniform or developed profiles to study drops and , and in simulations, where inlets with specified initiate free layers for analyzing mixing and . These conditions help maintain overall conservation in closed domains when paired with appropriate outlets. A key challenge in applying boundary conditions is specifying realistic and profiles, as oversimplified uniform assumptions can lead to artificial upstream influence or unphysical reflections that distort the flow field, particularly in high-Reynolds-number simulations requiring detailed measurements or synthetic turbulence generation.

Outlet Boundary Conditions

Outlet boundary conditions in (CFD) are essential for modeling fluid exits in a manner that minimizes upstream influence and prevents unphysical reflections of waves or disturbances. These conditions are particularly crucial in simulations where the flow exits the into an unbounded or far-field region, ensuring and physical accuracy. They serve as a complement to inlet boundary conditions to fully specify the problem . Common types of outlet boundary conditions include pressure outlets and convective outlets. In a pressure outlet condition, the is fixed at a specified value, while velocities and other variables are extrapolated from the interior , often assuming zero normal for outflow. This approach is suitable for flows where the exit is known or ambient. Convective outlets, on the other hand, model the of disturbances out of the using a wave-like , given by the equation \frac{\partial u}{\partial t} + U \frac{\partial u}{\partial x} = 0, where u is the perturbation, t is time, U is the convective speed (typically the local speed), and x is the normal to the boundary. This condition allows outgoing to pass through without , making it ideal for unsteady flows. Mathematically, outlet conditions often invoke a zero-gradient for outflow variables, expressed as \frac{\partial \phi}{\partial n} = 0, where \phi represents quantities like components and \mathbf{n} is the outward . Alternatively, a fully developed may be applied, positing that the profile no longer evolves in the streamwise at the , which justifies zero gradients for all transported variables except pressure. These formulations ensure mass conservation and prevent by adjusting based on local . Outlet boundary conditions find applications in simulating exhaust pipes, where pressure outlets model discharge into ambient conditions while preventing recirculation; in wake regions behind bluff bodies, convective conditions handle without domain truncation errors; and in general CFD setups to avoid that could destabilize iterations. For numerical stability in compressible flows, non-reflecting boundaries based on characteristic variables are employed, decomposing the flow into acoustic, , and waves to specify only outgoing characteristics at the boundary, thus minimizing reflections. This approach enhances accuracy in high-speed simulations like jet exhausts. The development of outlet boundary conditions has evolved significantly, with the convective condition originating from Orlanski's 1976 work on simple treatments for unbounded flows, which introduced phase-speed-based to handle open boundaries effectively. Subsequent advancements, such as characteristic-based methods, built on this foundation to address viscous and reacting flows.

Constant Pressure Condition

The constant boundary condition in models interfaces where the domain connects to a large or ambient that maintains a uniform and fixed , such as an open outlet discharging into the atmosphere or a confined adjacent to a pressure-stabilized region. This condition assumes that the reservoir's size and dynamics impose negligible perturbations on the local at the boundary, effectively treating it as invariant. Physically, it arises in situations where pressure equalization occurs rapidly due to the incompressibility of the or the dominance of external forcing, ensuring that the boundary does not reflect waves or disturbances back into the domain. Mathematically, the condition is formulated as a Dirichlet boundary condition for pressure, specified as p = p_0 on the boundary surface, where p_0 is the prescribed constant value. For velocity, it is typically paired with a zero normal gradient condition, \frac{\partial (\mathbf{u} \cdot \mathbf{n})}{\partial n} = 0, where \mathbf{u} is the velocity vector and \mathbf{n} is the outward unit normal, allowing fluid to exit naturally without imposed resistance. This combination prevents artificial reflections and maintains physical consistency at open boundaries. In applications, the constant pressure condition is widely used for simulating free jets exiting into ambient air at fixed atmospheric pressure, atmospheric boundaries in geophysical flows where external pressure dominates, and outlets in incompressible pipe or channel flow computations to model unrestricted discharge. For instance, in buoyancy-driven cavity flows or multi-outlet internal systems, it facilitates accurate representation of pressure-driven mass outflow without specifying detailed downstream geometry. Within the incompressible Navier-Stokes equations, the condition plays a critical role in ensuring the solvability of the pressure Poisson equation, which enforces the divergence-free velocity constraint \nabla \cdot \mathbf{u} = 0. By providing Dirichlet data for on outflow portions of the boundary, it completes the for the elliptic pressure equation, yielding a unique solution and preventing ill-posedness in numerical schemes like methods. Without such specification, the field could admit non-unique solutions due to the elliptic nature of the equation. Variations of the condition include setting the reference p_0 = 0 to define a gauge field, which eliminates arbitrary additive constants and mitigates numerical singularities in relative formulations common to incompressible flows. This approach is particularly useful in closed domains with partial open boundaries, where absolute levels are irrelevant but relative gradients drive the flow.

Specified Pressure Gradient Condition

The specified pressure gradient condition constitutes a Neumann-type boundary condition in fluid dynamics, where the normal derivative of the pressure field is prescribed at the boundary to enforce a particular pressure flux or flow acceleration. This approach is particularly useful for scenarios requiring controlled mass flow rates without directly specifying inlet or outlet velocities, as it allows the pressure to evolve naturally while imposing a directional force on the fluid. Physically, it represents situations where an external mechanism, such as a pump, gravity, or imposed shear, drives the flow through a consistent change in pressure across the domain boundary, thereby influencing the momentum balance at the interface. Mathematically, the condition is formulated as \frac{\partial p}{\partial n} = \frac{dp}{dn}, where p is the , n denotes the outward direction to the , and \frac{dp}{dn} is the specified value. A common example arises in , where the vertical balances gravitational forces, given by \frac{\partial p}{\partial z} = -\rho g, with \rho as the and g as ; this ensures the pressure increases hydrostatically with depth without inducing spurious velocities. In numerical implementations, such as finite volume methods, this is discretized to compute the through faces, maintaining consistency with the laws. This boundary condition finds applications in gravity-driven flows, where the hydrostatic gradient simulates natural descent under body forces, as in thin-film coatings down inclined surfaces; in porous media outlets, it models steady seepage with Darcy's law by linking pressure drop to permeability and flow rate; and in wind tunnel tests, where a uniform streamwise gradient replicates controlled free-stream conditions for aerodynamic validation. In fully developed channel flows, a constant adverse or favorable pressure gradient drives turbulence, enabling direct numerical simulations to study wall-bounded shear layers. The coupling between the specified pressure gradient and velocity arises directly from the momentum equation in the Navier-Stokes system, where the pressure term -\frac{1}{\rho} \nabla p acts as a driving force. At the boundary, this yields \frac{\partial p}{\partial n} = -\rho \left( \frac{\partial \mathbf{u}}{\partial t} + (\mathbf{u} \cdot \nabla) \mathbf{u} \right) \cdot \mathbf{n}, with \mathbf{u} as the vector; in steady, fully developed flows, it simplifies to balancing viscous against the imposed , ensuring compatibility with no-slip walls or periodic domains. This derivation guarantees well-posedness for incompressible flows by projecting the boundary-normal onto the . Compared to constant pressure conditions, the specified offers advantages in non-uniform , such as those with varying or multiphase interfaces, by permitting local adjustments that prevent artificial reflections or instabilities while precisely controlling bulk rates. It is especially effective in periodic setups, reducing computational size without compromising physical realism.

Symmetry and Geometric Conditions

Symmetric Boundary Condition

The in exploits mirror-like reflection across a of to model behavior, ensuring that the computational domain represents only a portion of the full physical setup while maintaining equivalence to the complete system. This condition is physically based on the assumption that the and field exhibit perfect bilateral , preventing any penetration of across the and enforcing reflective properties for velocities and scalars. It is particularly useful in reducing computational costs by simulating half or quarter domains for symmetric configurations, such as in external or problems where the remains invariant under reflection. Mathematically, the symmetric imposes a zero component at the , expressed as \mathbf{u} \cdot \mathbf{n} = 0, where \mathbf{u} is the and \mathbf{n} is the unit to the plane. Additionally, it requires zero gradients for all transported variables, including scalars like or species concentration, given by \frac{\partial \phi}{\partial n} = 0, where \phi represents any ; this ensures zero convective and diffusive across the plane. For viscous flows, the also implies zero in the tangential directions, effectively mimicking an inviscid slip wall without frictional effects to the . These formulations combine Dirichlet conditions for the and conditions for tangential components and gradients, preserving the of the solution. In applications, symmetric boundary conditions are commonly applied to half-domain simulations of axisymmetric bodies, such as aircraft wings, where the root plane serves as the symmetry boundary to capture and characteristics without modeling the full span. For instance, in (CFD) analyses of symmetric airfoils like the NACA 0015 or high-aspect-ratio wings, this approach halves the mesh size while accurately predicting symmetric flow fields around the mid-span. It is also employed in internal flows, like pipe bends or heat exchangers with symmetric inlets, to streamline simulations of laminar or turbulent regimes. An extension to axisymmetric geometries can be considered for cylindrical setups, but planar symmetry remains distinct for reductions. In CFD implementation, symmetric boundaries are typically treated as special wall types that enforce slip conditions, where the solver applies the zero-normal-velocity constraint and extrapolates tangential velocities without viscous drag, often using finite volume or finite element methods to maintain balance. Software like Fluent or configures these via boundary zone definitions, automatically handling the reflective mapping for variables during iterations. This setup is computationally efficient, as it avoids resolving the mirrored domain, but requires precise alignment of the plane with the geometric symmetry axis. Limitations arise when the flow field deviates from perfect , such as in cases involving asymmetric perturbations like trailing wakes, buoyancy-driven flows, or swirl, where the condition would artificially constrain the solution and lead to inaccuracies. For example, in wing simulations, it fails downstream of the trailing edge if introduces asymmetry, necessitating full-domain modeling instead. Thus, validation against experimental data or full simulations is essential to confirm the holds throughout the .

Axisymmetric Boundary Condition

The axisymmetric boundary condition is employed in fluid dynamics for flows that possess rotational symmetry about a central axis, where physical quantities such as velocity and pressure show no dependence on the azimuthal angle θ in cylindrical coordinates (r, θ, z). This assumption simplifies three-dimensional problems to two-dimensional computations in the r-z plane, reducing computational cost while maintaining accuracy for rotationally invariant geometries. The condition physically implies that the flow is invariant under rotation around the axis, eliminating azimuthal gradients and setting the azimuthal velocity component u_θ to zero for non-swirling cases. Mathematically, the axisymmetric condition enforces ∂/∂θ = 0 for all flow variables, with boundary conditions specifying the u_r and axial velocity u_z on surfaces while u_θ = 0. On the symmetry axis (r = 0), additional constraints apply, such as u_r = 0 and ∂u_z/∂r = 0, to ensure physical consistency and prevent singularities. The governing equations are the axisymmetric forms of the Navier-Stokes equations in cylindrical coordinates, incorporating geometric factors like 1/r terms. For , the becomes \frac{1}{r} \frac{\partial (r u_r)}{\partial r} + \frac{\partial u_z}{\partial z} = 0, while the momentum equations include centrifugal and Coriolis-like terms absent in Cartesian forms, such as -u_θ²/r in the radial direction for swirling flows. These modifications account for the curvature of the coordinate system and are derived from the full three-dimensional equations by applying the azimuthal invariance. This boundary condition finds applications in modeling internal flows through pipes, where fully developed laminar flow (Hagen-Poiseuille) exhibits perfect axisymmetry, enabling efficient prediction of pressure drops and velocity profiles. It is also widely used for nozzle flows, such as in rocket or jet propulsion systems, to simulate compressible expansion without azimuthal variations. Additionally, axisymmetric conditions apply to bubble dynamics, including the rise of spherical or oblate bubbles in liquids, where symmetry simplifies tracking of wake structures and deformation. In numerical solvers, the axisymmetric involves transforming Cartesian-based codes to cylindrical coordinates or using dedicated axisymmetric modules, where the is treated as a meridional slice with conditions enforcing . This transformation preserves the physics by integrating the 2π azimuthal periodicity implicitly, avoiding full meshes for symmetric problems.

Cyclic and Repeating Conditions

Periodic Boundary Condition

Periodic boundary conditions in fluid dynamics are employed to simulate flows that exhibit repeating patterns across the boundaries of the computational domain, effectively modeling infinite or extended domains without explicitly resolving the entire extent. Physically, these conditions assume that the flow field repeats identically after a spatial period L, such that the state at one boundary matches the state at the opposite boundary, as if the fluid exiting one side re-enters the other unchanged. This approach is particularly suited to scenarios like periodic geometries in systems, where the flow in one or module mirrors adjacent ones, such as in crystal-like structures for porous media flows or repeating stages in . Mathematically, periodic boundary conditions enforce equality of flow variables across paired boundaries, expressed as \mathbf{u}(\mathbf{x}) = \mathbf{u}(\mathbf{x} + \mathbf{L}) for \mathbf{u}, and similarly for , , and other scalars, where \mathbf{L} is the periodicity vector defining the repeat distance. In numerical implementations, such as finite volume or methods, this is achieved by pairing inlet-outlet boundaries and copying values from one to the other, often using ghost cells outside the domain to maintain continuity during operations. For instance, ghost cells adjacent to the boundary are populated with values mirrored from the opposite side, ensuring smooth flux calculations without discontinuities. This method preserves the periodic structure while facilitating efficient computation on reduced domains. Two primary types of periodic boundary conditions are distinguished: simple periodic, which applies direct translation in the same direction (e.g., streamwise periodicity in flows), and spanwise periodic, which enforces repetition in orthogonal directions (e.g., across the span in three-dimensional simulations). Simple periodic conditions are straightforward for one-dimensional repeats, while spanwise variants handle multi-directional periodicity, such as in simulations of cross-flow over periodic arrays. These types allow flexibility in modeling anisotropic repeats without altering the core enforcement mechanism. Applications of periodic boundary conditions are widespread in computational fluid dynamics, including direct numerical simulations of homogeneous isotropic within cuboidal domains to study statistical properties over long times, heat exchanger modules where flow through repeating tube banks is analyzed for convective , and fully developed flows to isolate effects from entrance behaviors. In multi-stage compressors, they enable modeling of repeating passages at various Reynolds numbers, capturing unsteady flow features like wakes without simulating the full machine. These conditions ensure , , and across the paired boundaries, as the net flux into the domain equals the net flux out, maintaining global balances essential for accurate long-term simulations. Periodic boundary conditions share conceptual similarities with cyclic conditions for handling mapped interfaces but assume direct, unmapped translations.

Cyclic Boundary Condition

The cyclic boundary condition addresses flows with non-uniform periodicity, where the flow field repeats across boundaries subject to a , such as or scaling, rather than pure . This condition is particularly relevant in scenarios involving , where the physical flow in adjacent sectors of a is identical up to a transformation around an axis of . By applying this boundary condition, computational can be reduced to a single or few sectors, significantly lowering computational cost while preserving the essential flow physics. The physical basis stems from the invariance of the under the specified ; for instance, in rotating flows, the and at a point on one match the transformed values at the corresponding point on the paired , ensuring of the across the . This accounts for mismatches in geometry or motion, such as differing angular positions or scales between connected regions. In rotor-stator configurations, the condition enforces that the relative motion between rotating and stationary components is captured without simulating the full 360-degree annulus. Mathematically, the cyclic boundary condition is expressed as \mathbf{u}(\mathbf{x}) = R \mathbf{u}(\mathbf{x} + \mathbf{L}), where \mathbf{u} is the , \mathbf{x} is the , R is the corresponding to the periodicity angle, and \mathbf{L} represents the displacement vector along the periodicity direction. For scalar fields like p, the condition simplifies to p(\mathbf{x}) = p(\mathbf{x} + \mathbf{L}), while vector quantities undergo the full transformation. At interfaces, flux matching ensures , such that the and fluxes across the boundaries are equal after transformation, preventing artificial discontinuities. In applications, cyclic boundary conditions are widely used for rotor-stator interfaces in , such as axial compressors and turbines, where only a single blade passage is simulated by pairing azimuthal boundaries separated by the angle. They also apply to helical flows in devices like screw pumps or marine propellers, capturing the combined axial translation and rotational periodicity. These conditions enable accurate prediction of unsteady interactions, such as wake passing, in reduced-domain simulations. Implementation in (CFD) relies on general grid interface (GGI) methods to handle non-conformal meshes, where grid points on paired boundaries do not align. The GGI approach uses weighted least-squares to transfer variables and fluxes between mismatched faces, maintaining second-order accuracy and properties. This is particularly essential for dynamic simulations involving relative motion, as in sliding rotor-stator setups. The method was developed for to facilitate analyses, allowing seamless coupling of rotating and stationary domains. Unlike simple periodic conditions that require conformal grids and pure translational repetition, cyclic boundary conditions accommodate geometric mismatches and transformations like , enabling more flexible meshing and realistic modeling of complex, non-uniform periodic flows.

Advanced and Specialized Conditions

Far-Field Boundary Condition

In simulations of unbounded domains, far-field conditions are imposed at artificial outer to approximate free-stream or quiescent conditions sufficiently distant from the primary flow features, thereby enabling the modeling of infinite extents within finite computational grids. This approach is essential for capturing the asymptotic behavior of disturbances that decay away from sources, such as wakes or pressure waves, without introducing spurious reflections that could contaminate the solution. The physical basis relies on the principle that, at large distances, the flow approaches a uniform free-stream state \mathbf{U}_\infty for external or a quiescent state () for low-speed or diffusive processes, ensuring that the minimally influences the near-field dynamics. Mathematically, far-field conditions often prescribe the velocity as \mathbf{u} = \mathbf{U}_\infty + O(1/r), where perturbations decay inversely with the radial distance r from the origin, reflecting the far-field expansion of the governing equations for compressible flows. For characteristic-based formulations, Riemann invariants are employed to decompose the flow into incoming and outgoing waves along characteristic lines, allowing non-reflecting treatment by extrapolating incoming invariants from the free-stream while advecting outgoing ones. A prominent type is the , adapted for hyperbolic flows as \frac{\partial \phi}{\partial t} + c \frac{\partial \phi}{\partial x} = 0, where \phi is a , c is the wave speed, and the form ensures outgoing waves propagate without reflection; this is particularly effective for time-dependent wave-like phenomena in external flows. These formulations derive from of the linearized Euler or Navier-Stokes equations, prioritizing low-order accuracy to balance computational cost and . Far-field boundary conditions find primary applications in external aerodynamic simulations, such as flows around wings or vehicles, where the domain must extend to to resolve wake decay and far-wake interactions accurately. They are also crucial for modeling dispersion in atmospheric or environmental flows, simulating the transport and dilution of contaminants over large, open without boundary-induced recirculation. For instance, in air quality assessments, these conditions enable the representation of unbounded and processes. Numerically, implementing far-field conditions often involves auxiliary techniques like sponge layers, which introduce artificial damping terms—typically Rayleigh or artificial viscosity—into the momentum and energy equations within a buffer zone adjacent to the boundary. This absorption mechanism gradually forces the solution toward the target free-stream or quiescent state, effectively dissipating outgoing waves and minimizing reflections, with optimal layer thickness and damping profiles tuned for broadband frequency response in compressible simulations. Such methods enhance stability in high-fidelity large-eddy or direct numerical simulations of turbulent external flows.

Free-Surface Boundary Condition

The free-surface boundary condition governs the behavior at deformable interfaces between a and another , such as gas or immiscible , where no solid constraint is present. Physically, this condition arises from the balance of stresses across the and the requirement that the surface evolves kinematically with the fluid motion, allowing for deformation under forces like , , and . Unlike rigid boundaries, the free surface permits normal displacement and tangential slip, enabling phenomena where the shape changes dynamically. Mathematically, the kinematic condition ensures that fluid particles on the remain there, expressed as the of a defining the surface being zero: \frac{D F}{D t} = 0, where F = 0 delineates the (e.g., F(\mathbf{x}, t) = \eta(\mathbf{x}, t) - z for a surface at \eta). This enforces that the normal component of the fluid velocity matches the interface velocity. The dynamic condition balances the stresses, given by the jump in the stress tensor across the : [\![\boldsymbol{\tau} \cdot \mathbf{n}]\!] = \sigma \kappa \mathbf{n}, where \boldsymbol{\tau} is the , \mathbf{n} is the unit normal, \sigma is the coefficient, and \kappa is the (from the Young-Laplace equation). For inviscid flows without surface tension, this simplifies to constant at the . These conditions find applications in modeling ocean waves, where gravity-driven surface deformations propagate energy across large scales; droplet dynamics, involving capillary-driven motion and breakup in multiphase flows; and boiling flows, capturing bubble growth and departure at heated surfaces through phase-change interfaces. In ocean wave simulations, the conditions enable prediction of wave propagation and breaking, critical for . For droplets, they describe impact, spreading, and ejection in or spray combustion. In boiling, they model vapor-liquid interfaces during , influencing efficiency in thermal systems. In computational fluid dynamics (CFD), free-surface conditions are coupled with interface-tracking methods like the level-set or volume-of-fluid (VOF) approaches to resolve the deformable boundary. The level-set method embeds the interface as the zero level of a signed distance function, advected by the flow while enforcing the boundary conditions through reinitialization and stress projections. VOF tracks the interface via volume fractions in Eulerian grids, applying the conditions at reconstructed surfaces to maintain sharpness and mass conservation. These methods, widely adopted since their development in the 1980s and 1990s, facilitate simulations of complex free-surface evolutions. For small-amplitude perturbations, a linearized version simplifies analysis, particularly in wave problems: \frac{\partial \eta}{\partial t} = u_n at z = 0, where u_n is the normal velocity component and the interface is approximated as flat. This form, derived via Taylor expansion, decouples the conditions and enables analytical solutions for linear wave dispersion.

References

  1. [1]
    None
    ### Summary of Boundary Conditions in Fluid Mechanics
  2. [2]
    [PDF] Fundamentals of Computational Fluid Dynamics
    Aug 26, 1999 · ... Boundary Con- ditions. Once the problem has been specified, an appropriate set of governing equations and boundary conditions must be selected.
  3. [3]
    [PDF] Boundary and Initial Condition Issues for Computational Fluid ...
    Jun 1, 2021 · Initial and boundary conditions are a source of uncertainty in Mod-Sim, and their physical specification is important for sensitive problems. ...
  4. [4]
    [PDF] Common Boundary Conditions in Fluid Mechanics
    Common boundary conditions include no-slip at the wall, symmetry, stress continuity, and velocity continuity.
  5. [5]
    [PDF] Initial and Boundary Conditions for Viscous-Flow Problems - Stern Lab
    Initial conditions. Boundary conditions: definition and types. Examples of Boundary conditions. ▫ Solid wall: no slip, FSI, moving contact line. ▫ Single ...
  6. [6]
    [PDF] Code verification of boundary conditions for compressible and ...
    Oct 15, 2014 · This work presents techniques for code verification of boundary conditions in CFD codes, including subsonic inflow/outflow, supersonic outflow, ...
  7. [7]
    What Are Boundary Conditions? Numerics Background | SimScale
    Aug 11, 2023 · Boundary conditions are constraints necessary for the solution of a boundary value problem. It is opposed to the initial value problem.Missing: Stokes | Show results with:Stokes
  8. [8]
    [PDF] The First Five Births of the Navier-Stokes Equation
    The Navier-Stokes equation is now regarded as the universal basis of fluid ... the boundary conditions. For simplicity Navier assumed that the local ...
  9. [9]
    [PDF] Computational Fluid Dynamics
    Boundary Conditions: Real physical boundaries are not a major problem. ... “Open” or artificial boundaries occur whenever a domain is truncated. This ...
  10. [10]
    What Are the Navier-Stokes Equations? - COMSOL
    Jan 15, 2015 · The fluid velocity is specified at the inlet and pressure prescribed at the outlet. A no-slip boundary condition (i.e., the velocity is set to ...
  11. [11]
  12. [12]
  13. [13]
    [PDF] 3 The Navier-Stokes Equation - DAMTP
    The no-slip condition (3.11) doesn't follow from the Navier-Stokes equation. Instead, it is something additional that we assert. It is, however, physically ...
  14. [14]
    Linear stability analysis of Navier–Stokes–Voigt fluid with slip ...
    Jan 21, 2025 · Navier41 in 1823 proposed the model describing slip at a wall wherein he stated that the tangential fluid velocity directly adjacent to the ...
  15. [15]
    The no-slip condition of fluid dynamics | Erkenntnis
    Jun 9, 1989 · The historical perspective is based on the works of George Stokes and serves as the foundation for the philosophical perspective.
  16. [16]
    [PDF] Ludwig Prandtl and Boundary Layers in Fluid Flow
    By the no-slip condition, the fluid next to the lower plate has zero velo.city and the fluid next to the upper plate has velocity Uo; if the gap between the ...Missing: citation | Show results with:citation
  17. [17]
  18. [18]
    [PDF] Prandtl's Boundary Layer Theory - UC Davis Math
    The appropriate boundary condition for the Euler equations is the 'no-flow' conditions, u(x, t) . n(x) = 0 for x e dΩ, where n is the unit normal vector to ...Missing: citation | Show results with:citation
  19. [19]
  20. [20]
    [PDF] On the Wall Boundary Condition for Turbulence Models
    ABSTRACT. This report explains and discuss two main boundary condi- tions for turbulence models: the integrating and the wall fun- ction approach.
  21. [21]
    [PDF] The No-Slip Boundary Condition in Fluid Mechanics
    They derived the well- known equations of motion known now as the N avier-Stokes equations which relate the acceleration of the fluid element to the net force ...<|control11|><|separator|>
  22. [22]
    Inlet Boundary Condition - an overview | ScienceDirect Topics
    Inlet boundary conditions refer to a type of Dirichlet boundary condition where specific dependent variables such as velocity and turbulence quantities are ...
  23. [23]
    7.4. Boundary Conditions - ANSYS Help
    Velocity inlet boundary conditions are used to define the velocity and scalar properties of the flow at inlet boundaries. Pressure inlet boundary conditions are ...
  24. [24]
    Boundary Conditions in CFD - The Common Types
    An Inlet boundary condition specifies how the fluid enters the computational domain. ... Outlet boundary conditions specify how the fluid leaves the domain. The ...Missing: classification dynamics
  25. [25]
    Pressure Inlet and Outlet | Boundary Conditions - SimScale
    Aug 11, 2025 · The pressure inlet boundary condition defines an inflow condition based on the known pressure value (P) at the boundary.Pressure Inlet And Outlet... · 2. Fixed Value · TurbulenceMissing: wall | Show results with:wall<|control11|><|separator|>
  26. [26]
    Estimating the turbulence length scale - CFD Online
    Dec 7, 2024 · Setting a very low turbulence length scale at the inlet will quickly dissipate all turbulent energy and give a reduced turbulence intensity.
  27. [27]
    Inlet Values for the Turbulence Length Scale and Turbulent Intensity
    Fully turbulent flows usually have intensities between 0.05 (5%) and 0.1 (10%). The turbulence length scale LT is a measure of the size of the eddies that are ...Missing: specification | Show results with:specification
  28. [28]
    7.3 Inlet turbulence - Notes on CFD: General Principles
    Apr 11, 2022 · The values of eqn and eqn at domain inlets depend on the flow conditions upstream of the inlet. The figure below shows typical ranges of intensity eqn for ...
  29. [29]
    Exploring Inlet Boundary Conditions in CFD: Types and Applications
    Inlet boundary conditions in CFD determine how fluid enters the domain, crucial for accurate simulations. Options include velocity inlets for known speeds, ...
  30. [30]
    Assessing the impact of turbulent kinetic energy boundary ... - Nature
    Sep 5, 2023 · The results demonstrated that the sub-model simulations with an accurate boundary conditions successfully reconstruct the turbulent jet flow as ...Results · Methods · Computational Fluid Dynamics...<|separator|>
  31. [31]
    [PDF] Formulation and Implementation of Inflow/Outflow Boundary ...
    Boundary conditions for simulating flow entering or exiting a domain are implemented, using an algorithm to control mass flow rate through outflow and specify ...
  32. [32]
    Reconstruction of numerical inlet boundary conditions using ...
    Aug 24, 2021 · Therefore, the real challenge consists in defining a method capable of determining the appropriate mean and fluctuating fields that should be ...
  33. [33]
    Reconstruction of proper numerical inlet boundary conditions for ...
    Mar 30, 2023 · This paper discusses an innovative strategy to determine appropriate mean and fluctuating inlet boundary conditions for the numerical flow ...<|control11|><|separator|>
  34. [34]
    Boundary conditions for direct simulations of compressible viscous ...
    The emphasis is on deriving boundary conditions compatible with modern non-dissipative algorithms used for direct simulations of turbulent flows. These methods ...
  35. [35]
    ANSYS FLUENT 12.0 User's Guide - 7.3.11 Outflow Boundary ...
    The zero diffusion flux condition applied by ANSYS FLUENT at outflow boundaries is approached physically in fully-developed flows. Fully-developed flows are ...Missing: dynamics | Show results with:dynamics
  36. [36]
    Predictions and experimental studies of the tail pipe noise of an ...
    The tail pipe noise from a commercial automotive muffler was studied experimentally and numerically under the condition of wide open throttle acceleration.
  37. [37]
    Pressure Boundary Conditions | CFD-101 - FLOW-3D
    Pressure boundary conditions specify pressure at boundaries. There are two types: static, where pressure is continuous, and stagnation, where pressure drops to ...
  38. [38]
    [PDF] Pressure boundary condition for the time-dependent incompressible ...
    Sep 30, 2005 · On the pressure boundary conditions for the incompressible Navier–Stokes equations. International Journal for Numerical Methods in Fluids ...
  39. [39]
    [PDF] Pressure boundary conditions for computing incompressible flows ...
    We show that this scheme com- putes the pressure and velocity accurately near open boundaries and solid objects, and extends the scope of SPH simulation beyond ...
  40. [40]
    Turbulence statistics in fully developed channel flow at low ...
    Apr 21, 2006 · A direct numerical simulation of a turbulent channel flow is performed. The unsteady Navier-Stokes equations are solved numerically at a Reynolds number of ...Missing: conditions | Show results with:conditions
  41. [41]
    7.3.15 Symmetry Boundary Conditions - AFS ENEA
    Symmetry boundary conditions are used when the physical geometry of interest, and the expected pattern of the flow/thermal solution, have mirror symmetry.
  42. [42]
    Symmetric Boundary Condition - an overview | ScienceDirect Topics
    Symmetric boundary condition, a combination of Dirichlet and Neumann boundary conditions, is the condition that no fluid penetrates across the boundary.
  43. [43]
    Symmetry
    The Symmetry boundary condition prescribes no penetration and vanishing shear stresses. The boundary condition is a combination of a Dirichlet condition and a ...
  44. [44]
    [PDF] Application of Wray-Agarwal Turbulence Model for Accurate ...
    Apr 28, 2016 · A symmetry boundary condition is used for the symmetry plane of the wing/body. A no-slip wall boundary condition is used for the wing-body ...
  45. [45]
    [PDF] Computational Fluid Dynamics Analysis of the Stall Characteristics ...
    • Optimized high aspect ratio wing and high lift nacelle. • Tip cruise ... • Half-span with symmetry boundary condition for symmetric flow, full-span ...
  46. [46]
    Boundary conditions for CFD analysis | Download Table
    An aerofoil is a streamline body. Symmetric aerofoil (NACA 0015) is used in many applications such as in aircraft submarine fins, rotary and some fixed wings.
  47. [47]
  48. [48]
    Symmetry | Boundary Conditions - SimScale
    Feb 5, 2025 · The symmetry boundary condition defines a mirror face/surface. It should only be used if the physical object or geometry and the expected flow field pattern
  49. [49]
    [PDF] An Efficient Spectral-Projection Method for the Navier–Stokes ...
    An efficient and accurate numerical scheme is presented for the axisymmetric. Navier–Stokes equations in primitive variables in a cylinder.
  50. [50]
    Axisymmetric Flow Field - an overview | ScienceDirect Topics
    The boundary condition is the condition that the governing equation should satisfy on the boundary of fluid motion, typically including the inlet, outlet, and ...
  51. [51]
    [PDF] Dimensional/Axisymmetric - NASA Technical Reports Server
    Its boundary conditions for k and e at the wall are simple. • It gives reasonable turbulent and mean flow results near a solid wall. • It has been used by other ...
  52. [52]
    [PDF] Navier-Stokes Computational Study of Axisymmetric ... - DTIC
    Apr 16, 1985 · This study uses a Navier-Stokes code to compute turbulent flow over axisymmetric bodies at transonic speeds, using a k-epsilon two-equation ...
  53. [53]
    Axisymmetric flow of ideal fluid past a gas bubble | Fluid Dynamics
    A generalization of the analytic functions is proposed to solve the problem of axisymmetric flow past a gas bubble. By the use of an integral transformatio.
  54. [54]
    [PDF] An Introduction to Computational Fluid Dynamics
    It must be stressed that the continuative boundary condition has no physical basis, rather it is a mathematical statement that may or may not provide the ...
  55. [55]
    [PDF] Boundary Conditions - NASA
    The pressure values are assigned in the same way. The velocity components at the ghost cells are obtained as follows. ... The periodic boundary conditions are ...
  56. [56]
    How to implement boundary conditions in CFD
    Then, we reviewed the three fundamental types that exist: Dirichlet, Neumann, and Robin-type boundary conditions. We also looked at the two different ways ...
  57. [57]
    Direct Numerical Simulation of Flow and Heat Transfer in a ... - MDPI
    This study employs Direct Numerical Simulation (DNS) to investigate the flow and heat transfer characteristics in a compressor blade passage at five Reynolds ...<|separator|>
  58. [58]
    7.18. Periodic Condition - ANSYS Help
    Periodic boundary conditions are used when the physical geometry of interest and the expected flow pattern have a periodically repeating nature.
  59. [59]
    Rotationally periodic boundary conditions - CFD Online
    Mar 19, 2021 · Rotationally periodic boundary conditions are like translational ones, but with the angle of the face pair taken into account, changing the ...Rotational Periodic Boundary Conditions - CFD OnlinePeriodic Boundary Conditions -- CFD Online Discussion ForumsMore results from www.cfd-online.com
  60. [60]
    [PDF] UNHB894372.Online. 526..543 - AUB
    Consider the cyclic boundary condition illustrated in Figure 2. Two patches separated by some angle a are designated to form a periodic connection. For this.Missing: fluid dynamics
  61. [61]
    CFD Optimization of Cooling Air Offtake Passages Within Rotor ...
    In the tangential direction, a single ∼2.4 deg pitch of the 152 offtake passages was represented, thereby allowing cyclic boundary conditions to be imposed ...
  62. [62]
    Development of a Generalized Grid Interface for Turbomachinery ...
    PDF | This paper presents an implementation of the Generalized Grid Interface (GGI) for OpenFOAM. It is used to couple multiple non-conformal regions.
  63. [63]
    Far Field Boundary Conditions for Compressible Flows - ADS
    Far Field Boundary Conditions for Compressible Flows ... Numerical Flow Visualization;; Steady State;; Airfoils;; Computational Fluid Dynamics;; Computer ...
  64. [64]
    Far field boundary conditions for compressible flows - ScienceDirect
    Far field boundary conditions for compressible flows☆. Author links open ... Bayliss, E. Turkel. Commun. Pure Appl. Math., 33 (1980), p. 707. Crossref ...
  65. [65]
    Effect of Far-Field Boundary Conditions on Boundary-Layer Transition
    The effect of far-field boundary conditions on the evolution of a finite-amplitude two-dimensional wave in the Blasius boundary layer is assessed.Missing: dynamics | Show results with:dynamics
  66. [66]
    [2411.13077] Far-field Boundary Conditions for Airfoil Simulation at ...
    Nov 20, 2024 · This study concerns the far-field boundary conditions (BCs) for airfoil simulations at high incidence where the lift and drag are comparable in magnitude and ...
  67. [67]
    Air Pollution Dispersion Modelling in Urban Environment Using CFD
    This article presents a detailed review of the use of CFD to model air pollution dispersion in an urban environment over the last decade.
  68. [68]
    Analysis of sponge zones for computational fluid mechanics
    Sponge zones add a forcing term to equations, minimizing reflections and forcing solutions towards a target state, often used to introduce disturbances.
  69. [69]
    [PDF] Boundary and initial conditions
    Initial conditions specify velocity at t=0. Boundary conditions include rigid, free surface, inflow, and elastic boundaries. Free surfaces require kinematic ...
  70. [70]
    [PDF] Free surface flows in fluid dynamics - Thomas Alazard
    Feb 7, 2021 · The purpose of this chapter is to introduce the free surface Euler equation. We will begin by briefly explaining the meaning of these ...
  71. [71]
    On the initial values of kinematic and dynamic free-surface boundary ...
    Sep 3, 2020 · Abstract. The kinematic and dynamic boundary conditions on the free surface of a fluid should be posed for water wave problems.
  72. [72]
    [PDF] Free-Surface Waves - MIT
    Boundary Conditions. • In order to solve the boundary value problem for free surface waves we need to understand the boundary conditions on the free surface, ...
  73. [73]
    (PDF) Free-surface flows under impacting droplets - ResearchGate
    Aug 5, 2025 · A numerical method which fulfils the free-surface boundary conditions and extrapolates the fluid velocity into empty grid cells outside the ...<|control11|><|separator|>
  74. [74]
    Computations of boiling flows - ScienceDirect.com
    In this paper we present a two-dimensional front tracking/finite difference method for fluid flow with phase change that enables the simulation of problems with ...
  75. [75]
    Enhanced Level-Set Method for free surface flow applications - arXiv
    Nov 12, 2024 · This publication presents a solver using the Level-Set method ( Sussman et al. [1994] ) for incompressible two phase flows with surface tension.2 Mathematical And Numerical... · 2.3 Ghost Fluid Method For... · 3 Test Cases