# How to implement Stress-Free BC

**URL:** <https://community.freefem.org/t/how-to-implement-stress-free-bc/516>\
**Category:** General Discussion\
**Created:** [July 23, 2020, 7:44am UTC](https://community.freefem.org/t/how-to-implement-stress-free-bc/516 "2020-07-23T07:44:43Z")\
**Posts on this page:** 2\
**Page:** 1

<div class="post-metadata">

**Author:** ![Hanen](https://avatars.discourse-cdn.com/v4/letter/h/5f8ce5/32.png) [@Hanen](https://community.freefem.org/u/Hanen)\
**Post date:** [July 23, 2020, 7:44am UTC](https://community.freefem.org/t/how-to-implement-stress-free-bc/516/1 "2020-07-23T07:44:43Z")

</div>

Hello,  
I have a question about the implementation of stress free BC, actually i have a system of equations with 6 unknown (u1,u2…u6), i want to impose this BC on the u2 and u4.

dn(u1)=-u2/(4\*r0) with dn the normal derivative of u1  
and dn(u4)=2(r-r0)\*u4/(r \*r0)

Is it possible to implement this two conditions at the same time with FreeFem++ ?  
Any help,  
Thank you,

---

<div class="post-metadata">

**Author:** ![cmd](https://yyz2.discourse-cdn.com/flex030/user_avatar/community.freefem.org/cmd/32/67_2.png) [@cmd](https://community.freefem.org/u/cmd)\
**Post date:** [July 23, 2020, 12:06pm UTC](https://community.freefem.org/t/how-to-implement-stress-free-bc/516/2 "2020-07-23T12:06:23Z")

</div>

The short answer is yes, but the way you enforce such conditions depends on your weak formulation. Usually, you would try to apply the Divergence Theorem in deriving your equation such that a Neumann-type boundary integral term appears. Then homogeneous Neumann conditions will be naturally enforced by skipping this boundary integral in your variational form (hence the nickname “do-nothing” boundary condition). In your case, you will need to add an integral along the boundary to account for the inhomogeneous part.